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THE HIGH BONE MASS GENE OF 1 lql3.3 
5 RELATED APPLICATIONS 

This application is a continuation-in-part of Application No. 09/229,319, filed January 
13, 1999, which claims benefit of U.S. Provisional Application No. 60/071,449, filed January 

A 

13, 1998, and U.S. Provisional Application No. 60/1 05,5 11, filed October 23, 1998, all of 
which are herein incorporated by reference in their entirety. 

10 FIELD OF THE INVENTION 

The present invention relates generally to the field of genetics, genomics and 
molecular biology. More particularly, the invention relates to methods and materials used to 
isolate, detect and sequence a high bone mass gene and corresponding wild-type gene, and 
mutants thereof The present invention also relates to the high bone mass gene, the 

15 corresponding wild-type gene, and mutants thereof The genes identified in the present 

invention are implicated in the ontology and physiology of bone development. The invention 
also provides nucleic acids, proteins, cloning vectors, expression vectors, transformed hosts, 
methods of developing pharmaceutical compositions, methods of identifying molecules 
involved in bone development, and methods of diagnosing and treating diseases involved in 

20 bone development. In preferred embodiments, the present invention is directed to methods 
for treating, diagnosing, preventing and screening for normal and abnormal conditions of 
bone, including metabolic bone diseases such as osteoporosis. 

BACKGROUND OF THE INVENTION 
Two of the most common types of osteoporosis are postmenopausal and senile 

25 osteoporosis. Osteoporosis affects men as well as women, and, taken with other 

abnormalities of bone, presents an ever-increasing health risk for an aging population. The 
most common type of osteoporosis is that associated with menopause. Most women lose 
between 20-60% of the bone mass in the trabecular compartment of the bone within 3-6 years 
after the cessation of menses. This rapid loss is generally associated with an increase of bone 
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resorption and formation. However, the resorptive cycle is more dominant and the result is a 
net loss of bone mass. Osteoporosis is a common and serious disease among postmenopausal 
women. There are an estimated 25 million women in the United States alone who are 
afflicted with this disease. The results of osteoporosis are both personally harmful, and also 
5 account for a large economic loss due to its chronicity and the need for extensive and long- 
term support (hospitalization and nursing home care) from the disease sequelae. This is 
especially true in more elderly patients. Additionally, while osteoporosis is generally not 
thought of as a life-threatening condition, a 20-30% mortality rate is related to hip fractures in 
elderly women. A large percentage of this mortality rate can be directly associated with 

10 postmenopausal osteoporosis. 

The most vulnerable tissue in the bone to the effects of postmenopausal osteoporosis 
is the trabecular bone. This tissue is often referred to as spongy bone and is particularly 
concentrated near the ends of the bone near the joints and in the vertebrae of the spine. The 
trabecular tissue is characterized by small structures which inter-connect with each other as 

15 well as the more solid and dense cortical tissue which makes up the outer surface and central 
shaft of the bone. This criss-cross network of trabeculae gives lateral support to the outer 
cortical structure and is critical to the biomechanical strength of the overall structure. In 
postmenopausal osteoporosis, it is primarily the net resorption and loss of the trabeculae 
w hich lead to the failure and fracture of the bone. In light of the loss of the trabeculae in 

20 postmenopausal women, it is not surprising that the most common fractures are those 

associated with bones which are highly dependent on trabecular support, e.g., the vertebrae, 
the neck of the femur, and the forearm. Indeed, hip fracture, Colle's fractures, and vertebral 
crush fractures are indicative of postmenopausal osteoporosis. 

One of the earliest generally accepted methods for treatment of postmenopausal 

25 osteoporosis was estrogen replacement therapy. Although this therapy frequently is 

successful, patient compliance is low, primarily due to the undesirable side-effects of chronic 
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estrogen treatment. Frequently cited side-effects of estrogen replacement therapy include 
reinitiation of menses, bloating, depression, and fear of breast or uterine cancer. In order to 
limit the known threat of uterine cancer in those women who have not undergone a 
hysterectomy, a protocol of estrogen and progestin cyclic therapy is often employed. This 
5 protocol is similar to that which is used in birth control regimens, and often is not tolerated by 
women because of the side-effects characteristic of progestin. More recently, certain 
antiestrogens, originally developed for the treatment of breast cancer, have been shown in 
experimental models of postmenopausal osteoporosis to be efficacious. Among these agents 
is raloxifene (See, U.S. Patent No. 5,393,763, and Black et al, J. Clin. Invest., 93:63-69 

10 (1994)). In addition, tamoxifene, a widely used clinical agent for the treatment of breast 

cancer, has been shown to increase bone mineral density in post menopausal women suffering 
from breast cancer (Love et al, N. Engl. J. Med., 326:852-856 (1992)). 

Another therapy for the treatment of postmenopausal osteoporosis is the use of 
calcitonin. Calcitonin is a naturally occurring peptide which inhibits bone resorption and has 

15 been approved for this use in many countries (Overgaard et al, Br. Med. J., 305:556-561 
(1992)). The use of calcitonin has been somewhat limited, however. Its effects are very 
modest in increasing bone mineral density and the treatment is very expensive. Another 
therapy for the treatment of postmenopausal osteoporosis is the use of bis-phosphonates. 
These compounds were originally developed for use in Paget's disease and malignant 

20 hypercalcemia. They have been shown to inhibit bone resorption. Alendronate, one 

compound of this class, has been approved for the treatment of postmenopausal osteoporosis. 
These agents may be helpful in the treatment of osteoporosis, but these agents also have 
potential liabilities which include osteomalacia, extremely long half-life in bone (greater than 
2 years), and possible "frozen bone syndrome," e.g., the cessation of normal bone remodeling. 

25 Senile osteoporosis is similar to postmenopausal osteoporosis in that it is marked by 

the loss of bone mineral density and resulting increase in fracture rate, morbidity, and 
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associated mortality. Generally, it occurs in later life, i.e., after 70 years of age. Historically, 
senile osteoporosis has been more common in females, but with the advent of a more elderly 
male population, this disease is becoming a major factor in the health of both sexes. It is not 
clear what, if any, role hormones such as testosterone or estrogen have in this disease, and its 
5 etiology remains obscure. Treatment of this disease has not been very satisfactory. Hormone 
therapy, estrogen in women and testosterone in men, has shown equivocal results; calcitonin 
and bis-phosphonates may be of some utility. 

The peak mass of the skeleton at maturity is largely under genetic control. Twin 
studies have shown that the variance in bone mass between adult monozygotic twins is 

10 smaller than between dizygotic twins (Slemenda et al, J. Bone Miner. Res., 6:561-567 (1991); 
Young et al, J. Bone Miner. Res., 6:561-567 (1995); Pocock et al, J. Clin. Invest., 80:706-710 
(1987); Kelly et al, J. Bone Miner. Res., 8:1 1-17 (1993)), and it has been estimated that up to 
60% or more of the variance in skeletal mass is inherited (Krall et al, J. Bone Miner. Res., 
10:S367 (1993)). Peak skeletal mass is the most powerful determinant of bone mass in 

15 elderly years (Hui et al, Ann. Int. Med., 1 1 1 :355-361 (1989)), even though the rate of age- 
related bone loss in adult and later life is also a strong determinant (Hui et al, Osteoporosis 
Int., 1 :30-34 (1995)). Since bone mass is the principal measurable determinant of fracture 
risk, the inherited peak skeletal mass achieved at maturity is an important determinant of an 
individual's risk of fracture later in life. Thus, study of the genetic basis of bone mass is of 

20 considerable interest in the etiology of fractures due to osteoporosis. 

Recently, a strong interest in the genetic control of peak bone mass has developed in 
the field of osteoporosis. The interest has focused mainly on candidate genes with suitable 
polymorphisms to test for association with variation in bone mass within the normal range, or 
has focused on examination of genes and gene loci associated with low bone mass in the 

25 range found in patients with osteoporosis. The vitamin D receptor locus (VDR) (Morrison et 
ah Nature, 367:284-287 (1994)), PTH gene (Howard et al, J. Clin. Endocrinol. Metab., 
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80:2800-2805 (1995); Johnson ct al, J. Bone Miner. Res., 8:1 1-17 (1995); Gong ct al, J. Bone 
Miner. Res., 10:S462 (1995)) and the estrogen receptor gene (Hosoi et al, J. Bone Miner. 
Res., 1():S 170 (1995); Morrison et al, Nature, 367:284-287 (1994)) have figured most 
prominently in this work. These studies are difficult because bone mass (the phenotype) is a 

5 continuous, quantitative, polygenic trait, and is confounded by environmental factors such as 
nutrition, co-morbid disease, age, physical activity, and other factors. Also, this type of study 
design requires large numbers of subjects. In particular, the results of VDR studies to date 
have been confusing and contradictory (Garnero et al, J. Bone Miner. Res., 10:1283-1288 
(1995); Eisman et al,./ Bone. Miner. Res., 10:1289-1293 (1995); Peacock, J. Bone Miner. 

10 Res., 10:1294-1297 (1995)). Furthermore, the work thus far has not shed much light on the 
mechanism(s) whereby the genetic influences might exert their effect on bone mass. 

While it is well known that peak bone mass is largely determined by genetic rather 
than environmental factors, studies to determine the gene loci (and ultimately the genes) 
linked to variation in bone mass are difficult and expensive. Study designs which utilize the 

15 power of linkage analysis, e.g., sib-pair or extended family, are generally more informative 
than simple association studies, although the latter do have value. However, genetic linkage 
studies involving bone mass are hampered by two major problems. The first problem is the 
phenotype, as discussed briefly above. Bone mass is a continuous, quantitative trait, and 
establishing a discrete phenotype is difficult. Each anatomical site for measurement may be 

20 influenced by several genes, many of which may be different from site to site. The second 
problem is the age component of the phenotype. By the time an individual can be identified 
as having low bone mass, there is a high probability that their parents or other members of 
prior generations will be deceased and therefore unavailable for study, and younger 
generations may not have even reached peak bone mass, making their phenotyping uncertain 

25 for genetic analysis. 

Regardless, linkage analysis can be used to find the location of a gene causing a 
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hereditary "disorder" and does not require any know ledge of the biochemical nature of the 
disorder, i.e., a mutated protein that is believed to cause the disorder does not need to be 
known. Traditional approaches depend on assumptions concerning the disease process that 
might implicate a known protein as a candidate to be evaluated. The genetic localization 
5 approach using linkage analysis can be used to first find the general chromosomal region in 
which the defective gene is located and then to gradually reduce the size of the region in order 
to determine the location of the specific mutated gene as precisely as possible. After the gene 
itself is discovered w ithin the candidate region, the messenger RNA and the protein are 
identified and, along with the DNA, are checked for mutations. 
10 The genetic localization approach has practical implications since the location of the 

disease can be used for prenatal diagnosis even before the altered gene that causes the disease 
is found. Linkage analysis can enable families, even many of those that do not have a sick 
child, to know whether they are carriers of a disease gene and to evaluate the condition of an 
unborn child through molecular diagnosis. The transmission of a disease within families, 
15 then, can be used to find the defective gene. As used herein, reference to "high bone mass" 
(HBM) is analogous to reference to a disease state, although from a practical standpoint high 
bone mass can actually help a subject avoid the disease known as osteoporosis. 

Linkage analysis is possible because of the nature of inheritance of chromosomes 
from parents to offspring. During meiosis, the two parental homologues pair to guide their 
20 proper separation to daughter cells. While they are lined up and paired, the two homologues 
exchange pieces of the chromosomes, in an event called "crossing over" or "recombination." 
The resulting chromosomes are chimeric, that is, they contain parts that originate from both 
parental homologues. The closer together two sequences are on the chromosome, the less 
likely that a recombination event will occur between them, and the more closely linked they 
25 are. In a linkage analysis experiment, two positions on the chromosomes are followed from 
one generation to the next to determine the frequency of recombination between them. In a 
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study of an inherited disease, one of the chromosomal positions is marked by the disease gene 
or its normal counterpart, i.e., the inheritance of the chromosomal region can be determined 
by examining whether the individual displays symptoms of the disorder or not. The other 
position is marked by a DNA sequence that shows natural variation in the population such 
5 that the two homologies can be distinguished based on the copy of the "marker" sequence 
that they possess. In every family, the inheritance of the genetic marker sequence is 
compared to the inheritance of the disease state. If, w ithin a family carrying an autosomal 
dominant disorder such as high bone mass, every affected individual carries the same form of 
the marker and all the unaffected individuals carry at least one different form of the marker, 

10 there is a great probability that the disease gene and the marker are located close to each 
other. In this way, chromosomes may be systematically checked with known markers and 
compared to the disease state. The data obtained from the different families is combined, and 
analyzed together by a computer using statistical methods. The result is information 
indicating the probability of linkage between the genetic marker and the disease allowing 

15 different distances between them. A positive result can mean that the disease is very close to 
the marker, while a negative result indicates that it is far away on that chromosome, or on an 
entirely different chromosome. 

Linkage analysis is performed by typing all members of the affected family at a given 
marker locus and evaluating the co-inheritance of a particular disease state with the marker 

20 probe, thereby determining how often the two of them are co-inherited. The recombination 
frequency can be used as a measure of the genetic distance between two gene loci. A 
recombination frequency of 1% is equivalent to 1 map unit, or 1 centiMorgan (cM), which is 
roughly equivalent to 1,000 kb of DNA. This relationship holds up to frequencies of about 
20% or20cM. 

25 The entire human genome is 3,300 cM long. In order to find an unknown disease 

uene within 5-10 cM of a marker locus, the w hole human genome can be searched with 
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roughly 330 informative marker loci spaced at approximately 10 cM intervals (Botstein et al, 
Am. J. Hum. Genet., 32:314-331 (1980)). The reliability of linkage results is established by 
using a number of statistical methods. The method most commonly used for the analysis of 
linkage in humans is the LOD score method (Morton, Prog. Clin. Biol. Res., 147:245-265 
5 (1984), Morton et 'd\,Am. J. Hum. Genet., 38:868-883 (1986)) which was incorporated into 
the computer program LIPED by Ott, Am. J. Hum. Genet., 28:528-529 ( 1 976). LOD scores 
are the logarithm of the ratio of the likelihood that two loci are linked at a given distance to 
that they are not linked (>50 cM apart). The advantage of using logarithmic values is that 
they can be summed among families with the same disease. This becomes necessary given 
10 the relatively small size of human families. 

By convention, a total LOD score greater than + 3.0 (that is, odds of linkage at the 
specified recombination frequency being 1000 times greater than odds of no linkage) is 
considered to be significant evidence for linkage at that particular recombination frequency. 
A total LOD score of less than - 2.0 (that is, odds of no linkage being 100 times greater than 
15 odds of linkage at the specified frequency) is considered to be strong evidence that the two 
loci under consideration are not linked at that particular recombination frequency. Until 
recently, most linkage analyses have been performed on the basis of two-point data, which is 
the relationship between the disorder under consideration and a particular genetic marker. 
However, as a result of the rapid advances in mapping the human genome over the last few 
20 years, and concomitant improvements in computer methodology, it has become feasible to 
carry out linkage analyses using multi-point data. Multi-point analysis provide a 
simultaneous analysis of linkage between the disease and several linked genetic markers, 
when the recombination distance among the markers is known. 

Multi-point analysis is advantageous for two reasons. First, the informativeness of the 
25 pedigree is usually increased. Each pedigree has a certain amount of potential information, 
dependent on the number of parents heterozygous for the marker loci and the number of 
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affected individuals in the family. However, few markers are sufficiently polymorphic as to 
be informative in all those individuals. If multiple markers are considered simultaneously, 
then the probability of an individual being heterozygous for at least one of the markers is 
greatly increased. Second, an indication of the position of the disease gene among the 
5 markers may be determined. This allows identification of flanking markers, and thus 

eventually allow s isolation of a small region in which the disease gene resides. Lathrop et al, 
Proc. Nail. Acad. Sci. USA, 81:3443-3446 (1984) have written the most widely used 
computer package, LINKAGE, for multi-point analysis. 

There is a need in the art for identifying the gene associated with a high bone mass 
10 phenotype. The present invention is directed to this, as well as other, important ends. 

SUMMARY OF THF INVENTION 
The present invention describes the Zmaxl gene and the HBM gene on chromosome 

I lql3.3 by genetic linkage and mutation analysis. The use of additional genetic markers 
linked to the genes has aided this discovery. By using linkage analysis and mutation analysis, 

15 persons predisposed to HBM may be readily identified. Cloning methods using Bacterial 
Artificial Chromosomes have enabled the inventors to focus on the chromosome region of 

I I ql 3.3 and to accelerate the sequencing of the autosomal dominant gene. In addition, the 
invention identifies the Zmaxl gene and the HBM gene, and identifies the guanine-to- 
thymine polymorphism mutation at position 582 in the Zmaxl gene that produces the HBM 

20 gene and the HBM phenotype. 

The present invention identifies the Zmaxl gene and the HBM gene, which can be 
used to determine if people are predisposed to HBM and, therefore, not susceptible to 
diseases characterized by reduced bone density, including, for example, osteoporosis, or are 
predisposed and susceptible to diseases characterized by abnormally high bone density, such 
25 as. for example, osteoporosis. Older individuals carrying the HBM gene express the HBM 
protein, and, therefore, do not develop osteoporosis. In other words, the HBM gene is a 
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suppressor of osteoporosis. This /// vivo observation is a strong evidence that treatment of 
normal individuals with the HBM gene or protein, or fragments thereof, will ameliorate 

osteoporosis. 

Moreover, such treatment will be indicated in the treatment of bone lesions, 
particularly bone fractures, for bone remodeling in the healing of such lesions. For example, 
persons predisposed to or suffering from stress fractures (i.e., the accumulation of stress- 
induced microfractures, eventually resulting in a true fracture through the bone cortex) may 
be identified and/or treated by means of the invention. Moreover, the methods and 
compositions of the invention will be of use in the treatment of secondary osteoporosis, 
w here the course of therapy involves bone remodeling, such as endocrine conditions 
accompanying corticosteroid administration, hyperthyroidism, hypogonadism, hematologic 
malignancies, malabsorption and alcoholism, as well as disorders associated with vitamin D 
and/or phosphate metabolism, such as osteomalacia and rickets, and diseases characterized by 
abnormal or disordered bone remodeling, such as Paget's disease, and in neoplasms of bone, 
which may be benign or malignant. 

In various embodiments, the present invention is directed to nucleic acids, proteins, 
vectors, and transformed hosts of HBM and Zmaxl. 

Additionally, the present invention is directed to applications of the above 
embodiments of the invention including, for example, gene therapy, pharmaceutical 
development, and diagnostic assays for bone development disorders. In preferred 
embodiments, the present invention is directed to methods for treating, diagnosing, 
preventing and screening for osteoporosis. 

These and other aspects of the present invention are described in more detail below. 



BRIEF DESCRIPTION OF THE FIGURES 
Fig. 1 shows the pedigree of the individuals used in the genetic linkage studies. 
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Under each individual is an ID number, the /-score for spinal BMD, and the allele calls for 
the critical markers on chromosome 11. Solid symbols represent "affected" individuals. 
Symbols containing "N" are "unaffected" individuals. DNA from 37 individuals was 
genotyped. Question marks denote unknown genotypes or individuals who were not 
5 genotyped. 

Fig. 2 depicts the BAC/STS content physical map of the HBM region in 1 1 q 1 3.3. 
STS markers derived from genes, ESTs, microsatellites, random sequences, and BAC 
endsequences are denoted above the long horizontal line. For markers that are present in 
GDB the same nomenclature has been used. Locus names (Dl 1S####) are listed in 

10 parentheses after the primary name if available. STSs derived from BAC endsequences are 
listed with the BAC name first followed by L or R for the left and right end of the clone, 
respectively. The two large arrows indicate the genetic markers that define the HBM critical 
region. The horizontal lines below the STSs indicate BAC clones identified by PCR-based 
screening of a nine-fold coverage BAC library. Open circles indicate that the marker did not 

15 amplify the corresponding BAC library address during library screening. Clone names use 
the following convention: B for BAC, the plate, row and column address, followed by -H 
indicating the HBM project (i.e., B36F16-H). 

Figs. 3A-3F show the genomic structure of Zmaxl with flanking intron sequences. 
Translation is initiated by the underlined "ATG" in exon 1 . The site of the polymorphism in 

20 the HBM gene is in exon 3 and is represented by the underlined "G," whereby this nucleotide 
is a "T" in the HBM gene. The 3' untranslated region of the mRNA is underlined within exon 
23 (exon 1 , SEQ ID NO:40; exon 2, SEQ ID NO:41 ; exon 3, SEQ ID NO:42; exon 4, SEQ ID 
NO:43; exon 5, SEQ ID NO:44; exon 6, SEQ ID NO:45; exon 7, SEQ ID NO:46; exon 8, 
SEQ ID NO:47; exon 9, SEQ ID NO:48; exon 10, SEQ ID NO:49; exon 1 1 , SEQ ID NO:50; 

25 exon 12, SEQ ID NO:51; exon 13, SEQ ID NO:52; exon 14, SEQ ID NO:53; exon 15, SEQ 
ID NO:54; exon 16, SEQ ID NO:55; exon 17, SEQ ID NO:56; exon 18, SEQ ID NO:57; 
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exon 19, SEQ ID NO:58; exon 20, SEQ ID NO:59; exon 21, SEQ ID NO:60; exon 22, SEQ 
ID NO:61; and exon 23; SEQ ID NO:62). 

Fig. 4 shows the domain organization of Zmaxl, including the YWTD spacers, the 
extracellular attachment site, the binding site for LDL and calcium, the cysteine-rich growth 
5 factor repeats, the transmembrane region, the ideal PEST region w ith the CK-II 

phosphorylation site and the internalization domain. Fig. 4 also shows the site of the glycine 
to valine change that occurs in the HBM protein. The signal peptide is located at amino acids 
1-22, the extracellular domain is located at amino acids 23-1385, the transmembrane segment 
is located at amino acids 1386-1413, and the cytoplasmic domain is located at amino acids 
10 1414-1615. 

Fig. 5 is a schematic illustration of the BAC contigs B527D12 and B200E21 in 
relation to the HBM gene. 
^ r ^ Figs. 6A-6E are the nucleotide and amino acid sequences of the wild-type gene, 

Zmax 1 . The location for the base pair substitution at nucleotide 582, a guanine to thymine, is 
15 underlined. This allelic variant is the HBM gene. The HBM gene encodes for a protein with 
an ammo acid substitution of glycine to valine at position 171. The 5' untranslated region 
(UTR) boundaries bases 1 to 70, and the 3' UTR boundaries bases 4916-5120. 

Figs. 7A and 7B are northern blot analyses showing the expression of Zmaxl in 
various tissues. 
20 Fig. 8 is a PCR product analysis. 

Fig. 9 is allele specific oligonucleotide detection of the Zmaxl exon 3 mutation. 
Fig. 10 is the cellular localization of mouse Zmaxl by in situ hybridization at 100X 
magnification using sense and antisense probes. 

Fie. 11 is the cellular localization of mouse Zmaxl by in situ hybridization at 400X 
25 magnification using sense and antisense probes. 

; M r - Fig. 12 is the cellular localization of mouse Zmaxl by in situ hybridization of 
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osteoblasts in the endosteum at 400X magnification using sense and antisense probes. 
Fig. 13 shows antisense inhibition of Zmaxl expression in MC-3T3 cells. 



OFT All .ED DESCRIPTION OF THE IN VENTION 

S To aid in the understanding of the specification and claims, the following definitions 

are provided. 

"Gene" refers to a DNA sequence that encodes through its template or messenger 
RNA a sequence of amino acids characteristic of a specific peptide. The term "gene" includes 
intervening, non-coding regions, as well as regulatory regions, and can include 5* and 3' ends. 

10 "Gene sequence" refers to a DNA molecule, including both a DNA molecule which 

contains a non-transcribed or non-translated sequence. The term is also intended to include 
any combination of gene(s), gene fragment(s), non-transcribed sequence(s) or non-translated 
sequence(s) which are present on the same DNA molecule. 

The sequences of the present invention may be derived from a variety of sources 

15 including DNA, cDNA, synthetic DNA, synthetic RNA or combinations thereof. Such 

sequences may comprise genomic DNA which may or may not include naturally occurring 
introns. Moreover, such genomic DNA may be obtained in association with promoter regions 
or poly (A) sequences. The sequences, genomic DNA or cDNA may be obtained in any of 
several ways. Genomic DNA can be extracted and purified from suitable cells by means well 

20 known in the art. Alternatively, mRNA can be isolated from a cell and used to produce 
cDNA by reverse transcription or other means. 

"cDNA" refers to complementary or copy DNA produced from an RNA template by 
the action of RNA-dependent DNA polymerase (reverse transcriptase). Thus, a "cDNA 
clone" means a duplex DNA sequence complementary to an RNA molecule of interest, 
25 carried in a cloning vector or PCR amplified. This term includes genes from which the 
intervening sequences have been removed. 
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"Recombinant DNA" means a molecule that has been recombined by /// vitro splicing 
cDNA or a genomic DNA sequence. 

"Cloning" refers to the use of//; vitro recombination techniques to insert a particular 
gene or other DNA sequence into a vector molecule. In order to successfully clone a desired 
5 gene, it is necessary to use methods for generating DNA fragments, for joining the fragments 
to vector molecules, for introducing the composite DNA molecule into a host cell in which it 
can replicate, and for selecting the clone having the target gene from amongst the recipient 
host cells. 

"cDNA library" refers to a collection of recombinant DNA molecules containing 
10 cDNA inserts which together comprise the entire genome of an organism. Such a cDNA 
library can be prepared by methods known to one skilled in the art and described by, for 
example, Cowell and Austin, "cDNA Library Protocols," Methods in Molecular Biology 
(1997). Generally, RNA is first isolated from the cells of an organism from whose genome it 
is desired to clone a particular gene. 
15 "Cloning vehicle" refers to a plasmid or phage DNA or other DNA sequence which is 

able to replicate in a host cell. The cloning vehicle is characterized by one or more 
endonuclease recognition sites at which such DNA sequences may be cut in a determinable 
fashion without loss of an essential biological function of the DNA, which may contain a 
marker suitable for use in the identification of transformed cells. 
20 "Expression control sequence" refers to a sequence of nucleotides that control or 

regulate expression of structural genes when operably linked to those genes. These include, 
for example, the lac systems, the tip system, major operator and promoter regions of the 
phage lambda, the control region of fd coat protein and other sequences known to control the 
expression of genes in prokaryotic or eukaryotic cells. Expression control sequences will 
25 vary depending on whether the vector is designed to express the operably linked gene in a 
prokaryotic or eukaryotic host, and may contain transcriptional elements such as enhancer 
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elements, termination sequences, tissue-specificity elements and/or translational initiation and 
termination sites. 

"Expression vehicle" refers to a vehicle or vector similar to a cloning vehicle but 
which is capable of expressing a gene which has been cloned into it, after transformation into 
a host. The cloned gene is usually placed under the control of (i.e., operably linked to) an 
expression control sequence. 

"Operator" refers to a DNA sequence capable of interacting with the specific 
repressor, thereby controlling the transcription of adjacent gene(s). 

"Promoter" refers to a DNA sequence that can be recognized by an RNA polymerase. 
The presence of such a sequence permits the RNA polymerase to bind and initiate 
transcription of operably linked gene sequences. 

"Promoter region" is intended to include the promoter as well as other gene sequences 
which may be necessary for the initiation of transcription. The presence of a promoter region 
is sufficient to cause the expression of an operably linked gene sequence. 

"Operably linked" means that the promoter controls the initiation of expression of the 
gene. A promoter is operably linked to a sequence of proximal DNA if upon introduction 
into a host cell the promoter determines the transcription of the proximal DNA sequence(s) 
into one or more species of RNA. A promoter is operably linked to a DNA sequence if the 
promoter is capable of initiating transcription of that DNA sequence. 

"Prokaryote" refers to all organisms without a true nucleus, including bacteria. 
"Eukaryote" refers to organisms and cells that have a true nucleus, including 
mammalian cells. 

"Host" includes prokaryotes and eukaryotes, such as yeast and filamentous fungi, as 
well as plant and animal cells. The term includes an organism or cell that is the recipient of £ 
replicable expression vehicle. 

"Fragment" of a gene refers to any variant of the gene that possesses the biological 
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activ ity of that gene. 

"Variant" refers to a gene that is substantially similar in structure and biological 
activity or immunological characteristics to either the entire gene or to a fragment of the gene. 
Provided that the two genes possess a similar activity, they are considered variant as that 

5 term is used herein even if the sequence of amino acid residues is not identical. 

"Amplification of nucleic acids" refers to methods such as polymerase chain reaction 
(PCR), ligation amplification (or ligase chain reaction, LCR) and amplification methods 
based on the use of Q-beta replicase. These methods are well known in the art and described, 
for example, in U.S. Patent Nos. 4,683,195 and 4,683,202. Reagents and hardware for 

10 conducting PCR are commercially available. Primers useful for amplifying sequences from 
the HBM region are preferably complementary to, and hybridize specifically to sequences in 
the HBM region or in regions that flank a target region therein. HBM sequences generated by 
amplification may be sequenced directly. Alternatively, the amplified sequence(s) may be 
cloned prior to sequence analysis. 

15 "Antibodies" may refer to polyclonal and/or monoclonal antibodies and fragments 

thereof, and immunologic binding equivalents thereof, that can bind to the HBM proteins and 
fragments thereof or to nucleic acid sequences from the HBM region, particularly from the 
HBM locus or a portion thereof. The term antibody is used both to refer to a homogeneous 
molecular entity, or a mixture such as a serum product made up of a plurality of different 

20 molecular entities. Proteins may be prepared synthetically in a protein synthesizer and 

coupled to a carrier molecule and injected over several months into rabbits. Rabbit sera is 
tested for immunoreactivity to the HBM protein or fragment. Monoclonal antibodies may be 
made by injecting mice with the proteins, or fragments thereof. Monoclonal antibodies will 
be screened by ELISA and tested for specific immunoreactivity with HBM protein or 

25 fragments thereof. Harlow et al, Antibodies: A Laboratory Manual, Cold Spring Harbor 
Laboratory, Cold Spring Harbor, NY (1988). These antibodies will be useful in assays as 
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well as pharmaceuticals. 

"HBM" refers to high bone mass. 

"HBM protein" refers to a protein that is identical to a Zmaxl protein except that it 
contains an alteration of glycine 171 to valine. An HBM protein is defined for any organism 
that encodes a Zmaxl true homologue. For example, a mouse HBM protein refers to the 
mouse Zmaxl protein having the glycine 170 to valine substitution. 

"HBM gene" refers to the genomic DNA sequence found in individuals showing the 
HBM characteristic or phenotype, where the sequence encodes the protein indicated by SEQ 
ID NO: 4. The HBM gene and the Zmaxl gene are allelic. The protein encoded by the HBM 
gene has the property of causing elevated bone mass, while the protein encoded by the Zmaxl 
gene does not. The HBM gene and the Zmaxl gene differ in that the HBM gene has a 
thymine at position 582, while the Zmaxl gene has a guanine at position 582. The HBM 
gene comprises the nucleic acid sequence shown as SEQ ID NO: 2. The HBM gene may also 
be referred to as an "HBM polymorphism." 

"Normal," "wild-type," "unaffected" and "Zmaxl" all refer to the genomic DNA 
sequence that encodes the protein indicated by SEQ ID NO: 3. The Zmaxl gene has a 
guanine at position 582. The Zmaxl gene comprises the nucleic acid sequence shown as 
SEQ ID NO: 1. "Normal," "wild-type," "unaffected" and "Zmaxl" also refer to allelic 
variants of the genomic sequence that encodes proteins that do not contribute to elevated bone 
mass. The Zmaxl gene is common in the human population, while the HBM gene is rare. 

"5YWT+EGF" refers to a repeat unit found in the Zmaxl protein, consisting of five 
YWT repeats followed by an EGF repeat. 

"Bone development" generally refers to any process involved in the change of bone 
over time, including, for example, normal development, changes that occur during disease 
states, and changes that occur during aging. "Bone development disorder" particularly refers 
to any disorders in bone development including, for example, changes that occur during 
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disease states and changes that occur during aging. Bone development may be progressive or 
cyclical in nature. Aspects of bone that may change during development include, for 
example, mineralization, formation of specific anatomical features, and relative or absolute 
numbers of various cell types. 
5 "Bone modulation" or "modulation of bone formation" refers to the ability to affect 

any of the physiological processes involved in bone remodeling, as will be appreciated by one 
skilled in the art, including, for example, bone resorption and appositional bone growth, by, 
inter alia, osteoclastic and osteoblastic activity, and may comprise some or all of bone 
formation and development as used herein. 
10 "Normal bone density" refers to a bone density within two standard deviations of a Z 

score of 0. 

A "Zmaxl system" refers to a purified protein, cell extract, cell, animal, human or any 
other composition of matter in which Zmaxl is present in a normal or mutant form. 

A "surrogate marker" refers to a diagnostic indication, symptom, sign or other feature 
15 that can be observed in a cell, tissue, human or animal that is correlated with the HBM gene 
or elevated bone mass or both, but that is easier to measure than bone density. The general 
concept of a surrogate marker is well accepted in diagnostic medicine. 

The present invention encompasses the Zmaxl gene and Zmaxl protein in the forms 
indicated by SEQ ID NOS: 1 and 3, respectively, and other closely related variants, as well as 
20 the adjacent chromosomal regions of Zmaxl necessary for its accurate expression. In a 

preferred embodiment, the present invention is directed to at least 15 contiguous nucleotides 
of the nucleic acid sequence of SEQ ID NO: 1. 

The present invention also encompasses the HBM gene and HBM protein in the forms 
indicated by SEQ ID NO: 2 and 4, respectively, and other closely related variants, as well as 
25 the adjacent chromosomal regions of the HBM gene necessary for its accurate expression. In 
a preferred embodiment, the present invention is directed to at least 15 contiguous nucleotides 
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of the nucleic acid sequence of SEQ ID NO: 2. More preferably, the present invention is 
directed to at least 15 contiguous nucleotides of the nucleic acid sequence of SEQ ID NO: 2, 
wherein one of the 15 contiguous nucleotides is the thymine at nucleotide 582. 

The invention also relates to the nucleotide sequence of the Zmaxl gene region, as 
5 well as the nucleotide sequence of the HBM gene region. More particularly, a preferred 

embodiment are the BAC clones containing segments of the Zmaxl gene region B200E21-H 
and B527D12-H. A preferred embodiment is the nucleotide sequence of the BAC clones 
consisting of SEQ ID NOS: 5-12. 

The invention also concerns the use of the nucleotide sequence to identify DNA 
10 probes for the Zmaxl gene and the HBM gene, PCR primers to amplify the Zmaxl gene and 
the HBM gene, nucleotide polymorphisms in the Zmaxl gene and the HBM gene, and 
regulatory elements of the Zmaxl gene and the HBM gene. 

This invention describes the further localization of the chromosomal location of the 
Zmaxl gene and HBM gene on chromosome 1 lql3.3 between genetic markers Dl 1S987 and 
15 SNP CONTIG033-6, as well as the DNA sequences of the Zmaxl gene and the HBM gene. 
The chromosomal location was refined by the addition of more genetic markers to the 
mapping panel used to map the gene, and by the extension of the pedigree to include more 
individuals. The pedigree extension was critical because the new individuals that have been 
genotyped harbor critical recombination events that narrow the region. To identify genes in 
20 the region on 1 lq 13.3, a set of BAC clones containing this chromosomal region was 

identified. The BAC clones served as a template for genomic DNA sequencing, and also as a 
reagent for identifying coding sequences by direct cDNA selection. Genomic sequencing and 
direct cDNA selection were used to characterize more than 1 .5 million base pairs of DNA 
from 1 1 q 1 3.3. The Zmaxl gene was identified within this region and the HBM gene was 
25 then discovered after mutational analysis of affected and unaffected individuals. 

When a gene has been genetically localized to a specific chromosomal region, the 
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licnes in this region can be characterized at the molecular level by a series of steps that 
include: cloning of the entire region of DNA in a set of overlapping clones (physical 
mapping), characterization of genes encoded by these clones by a combination of direct 
cDNA selection, exon trapping and DNA sequencing (gene identification), and identification 
of mutations in these genes by comparative DNA sequencing of affected and unaffected 
members of the HBM kindred (mutation analysis). 

Physical mapping is accomplished by screening libraries of human DNA cloned in 
vectors that are propagated in E. coli or S. cereviseae using PCR assays designed to amplify 
unique molecular landmarks in the chromosomal region of interest. To generate a physical 
map of the HBM candidate region, a library of human DNA cloned in Bacterial Artificial 
Chromosomes (BACs) was screened with a set of Sequence Tagged Site (STS) markers that 
had been previously mapped to chromosome 1 Iql2-ql3 by the efforts of the Human Genome 
Project. 

STSs are unique molecular landmarks in the human genome that can be assayed by 
PCR. Through the combined efforts of the Human Genome Project, the location of thousands 
of STSs on the twenty-two autosomes and two sex chromosomes has been determined. For a 
positional cloning effort, the physical map is tied to the genetic map because the markers used 
for genetic mapping can also be used as STSs for physical mapping. By screening a BAC 
library with a combination of STSs derived from genetic markers, genes, and random DNA 
fragments, a physical map comprised of overlapping clones representing all of the DNA in a 
chromosomal region of interest can be assembled. 

BACs are cloning vectors for large (80 kilobase to 200 kilobase) segments of human 
or other DNA that are propagated in E. coli. To construct a physical map using BACs, a 
library of BAC clones is screened so that individual clones harboring the DNA sequence 
corresponding to a given STS or set of STSs are identified. Throughout most of the human 
genome, the STS markers are spaced approximately 20 to 50 kilobases apart, so that an 
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individual BAC clone typically contains at least two STS markers. In addition, the BAC 
libraries that were screened contain enough cloned DNA to cover the human genome six 
times over. Therefore, an individual STS typically identifies more than one BAC clone. By 
screening a six-fold coverage BAC library with a series of STS markers spaced 
5 approximately 50 kilobases apart, a physical map consisting of a series of overlapping BAC 
clones, i.e. BAC contigs, can be assembled for any region of the human genome. This map is 
closely tied to the genetic map because many of the STS markers used to prepare the physical 
map are also genetic markers. 

When constructing a physical map, it often happens that there are gaps in the STS 
10 map of the genome that result in the inability to identify BAC clones that are overlapping in a 
given location. Typically, the physical map is first constructed from a set of STSs that have 
been identified through the publicly available literature and World Wide Web resources. The 
initial map consists of several separate BAC contigs that are separated by gaps of unknown 
molecular distance. To identify BAC clones that fill these gaps, it is necessary to develop 
15 new STS markers from the ends of the clones on either side of the gap. This is done by 

sequencing the terminal 200 to 300 base pairs of the BACs flanking the gap, and developing a 
PCR assay to amplify a sequence of 100 or more base pairs. If the terminal sequences are 
demonstrated to be unique within the human genome, then the new STS can be used to screen 
the BAC library to identify additional BACs that contain the DNA from the gap in the 
20 physical map. To assemble a BAC contig that covers a region the size of the HBM candidate 
region (2,000,000 or more base pairs), it is often necessary to develop new STS markers from 
the ends of several clones. 

After building a BAC contig, this set of overlapping clones serves as a template for 
identifying the genes encoded in the chromosomal region. Gene identification can be 
25 accomplished by many methods. Three methods are commonly used: ( 1 ) a set of BACs 

selected from the BAC contig to represent the entire chromosomal region can be sequenced, 
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and computational methods can be used to identify all of the genes, (2) the BACs from the 
BAC contig can be used as a reagent to clone cDNAs corresponding to the genes encoded in 
the region by a method termed direct cDNA selection, or (3) the BACs from the BAC contig 
can be used to identify coding sequences by selecting for specific DNA sequence motifs in a 
5 procedure called exon trapping. The present invention includes genes identified by the first 
two methods. 

To sequence the entire BAC contig representing the HBM candidate region, a set of 
BACs was chosen for subcloning into plasmid vectors and subsequent DNA sequencing of 
these subclones. Since the DNA cloned in the BACs represents genomic DNA, this 

10 sequencing is referred to as genomic sequencing to distinguish it from cDNA sequencing. To 
initiate the genomic sequencing for a chromosomal region of interest, several non- 
overlapping BAC clones are chosen. DNA for each BAC clone is prepared, and the clones 
are sheared into random small fragments which are subsequently cloned into standard 
plasmid vectors such as pUC18. The plasmid clones are then grown to propagate the smaller 

15 fragments, and these are the templates for sequencing. To ensure adequate coverage and 
sequence quality for the BAC DNA sequence, sufficient plasmid clones are sequenced to 
yield six-fold coverage of the BAC clone. For example, if the BAC is 100 kilobases long, 
then phagemids are sequenced to yield 600 kilobases of sequence. Since the BAC DNA was 
randomly sheared prior to cloning in the phagemid vector, the 600 kilobases of raw DNA 

20 sequence can be assembled by computational methods into overlapping DNA sequences 
termed sequence contigs. For the purposes of initial gene identification by computational 
methods, six-fold coverage of each BAC is sufficient to yield ten to twenty sequence contigs 
of 1000 base pairs to 20,000 base pairs. 

The sequencing strategy employed in this invention was to initially sequence "seed" 
25 BACs from the BAC contig in the HBM candidate region. The sequence of the "seed" BACs 
was then used to identify minimally overlapping BACs from the contig, and these were 
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subsequently sequenced. In this manner, the entire candidate region was sequenced, with 
several small sequence gaps left in each BAC. This sequence served as the template for 
computational gene identification. One method for computational gene identification is to 
compare the sequence of BAC contig to publicly available databases of cDNA and genomic 
5 sequences, e.g. unigene, dbEST, genbank. These comparisons are typically done using the 
BLAST family of computer algorithms and programs (Altschul et al, J. Mol Biol, 215:403- 
410 (1990)). The BAC sequence can also be translated into protein sequence, and the protein 
sequence can be used to search publicly available protein databases, using a version of 
BLAST designed to analyze protein sequences (Altschul et al, Nucl Acids Res., 25:3389- 
10 3402 (1997)). Another method is to use computer algorithms such as MZEF (Zhang, Proc. 
Natl Acad. Srz., 94:565-568 (1997)) and GRAIL (Uberbacher et al, Methods Enzymol, 
266:259-281 (1996)), which predict the location of exons in the sequence based on the 
presence of specific DNA sequence motifs that are common to all exons, as well as the 
presence of codon usage typical of human protein encoding sequences. 
15 in addition to identifying genes by computational methods, genes were also identified 

by direct cDNA selection (Del Mastro et al, Genome Res. 5(2): 185- 194 (1995)). In direct 
cDN A selection, cDN A pools from tissues of interest are prepared, and the BACs from the 
candidate region are used in a liquid hybridization assay to capture the cDNAs which base 
pair to coding regions in the BAC. In the methods described herein, the cDNA pools were 
20 created from several different tissues by random priming the first strand cDNA from polyA 
RNA, synthesizing the second strand cDNA by standard methods, and adding linkers to the 
ends of the cDNA fragments. The linkers are used to amplify the cDNA pools. The BAC 
clones are used as a template for in vitro DNA synthesis to create a biotin labelled copy of the 
BAC DNA. The biotin labelled copy of the BAC DNA is then denatured and incubated with 
25 an excess of the PCR amplified, linkered cDNA pools which have also been denatured. The 
BAC DNA and cDNA are allowed to anneal in solution, and heteroduplexes between the 
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BAC and the cDNA are isolated using streptavidin coated magnetic beads. The cDNAs that 
are captured by the BAC are then amplified using primers complimentary to the linker 
sequences, and the hybridization/selection process is repeated for a second round. After two 
rounds of direct cDNA selection, the cDNA fragments are cloned, and a library of these direct 
5 selected fragments is created. 

The cDNA clones isolated by direct selection are analyzed by two methods. Since a 
pool of BACs from the HBM candidate region is used to provide the genomic DNA 
sequence, the cDNAs must be mapped to individual BACs. This is accomplished by arraying 
the BACs in microliter dishes, and replicating their DNA in high density grids. Individual 
10 cDNA clones are then hybridized to the grid to confirm that they have sequence identity to an 
individual BAC from the set used for direct selection, and to determine the specific identity of 
that BAC. cDNA clones that are confirmed to correspond to individual BACs are sequenced. 
To determine whether the cDNA clones isolated by direct selection share sequence identity 
or similarity to previously identified genes, the DNA and protein coding sequences are 
15 compared to publicly available databases using the BLAST family of programs. 

The combination of genomic DNA sequence and cDNA sequence provided by BAC 
sequencing and by direct cDNA selection yields an initial list of putative genes in the region. 
The genes in the region were all candidates for the HBM locus. To further characterize each 
gene, Northern blots were performed to determine the size of the transcript corresponding to 
20 each gene, and to determine which putative exons were transcribed together to make an 

individual gene. For Northern blot analysis of each gene, probes were prepared from direct 
selected cDNA clones or by PCR amplifying specific fragments from genomic DNA or from 
the BAC encoding the putative gene of interest. The Northern blots gave information on the 
size of the transcript and the tissues in which it was expressed. For transcripts which were 
25 not highly expressed, it was sometimes necessary to perform a reverse transcription PCR 
assay using RNA from the tissues of interest as a template for the reaction. 
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Gene identification by computational methods and by direct cDNA selection provides 
unique information about the genes in a region of a chromosome. When genes are identified, 
then it is possible to examine different individuals for mutations in each gene. 
I. Phenotyping using DXA Measurements 
5 Spinal bone mineral content (BMC) and bone mineral density (BMD) measurements 

performed at Creighton University (Omaha, Nebraska) were made by DXA using a Norland 
Instruments densitometer (Norland XR2600 Densitometer, Dual Energy X-ray 
Absorptiometry, DXA). Spinal BMC and BMD at other locations used the machinery 
available. There are estimated to be 800 DXA machines currently operating in the U.S. Most 
10 larger cities have offices or imaging centers which have DXA capabilities, usually a Lunar or 
Hologic machine. Each location that provided spine BMC and BMD data included copies of 
the printouts from their machines to provide verification that the regions of interest for 
measurement of BMD have been chosen appropriately. Complete clinical histories and 
skeletal radiographs were obtained. 
15 The HBM phenotype is defined by the following criteria: very high spinal BMD; a 

clinical history devoid of any known high bone mass syndrome; and skeletal radiographs 
showing a normal shape of the appendicular skeleton. 
II, Genotyping of Microsatellite Markers 

To narrow the genetic interval to a region smaller than that originally reported by 
20 Johnson et al. Am. J. Hum. Genet., 60:1326-1332 (1997), additional microsatellite markers on 
chromosome 1 1 q 1 2-1 3 were typed. The new markers included: Dl 1S4191, D11S1883, 
Dl 1S1785, D11S41 13, Dl 1S4136, Dl 1S4139, (Dib, et al, Nature, 380:152-154 (1996), 
FGF3 (Polymeropolous, et al, Nucl. Acid Res., 18:7468 (1990)), as well as 
GTC HBM Marker l , GTC HBM Marker 2, GTC HBM Marker 3, 
25 GTC HBM Marker 4, GTC HBM Marker 5, GTC^HBM_Marker_6, and 
GTC HBM Marker 7 (See Fig. 2). 
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Blood (20 ml) was drawn into lavender cap (EDTA containing) tubes by a certified 
phlebotomist. The blood was stored refrigerated until DNA extraction. DNA has been 
extracted from blood stored for up to 7 days in the refrigerator without reduction in the 
quality or quantity of yield. For those subjects that have blood drawn at distant sites, a 

5 shipping protocol was successfully used on more than a dozen occasions. Blood samples 
were shipped by overnight express in a styrofoam container with freezer packs to provide 
cooling. Lavender cap tubes were placed on individual plastic shipping tubes and then into 
"zip-lock" biohazard bags. When the samples arrived the next day, they were immediately 
processed to extract DNA. 

10 The DNA extraction procedure used a kit purchased from Gentra Systems, Inc. 

(Minneapolis, Minnesota). Briefly, the procedure involved adding 3 volumes of a red blood 
cell lysis buffer to the whole blood. After incubations for 10 minutes at room temperature, 
the solution was centrifuged in a Beckman tabletop centrifuge at 2,000 X g for 10 minutes. 
The white blood cell pellet was resuspended in Cell Lysis Buffer. Once the pellet was 

15 completely resuspended and free of cell clumps, the solution was digested with RNase A for 

GO 

(\ 1 5 minutes at 37 Proteins were precipitated by addition of the provided Protein 

Precipitation Solution and removed by centrifugation. The DNA was precipitated out of the 
supernatant by addition of isopropanol. This method was simple and fast, requiring only 1-2 
hours, and allowed for the processing of dozens of samples simultaneously. The yield of 
20 DNA was routinely >8 mg for a 20 ml sample of whole blood and had a MW of >50 kb. 
DNA was archived by storing coded 50 ug aliquots at -80^as an ethanol precipitate. 

DNA was genotyped using one fluorescently labeled oligonucleotide primer and one 
unlabeled oligonucleotide primer. Labeled and unlabeled oligonucleotides were obtained 
from Integrated DNA Technologies, Inc. (Coralville, Iowa). All other reagents for 
25 microsatellite genotyping were purchased from Perkin Elmer-Applied Biosystems, Inc. ("PE- 
ABI") (Norwalk, Connecticut). Individual PCR reactions were performed for each marker, as 
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described by PE-ABI using AmpliTag DNA Polymerase. The reactions were added to 3.5 ul 
ofloading buffer containing deionized formamide, blue dextran and TAMRA 350 size 
pi standards (PE-ABl). After heating at 95 jf for 5 minutes to denature the DNA, the samples 

were loaded and eleclrophoresed as described in the operator's manual for the Model 377 
5 DNA Sequencer (PE-ABl, Foster City, California). After gel electrophoresis, the data was 
analyzed using PE-ABl GENESCAN™ and GENOTYPER IM software. First, within the 
GENESCAN 1 M software, the lane tracking was manually optimized prior to the first step of 
analysis. After the gel lane data was extracted, the standard curve profiles of each lane were 
examined and verified for linearity and size calling. Lanes, which had problems with either 
10 of these parameters, were re-tracked and verified. Once all lanes were tracked and the size 
standards were correctly identified, the data were imported into GENOTYPER IM for allele 
identification To expedite allele calling (binning), the program Linkage Designer from the 
Internet web-site of Dr. Guy Van Camp (http://alt.www.uia.ac.be/u/dnalab/ld.html) was used. 
This program greatly facilitates the importing of data generated by GENOTYPER ,M into the 
15 pedigree drawing program Cyrillic (Version 2.0, Cherwell Scientific Publishing Limited, 
Oxford, Great Britain) and subsequent linkage analysis using the program LINKAGE 
(Lathrop et al, Am. J. Hum. Genet., 37:482-498 (1985)). 
III. Linkage Analysis 

Fig. 1 demonstrates the pedigree of the individuals used in the genetic linkage studies 
20 for this invention. Specifically, two-point linkage analysis was performed using the MLINK 
and LINKMAP components of the program LINKAGE (Lathrop et al, Am. J. Hum. Genet., 
37:482-498 (1985)). Pedigree/marker data was exported from Cyrillic as a pre-file into the 
Makeped program and converted into a suitable ped-file for linkage analysis. 

The original linkage analysis was performed using three models: (i) an autosomal 
25 dominant, fully penetrant model, (n) an autosomal dominant model with reduced penetrance, 
and (iii) a quantitative trait model. The HBM locus was mapped to chromosome 1 lql2-13 
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by analyzing DNA for linked markers from 22 members of a large, extended kindred. A 
highly automated technology was used with a panel of 345 fluorescent markers which 
spanned the 22 autosomes at a spacing interval ranging from 6-22 cM. Only markers from 
this region of chromosome 1 1 showed evidence of linkage (LOD score -3.0). The highest 
5 LOD score (5.74) obtained by two-point and multipoint analysis was Dl 1 S987 (map position 
55 in Fig. 2). The 95% confidence interval placed the HBM locus between markers Dl 1S905 
and D11S937 (map position 41-7 1 in Fig. 2). Haplotype analysis also places the Zmaxl gene 
in this same region. Further descriptions of the markers Dl 1S987, Dl 1S905, and Dl 1S937 
can be found in Gyapay et al, Nature Genetics, Vol. 7, (1994). 
10 In this invention, the inventors report the narrowing of the HBM interval to the region 

between markers Dl 1S987 and GTC HBM JVlarker_5 . These two markers lie between the 
delimiting markers from the original analysis (Dl 1S1 1S905 and Dl 1S937) and are 
approximately 3 cM from one another. The narrowing of the interval was accomplished 
using genotypic data from the markers Dl 1S4191, D11S1883, Dl 1 SI 785, Dl 1S41 13, 
15 Dl 1S4136, Dl 1S4139, (Dib et al, Nature, 380:152-154 (1996)), FGF3 (Polymeropolous et 
al, Nucl. Acid Res., 18:7468 (1990)) (information about the genetic markers can be found at 
the internet site of the Genome Database, http://gdbwww.gdb.org/), as well as the markers 
GTC_HBM_Marker_l , GTC_HBM_Marker_2, GTC HBM Marker 3, 
GTC_HBM Marker_4, GTC HBM JVlarker 5, GTC_HBM_Marker_6, and 
20 GTC HBM_Marker_7. 

As shown in Fig. 1, haplotype analysis with the above genetic markers identifies 
recombination events (crossovers) in individuals 9019 and 9020 that significantly refine the 
interval of chromosome 1 1 to which the Zmaxl gene is localized. Individual 9019 is an 
HBM-affected individual that inherits a portion of chromosome 1 1 from the maternal 
25 chromosome with the HBM gene, and a portion from the chromosome 1 1 homologue. The 
portion inherited from the HBM gene-carrying chromosome includes markers Dl 1S935, 
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Dl 1 SI 3 13, GTC_HBM_Marker_4, Dl 1S987, Dl 1S1296, GTC_HBM_Marker_6, 
GTC_HBM_Marker_2, Dl 1S970, GTC HBM Marker 3, Dl 1S41 13, 
GTC HBM Marker 1, GTC HBM Marker 7 and GTC HBM Marker 5. The portion 
from Dl 1S4136 and continuing in the telomeric direction is derived from the non-HBM 
5 chromosome. This data places the Zmaxl gene in a location centromeric to the marker 
GTC HBM Marker_5. Individual 9020 is an unaffected individual who also exhibits a 
critical recombination event. This individual inherits a recombinant paternal chromosome 1 1 
that includes markers Dl 1S935, Dl 1 SI 3 13, GTC HBM Marker_4, Dl 1S987, Dl IS 1296 
and GTC_HBM_Marker_6 from her father's (individual 01 15) chromosome 1 1 homologue 
10 that carries the HBM gene, and markers GTC HBM Marker_2, Dl 1S970, 
GTC HBM Marker 3, GTC_HBM_Marker_l , GTC_HBM_Marker_7, 
GTC HBM Marker 5, Dl 1S4136, Dl 1S4139, Dl 1S1314, and Dl 1S937 from her father's 
chromosome 1 1 that does not carry the HBM gene. Marker Dl 1 S41 1 3 is uninformative due 
to its homozygous nature in individual 0115. This recombination event places the 
15 centromeric boundary of the HBM region between markers Dl 1 SI 296 and Dl 1S987. 

Two-point linkage analysis was also used to confirm the location of the Zmaxl gene 
on chromosome 11. The linkage results for two point linkage analysis under a model of full 
penetrance are presented in Table 1 below. This table lists the genetic markers in the first 
column and the recombination fractions across the top of the table. Each cell of the column 
20 shows the LOD score for an individual marker tested for linkage to the Zmaxl gene at the 
recombination fraction shown in the first row. For example, the peak LOD score of 7.66 
occurs at marker Dl 1S970, which is within the interval defined by haplotype analysis. 

TABLE 1 



Marker 


0.0 


0.05 


0.1 


0.15 


0.2 


0.25 


0.3 


0.35 


0.4 
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D11S935 


- infinity 


0.39 


0.49 


0.47 


0.41 


0.33 


0.25 


0.17 


0.10 


Dl 1 SI 3 1 3 


- infinity 


2.64 


2.80 


2.80 


2.59 


2.30 


1 .93 


1.49 


1.00 


Dl 1S987 


- infinity 


5.49 


5.18 


4.70 


4.13 


3.49 


2.79 


2.03 


1.26 


Dl 1S41 13 


4.35 


3.99 


3.62 


3.24 


2.83 


2.40 


1.94 


1.46 


0.97 


Dl 1S1337 


2.29 


2.06 


1.81 


1.55 


1.27 


0.99 


0.70 


0.42 


0.18 


D11S970 


7.66 


6.99 


6.29 


5.56 


4.79 


3.99 
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4.52 
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3.51 
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A single nucleotide polymorphism (SNP) further defines the HBM region. This SNP 
is termed SNP_Contig033-6 and is located 25 kb centromeric to the genetic marker 
GTC_HBM Marker 5. This SNP is telomeric to the genetic marker GTC_HBM_Marker_7. 
SNP_Contig033-6 is present in HBM-affected individual 0113. However, the HBM-affected 
individual 9019, who is the son of 01 13, does not carry this SNP. Therefore, this indicates 
that the crossover is centromeric to this SNP. The primer sequence for the genetic markers 
GTC HBM_Marker_5 and GTC HBM Marker 7 is shown in Table 2 below. 

TABLE 2 



Marker 



Primer (Forward) 



Primer (Reverse) 
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GTCHBMMarker 

...5 


TTTTGGGTACACAATTCAGTC 
G 


AAAACTGTGGGTGCTTCT 
GG 


GTC HBM Marker 

7 


GTGATTGAGCCAATCCTGAGA 


TGAGCCAAATAAACCCCT 
TCT 



The kindred described have several features of great interest, the most important being 
that their bones, while very dense, have an absolutely normal shape. The outer dimensions of 
the skeletons of the HBM-affected individuals are normal, and, while medullary cavities are 
present, there is no interference with hematopoiesis. The HBM-affected members seem to be 
resistant to fracture, and there are no neurologic symptoms, and no symptoms of impairment 
of any organ or system function in the members examined. HBM-affected members of the 
kindred live to advanced age without undue illness or disability. Furthermore, the HBM 
phenotype matches no other bone disorders such as osteoporosis, osteoporosis pseudoglioma, 
Engelmann's disease, Ribbing's disease, hyperphosphatasemia, Van Buchem's disease, 
melorheostosis, osteopetrosis, pycnodysostosis, sclerostenosis, osteopoikilosis, acromegaly, 
Paget's disease, fibrous dysplasia, tubular stenosis, osteogenesis imperfecta, 
hypoparathyroidism, pseudohypoparathyroidism, pseudopseudohypoparathyroidism, primary 
and secondary hyperparathyroidism and associated syndromes, hypercalciuria, medullary 
carcinoma of the thyroid gland, osteomalacia and other diseases. Clearly, the HBM locus in 
this family has a very powerful and substantial role in regulating bone density, and its 
identification is an important step in understanding the pathway(s) that regulate bone density 
and the pathogenesis of diseases such as osteoporosis. 

In addition, older individuals carrying the HBM gene, and therefore expression of the 
HBM protein, do not show loss of bone mass characteristic of normal individuals. In other 
words, the HBM gene is a suppressor of osteoporosis. In essence, individuals carrying the 
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HBM gene are dosed with the HBM protein, and, as a result, do not develop osteoporosis. 
This //; vivo observation is strong evidence that treatment of normal individuals with the 
HBM gene or protein, or a fragment thereof, will ameliorate osteoporosis. 
IV. Physical Mapping 
5 To provide reagents for the cloning and characterization of the HBM locus, the 

genetic mapping data described above were used to construct a physical map of the region 
containing Zmaxl on chromosome 1 lql3.3. The physical map consists of an ordered set of 
molecular landmarks, and a set of BAC clones that contain the Zmaxl gene region from 
chromosome 1 1 q 1 3.3. 

10 Various publicly available mapping resources were utilized to identify existing STS 

markers (Olson et al, Science, 245:1434-1435 (1989)) in the HBM region. Resources 
included the GDB, the Whitehead Institute Genome Center, dbSTS and dbEST (NCBI), 
1 ldb, the University of Texas Southwestern GESTEC, the Stanford Human Genome Center, 
and several literature references (Courseaux et al, Genomics, 40:13-23 (1997), Courseaux et 

15 al, Genomics, 37:354-365 (1996), Guru et al, Genomics, 42:436-445 (1997), Hosoda et al, 

Genes Cells, 2:345-357 (1997), James et al, Nat. Genet., 8:70-76 (1994), Kitamura et al, DMA 
Research, 4:281-289 (1997), Lemmens et al, Genomics, 44:94-100 (1997), Smith et al, 
Genome Res., 7:835-842 (1997)). Maps were integrated manually to identify markers 
mapping to the region containing Zmaxl . 

20 Primers for existing STSs were obtained from the GDB or literature references are 

listed in Table 3 below. Thus, Table 3 shows the STS markers used to prepare the physical 
map of the Zmaxl gene region. 
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NDUFV1 1 


amplaxin ( EMS1 ) 


Nuclear mitotic apparatus protein 1 , NUMA 








LAR -interacting protein 1a mRNA 


Carnitine palmitoyi transferase I 


Beta-adrenergic receptor kinase 1 , ADRB1 




FGF4 




Human DNA helicase gen (SMBP2) 
















Aldehyde dehydrogenase (ALDH8) 


Human DNA helicase qen (SMBP2) 


Nuclear mitotic apparatus protein 1 , NUMA 


High sulphur keratin, KRN 




protein phosphatase 2A , PP2A 


Amplaxin 


Androgen Receptor 


51 C protein, inositol polyphosphate phosphatase-like 1 


N0F1 






NDUFV1 
















S13 avian erythroblastosis oncogene homolog 






Menin gene (MEN1) I 














Folate receptor 2 (FBP2) 


cGMP-stimulated 3', 5' -cyclic nucleotide phosphodiesterase PDE2A3 fPDE2A) 


Macrophage Migration Inhibitory factor 


P2U Purinoceptor | 








Folate receptor2 (FBP2) 




Folate receptor3 (FBP3) 






CATTCCACATGGATAGAC 


GCATTTTCTCATCATCCTTGC 


AGGTGTGTGTGCCAGGTTGA 


TTTG ATTTTGTGTC TCCC AAA 


CCAGTCACCTTTACTAGTCCTTTG I 


ACCAGGCATTGCACTAAAAG j 


ACATTTATATTTGGACATGCAACC ] 


ATGTGCTGGGCTGGAAAG | 


CTGCCTGTGTGGTGTCGC j 


GACCTCCTGTGACACCACG j 


GTAAGATTCTCCACTGTTGCACC | 


ACTCCTCATGTGAAGTCACCG 


CAGCATCTTCAGCACTTACC 


TGCTGCTGGGAGTCAGAGTC 


AAGGATCAAGCAGGCATTTG 


TGAACCCTGGAGGCAGAG 


GTGCTGGGATTACAGGTGT 


CCATGCTAGAGAAGCACAAC 


CAGACAGATAGCCCTGGGTTC 
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AGGATTCCTATCTGGGCTATG 


GAGAAGGCCGGGAGGCTCTG 
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Novel STSs were developed either from publicly available genomic sequence or from 
sequence-derived BAC insert ends. Primers were chosen using a script which automatically 
performs vector and repetitive sequence masking using Cross match (P. Green, U. of 
Washington) and subsequent primer picking using Primer3 (Ro/cn, Skaletsky (1996, 1997). 
5 Pnmer3 is available at www.genome.wi.mit. edu/genome_softwarc/other/primcr3.html. 

Polymerase chain reaction (PCR) conditions for each primer pair were initially 
optimized with respect to MgCU concentration. The standard buffer was 10 mM Tris-HCl 
(pH 8.3), 50 mM KC1, MgCU 0.2 mM each dNTP, 0.2 \xM each primer, 2.7 ng/|il human 
DNA, 0.25 units of AmpliTaq (Perkin Elmer) and MgCU concentrations of 1.0 mM, 1.5 mM, 
10 ft 2.0 mM or 2.4 mM. Cycling conditions included an initial denaturation at 94 j? for 2 minutes 
ft followed by 40 cycles at 94 g for 1 5 seconds, 55 ft for 25 seconds, and 12)C for 25 seconds 
jfa followed by a final extension at 72 £*for 3 minutes. Depending on the results from the initial 
round of optimization the conditions were further optimized if necessary. Variables included 
^increasing the annealing temperature to 58 £or 60 jZ^ increasing the cycle number to 42 and 
15 the annealing and extension times to 30 seconds, and using AmpliTaqGold (Perkin Elmer). 

BAC clones (Kim et al, Genomics, 32:213-218 (1996), Shizuya et al, Proc. Natl. 
Acad. Set. USA, 89:8794-8797 (1992)) containing STS markers of interest were obtained by 
PCR-based screening of DNA pools from a total human BAC library purchased from 
Research Genetics. DNA pools derived from library plates 1-596 were used corresponding to 
20 nine genomic equivalents of human DNA. The initial screening process involved PCR 

reactions of individual markers against superpools, i.e., a mixture of DNA derived from all 
BAC clones from eight 384-well library plates. For each positive superpool, plate (8), row 
(16) and column (24) pools were screened to identify a unique library address. PCR products 
were electrophoresed in 2% agarose gels (Sigma) containing 0.5 pg/ml ethidium bromide in 
25 IX TBE at 150 volts for 45 min. The electrophoresis units used were the Model A3-1 

systems from Ow l Scientific Products. Typically, gels contained 10 tiers of lanes with 50 
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wells tier. Molecular weight markers ( 100 bp ladder, Life Technologies, Bethesda, MD) 
w ere loaded at both ends of the gel. Images of the gels were captured with a Kodak DC40 
CCD camera and processed with Kodak 1 D software. The gel data were exported as tab 
delimited text files; names of the files included information about the library screened, the gel 
5 image files and the marker screened. These data were automatically imported using a 
customized Perl script into Filemaker PRO (Claris Corp.) databases for data storage and 
analysis. In cases where incomplete or ambiguous clone address information was obtained, 
additional experiments were performed to recover a unique, complete library address. 

Recovery of clonal BAC cultures from the library involved streaking out a sample 
10 from the library well onto LB agar (Maniatis et al, Molecular Cloning: A Laboratory 

Manual., Cold Spring Harbor Laboratory, Cold Spring Harbor, NY (1982)) containing 12.5 
jig/ml chloramphenicol (Sigma). Two individual colonies and a portion of the initial streak 
quadrant were tested with appropriate STS markers by colony PCR for verification. Positive 
clones were stored in LB broth containing 12.5 jig/ml chloramphenicol and 15% glycerol at - 

Several different types of DNA preparation methods were used for isolation of BAC 
DNA. The manual alkaline lysis miniprep protocol listed below (Maniatis et al, Molecular 
Cloning: A Laboratory Manual, Cold Spring Harbor Laboratory, Cold Spring Harbor, NY 
(1982)) was successfully used for most applications, i.e., restriction mapping, CHEF gel 
20 analysis, FISH mapping, but was not successfully reproducible in endsequencing. The 
Autogen and Qiagen protocols were used specifically for BAC DNA preparation for 
endsequencing purposes. 

Bacteria were grown in 15 ml Terrific Broth containing 12.5 ng/ml chloramphenicol 
(V in a 50 ml conical tube at 37 fc for 20 hrs with shaking at 300 rpm. The cultures were 
25 (\ centrifuged in a Sorvall RT 6000 D at 3000 rpm (-1800 g) at 4^for 15 min. The supernatant 
fa w as then aspirated as completely as possible. In some cases cell pellets w ere frozen at -20 
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ai this step tor up to 2 weeks. The pellet was then vortexed to homogenize the cells and 
minimize clumping. 250 (il of PI solution (50 mM glucose, 15 mM Tris-HCl, pH 8, 10 mM 
EDTA, and 100 ng/ml RNase A) was added and the mixture pipetted up and down to mix. 
The mixture was then transferred to a 2 ml Eppendorf tube. 350 jil of P2 solution (0.2 N 
5 NaOEf 1% SDS) was then added, the mixture mixed gently and incubated for 5 min. at room 
temperature. 350 jul of P3 solution (3M KOAc, pH 5.5) was added and the mixture mixed 
gently until a white precipitate formed. The solution was incubated on ice for 5 min. and then 
ft centrifuged at Ay in a microfuge for 10 min. The supernatant was transferred carefully 

(av oiding the white precipitate) to a fresh 2 ml Eppendorf tube, and 0.9 ml of isopropanol was 

10 added, the solution mixed and left on ice for 5 min. The samples were centrifuged for 10 
mm., and the supernatant removed carefully. Pellets were washed in 70% ethanol and air 
dried for 5 min. Pellets were resuspended in 200 \x\ of TE8 (10 mM Tris-HCl, pH 8.0, 1.0 
mM EDTA), and RNase A (Boehnnger Mannheim) added to 100 ng/ml. Samples were 
incubated at 37^2Tor 30 min., then precipitated by addition of C 2 H 3 0 2 Na 3H 2 0 to 0.5 M and 2 

15 volumes of ethanol. Samples were centrifuged for 10 min., and the pellets washed with 70% 
ethanol followed by air drying and dissolving in 50 |nl TE8. Typical yields for this DNA prep 
were 3-5 ng/15 ml bacterial culture. Ten to 15 ^1 were used for Hindlll restriction analysis; 5 
|il was used for NotI digestion and clone insert sizing by CHEF gel electrophoresis. 
BACs w ere inoculated into 15 ml of 2X LB Broth containing 12.5 |ig/ml 

20 chloramphenicol in a 50 ml conical tube. 4 tubes were inoculated for each clone. Cultures 
were grown overnight (-16 hr) at 37^with vigorous shaking (>300 rpm). Standard 
conditions for BAC DNA isolation were followed as recommended by the Autogen 740 
manufacturer. 3 ml samples of culture were placed into Autogen tubes for a total of 60 ml or 
20 tubes per clone. Samples were dissolved finally in 100 |il TE8 with 15 seconds of shaking 

25 as part of the Autogen protocol. After the Autogen protocol was finished DNA solutions 
w ere transferred from each individual tube and pooled into a 2 ml Eppendorf tube. Tubes 
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w ith large amounts of debris (carry over from the pelleting debris step) were avoided. The 
tubes were then rinsed with 0.5 ml of TE8 successively and this solution added to the pooled 
material. DNA solutions were stored at 4 J?; clumping tended to occur upon freezing at -20 {?. 
This DNA was either used directly for restriction mapping, CHEF gel analysis or FISH 
5 mapping or was further purified as described below for use in endsequencing reactions. 

The volume of DNA solutions was adjusted to 2 ml with TE8, samples were then 
A mixed gently and heated at 65 ^for 10 min. The DNA solutions were then centrifuged at 4^ 
for 5 min. and the supernatants transferred to a 15 ml conical tube. The NaCl concentration 
was then adjusted to 0.75 M (-0.3 ml of 5 M NaCl to the 2 ml sample). The total volume 
10 was then adjusted to 6 ml with Qiagen column equilibration buffer (Buffer QBT). The 
supernatant containing the DNA was then applied to the column and allowed to enter by 
gravity flow. Columns were washed twice with 10 ml of Qiagen Buffer QC. Bound DNA 
ft was then eluted with four separate 1 ml aliquots of Buffer QF kept at 65^?. DNA was 

precipitated with 0.7 volumes of isopropanol (-2.8 ml). Each sample was then transferred to 
15 4 individual 2.2 ml Eppendorf tubes and incubated at room temperature for 2 hr or overnight. 
(\ Samples were centrifuged in a microfuge for 10 min. at 4 ^ The supernatant was removed 
carefully and 1 ml of 70% ethanol was added. Samples were centrifuged again and because 
the DNA pellets were often loose at this stage, the supernatant removed carefully. Samples 
were centrifuged again to concentrate remaining liquid which was removed with a micropipet 
20 tip. DNA pellets were then dried in a desiccator for 10 min. 20 |il of sterile distilled and 
Q deionized H.O was added to each tube which was then placed at 4 ^overnight. The four 20 
ii 1 samples for each clone were pooled and the tubes rinsed with another 20 jal of sterile 

y 

(\ distilled and deionized H : 0 for a final volume of 100 jal. Samples were then heated at 65 J? 
for 5 min. and then mixed gently. Typical yields were 2-5 ng/60 ml culture as assessed by 
25 NotI digestion and comparison with uncut lambda DNA. 

3 ml of LB Broth containing 12.5 ng/ml of chloramphenicol was dispensed into 
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autoclaved Autogen tubes. A single tube was used for each clone. For inoculation, glycerol 
f\ stocks were removed from -70 storage and placed on dry ice. A small portion of the 
glycerol stock was removed from the original tube with a sterile toothpick and transferred 
into the Autogen tube; the toothpick was left in the Autogen tube for at least two minutes 
5 before discarding. After inoculation the tubes were covered with tape making sure the seal 
was tight. When all samples were inoculated, the tube units were transferred into an Autogen 
v ^ rack holder and placed into a rotary shaker at 37^for 16-17 hours at 250 rpm. Following 
growth, standard conditions for BAC DNA preparation, as defined by the manufacturer, were 
used to program the Autogen. Samples were not dissolved in TE8 as part of the program and 
10 DNA pellets were left dry. When the program was complete, the tubes were removed from 
the output tray and 30 |il of sterile distilled and deionized H 2 0 was added directly to the 
bottom of the tube. The tubes were then gently shaken for 2-5 seconds and then covered with 
parafilm and incubated at room temperature for 1-3 hours. DNA samples were then 
P( transferred to an Eppendorf tube and used either directly for sequencing or stored at 4 / (2 / for 
15 later use. 

V. BAC Clone Characterization for Physical Mapping 

DNA samples prepared either by manual alkaline lysis or the Autogen protocol were 
digested w ith Hindlll for analysis of restriction fragment sizes. This data were used to 
compare the extent of overlap among clones. Typically 1-2 jug were used for each reaction. 
20 Reaction mixtures included: IX Buffer 2 (New England Biolabs), 0.1 mg/ml bovine serum 

albumin (New England Biolabs), 50 |ig/ml RNase A (Boehringer Mannheim), and 20 units of^ 
(\ Hindlll (New England Biolabs) in a final volume of 25 |ul. Digestions were incubated at 37/^ 
for 4-6 hours. BAC DNA was also digested with NotI for estimation of insert size by CHEF 
gel analysis (see below). Reaction conditions were identical to those for Hindlll except that 
25 20 units of Notl were used. Six \x\ of 6X Ficoll loading buffer containing bromphenol blue 
and xylene cyanol was added prior to electrophoresis. 
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Hindi 1 1 digests were analyzed on 0.6% agarose (Seakem, FMC Bioproducts) in IX 
TBE containing 0.5 ng/ml ethidium bromide. Gels (20 cm X 25 cm) were electrophoresed in 
a Model A4 electrophoresis unit (Owl Scientific) at 50 volts for 20-24 hrs. Molecular weight 
size markers included undigested lambda DNA, Hindlll digested lambda DNA, and Haelll 
5 C\ diuested XI 74 DNA. Molecular weight markers were heated at 65 , 2 min. prior to 
loading the gel. Images were captured with a Kodak DC40 CCD camera and analyzed with 
Kodak ID software. 

NotI digests were analyzed on a CHEF DRJI (BioRad) electrophoresis unit according 
to the manufacturer's recommendations. Briefly, 1% agarose gels (BioRad pulsed field 
10 grade) were prepared in 0.5X TBE, equilibrated for 30 minutes in the electrophoresis unit at 
^ 14 jc, and electrophoresed at 6 volts/cm for 14 hrs with circulation. Switching times were 
ramped from 10 sec to 20 sec. Gels were stained after electrophoresis in 0.5 jig/ml ethidium 
bromide. Molecular weight markers included undigested lambda DNA, Hindlll digested 
lambda DNA, lambda ladder PFG ladder, and low range PFG marker (all from New England 
15 Biolabs). 

BAC DNA prepared either by the manual alkaline lysis or Autogen protocols were 
labeled for FISH analysis using a Bioprime labeling kit (BioRad) according to the 
manufacturer's recommendation with minor modifications. Approximately 200 ng of DNA 
was used for each 50 |al reaction. 3 were analyzed on a 2% agarose gel to determine the 
20 extent of labeling. Reactions were purified using a Sephadex G50 spin column prior to in situ 
hybridization. Metaphase FISH was performed as described (Ma et al, Cytogenet. Cell 
Genet., 74:266-271 (1996)). 
VI. BAC Endsequencing 

The sequencing of BAC insert ends utilized DNA prepared by either of the two 
25 methods described above. The DYEnamic energy transfer primers and Dynamic Direct cycle 
sequencing kits from Amersham were used for sequencing reactions. Ready made 

-43- 



Atty. Dkt. No.: 032796-014 
sequencing mix including the M13 -40 forward sequencing primer was used (Catalog # 
US79730) for the T7 BAC vector terminus; ready made sequencing mix (Catalog # 
US79530) was mixed with the M13 -28 reverse sequencing primer (Catalog # US79339) for 
the SP6 BAC vector terminus. The sequencing reaction mixes included one of the four 
5 fluorcscently labeled dye-primers, one of the four dideoxy termination mixes, dNTPs, 

reaction buffer, and Thermosequenase. For each BAC DNA sample, 3 jil of the BAC DNA 
sample was aliquoted to 4 PCR strip tubes. 2 jil of one of the four dye primer/termination 
mix combinations was then added to each of the four tubes. The tubes were then sealed and 
centrifuged briefly prior to PCR. Thermocychng conditions involved a 1 minute denaturation 

/* i. * 

10 # at 95 ^ 15 second annealing at 45 and extension for 1 minute at 70^'for 35 total cycles. 
After cycling the plates were centrifuged briefly to collect all the liquid to the bottom of the 
tubes. 5 \i\ of sterile distilled and deionized H.O was then added into each tube, the plates 
sealed and centrifuged briefly again. The four samples for each BAC were then pooled 
p^- together. DNA was then precipitated by adding 1.5 \x\ of 7.5 M NH 4 OAc and 100 jil 
15 100% ethanol to each tube. Samples were mixed by pipetting up and down once. The plates 
were then sealed and incubated on ice for 10 minutes. Plates were centrifuged in a table top 
Haraeus centrifuge at 4000 rpm (3,290 g) for 30 minutes at 4$fto recover the DNA. The 
supernatant was removed and excess liquid blotted onto paper towels. Pellets were washed 
^ by adding 100 \x\ of -20£ 70% ethanol into each tube and recentrifuging at 4000 rpm (3,290 
20 f\ g) for 10 minutes at 4^ The supernatant was removed and excess liquid again removed by 
blotting on a paper towel. Remaining traces of liquid were removed by placing the plates 
upside dow n over a paper tow el and centrifuging only until the centrifuge reached 800 rpm. 
Samples were then air dried at room temperature for 30 min. Tubes were capped and stored 
ft dry at -20 Cruntil electrophoresis. Immediately prior to electrophoresis the DNA was 
25 dissolved in 1 .5 jiil of Amersham loading dye. Plates were then sealed and centrifuged at 

2000 rpm (825 g). The plates were then vortexed on a plate shaker for 1-2 minutes. Samples 

-44- 




Atty. Dkt. No.: 032796-014 

f\ were then recentrifuged at 2000 rpm (825 g) briefly. Samples were then heated at 65 C for 2 
min. and immediately placed on ice. Standard gel electrophoresis was performed on AB1 377 
fluorescent sequencers according to the manufacturer's recommendation. 
VII. Sub-cloning and Sequencing of HBM BAC DNA 
5 The physical map of the Zmaxl gene region provides a set of BAC clones that contain 

within them the Zmaxl gene and the HBM gene. DNA sequencing of several of the BACs 
from the region has been completed. The DNA sequence data is a unique reagent that 
includes data that one skilled in the art can use to identify the Zmaxl gene and the HBM 
gene, or to prepare probes to identify the gene(s), or to identify DNA sequence 
10 polymorphisms that identify the gene(s). 

BAC DNA was isolated according to one of two protocols, either a Qiagen 
purification of BAC DNA (Qiagen, Inc. as described in the product literature) or a manual 
purification which is a modification of the standard alkaline lysis/Cesium Chloride 
preparation of plasmid DNA (see e.g., Ausubel et al, Current Protocols in Molecular Biology, 
15 John Wiley & Sons (1997)). Briefly for the manual protocol, cells were pelleted, 

resuspended in GTE (50 mM glucose, 25 mM Tris-Cl (pH 8), 10 mM EDTA) and lysozyme 
ft (50 mg/ml solution), followed by NaOH/SDS (l°/c^S©SA3N NaOH) and then an ice-cold 
solution of 3M KOAc (pH 4.5-4.8). RnaseA was added to the filtered supernatant, followed 
by Proteinase K and 20% SDS. The DNA was then precipitated with isopropanol, dried and 
20 resuspended in TE (10 mM Tris, 1 mM EDTA (pH 8.0)). The BAC DNA was further 
purified by Cesium Chloride density gradient centrifugation (Ausubel et al, Current 
Protocols in Molecular Biology, John Wiley & Sons (1997)). 

Following isolation, the BAC DNA was sheared hydrodynamically using an HPLC 
(Hengen, Trends in Biochem. Sci., 22:273-27 r 4 (1997)) to an insert size of 2000-3000 bp. 
25 After shearing, the DNA was concentrated and separated on a standard 1% agarose gel. A 

single fraction, corresponding to the approximate size, was excised from the gel and purified 
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bv elcctroelulion (Sambrook et al, Molecular Cloning: A Laboratory Manual, Cold Spring 



Harbor Laboratory, Cold Spring, NY (1989)). 

The purified DNA fragments were then blunt-ended using T4 DNA polymerase. The 
blunt-ended DNA was then ligated to unique BstXI-linker adapters (5' 
5 GTCTTCACCACGGGG and 5' GTGGTGAAGAC in 100-1000 fold molar excess). These 
linkers were complimentary to the BstXI-cut pMPX vectors (constructed by the inventors), 
while the overhang was not self-complimentary. Therefore, the linkers would not 
concatemerize nor would the cut-vector religate itself easily. The linker-adapted inserts were 
separated from the unincorporated linkers on a 1% agarose gel and purified using GeneClean 
10 (BIO 101, Inc.). The linker-adapted insert was then ligated to a modified pBlueScript vector 
to construct a "shotgun" subclone library. The vector contained an out-of-frame lacZ gene at 
the cloning site which became in-frame in the event that an adapter-dimer is cloned, allowing 
these to be avoided by their blue-color. 



15 sequencing methods. Only major modifications to the protocols are highlighted. Briefly, the 
library was then transformed into DH5 competent cells (Life Technologies, Bethesda, MD, 
DH5 transformation protocol). It was assessed by plating onto antibiotic plates containing 



f\ ampicillin and IPTG/Xgal. The plates were incubated overnight at 37£. Successful 

transformants were then used for plating of clones and picking for sequencing. The cultures 



20 A- were grown overnight at 37 f. DNA was purified using a silica bead DNA preparation (Ng et 
al, Nuci Acids Res,, 24:5045-5047 (1996)) method. In this manner, 25 |ig of DNA was 
obtained per clone. 

These purified DNA samples were then sequenced using ABI dye-terminator 
chemistry. The ABI dye terminator sequence reads were run on ABI377 machines and the 
25 data was directly transferred to UNIX machines following lane tracking of the gels. All reads 
were assembled using PHRAP (P. Green, Abstracts of DOE Human Genome Program 



All subsequent steps were based on sequencing by ABI377 automated DNA 
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Contractor-Grantee Workshop V, Jan. 1996, p. 157) with default parameters and quality 
scores. The initial assembly was done at 6-fold coverage and yielded an average of 8-1 5 
contigs. Following the initial assembly, missing mates (sequences from clones that only 
gave one strand reads) were identified and sequenced with ABI technology to allow the 
identification of additional overlapping contigs. Primers for walking were selected using a 
Genome Therapeutics program Pick primer near the ends of the clones to facilitate gap 
closure. These walks were sequenced using the selected clones and primers. Data were 
reassembled with PHRAP into sequence contigs. 
VIII. Gene Identification by Computational Methods 

Following assembly of the BAC sequences into contigs, the contigs were subjected to 
computational analyses to identify coding regions and regions bearing DNA sequence 
similarity to known genes. This protocol included the following steps. 

1. Degap the contigs: the sequence contigs often contain symbols (denoted by a 
period symbol) that represent locations where the individual ABI sequence reads have 
insertions or deletions. Prior to automated computational analysis of the contigs, the periods 
were removed. The original data was maintained for future reference. 

2. BAC vector sequences were "masked" within the sequence by using the 
program cross match (Phil Green, http:\\chimera.biotech.washington.edu\UWGC). Since the 
shotgun libraries construction detailed above leaves some BAC vector in the shotgun 
libraries, this program was used to compare the sequence of the BAC contigs to the BAC 
vector and to mask any vector sequence prior to subsequent steps. Masked sequences were 
marked by an "X" in the sequence files, and remained inert during subsequent analyses. 

3. E. coli sequences contaminating the BAC sequences were masked by 
comparing the BAC contigs to the entire E. coli DNA sequence. 

4. Repetitive elements known to be common in the human genome were masked 
using cross match. In this implementation of crossmatch, the BAC sequence was compared 
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to a database of human repetitive elements (Jerzy Jerka, Genetic Information Research 
Institute, Palo Alto, CA). The masked repeats were marked by X and remained inert during 
subsequent analyses. 

5. The location of exons within the sequence was predicted using the MZEF 
5 computer program (Zhang, Proc. Natl. Acad. Sci., 94:565-568 (1997)). 

6. The sequence was compared to the publicly available unigene database 
(National Center for Biotechnology Information, National Library of Medicine, 38A, 8N905, 
8600 Rockville Pike, Bcthesda, MD 20894; www.ncbi.nlm.nih.gov) using the blastn2 
algorithm (Altschul et al, NucL Acids Res., 25:3389-3402 ( 1997)). The parameters for this 

10 search were: E=0.05, v=50, B=50 (where E is the expected probability score cutoff, V is the 
number of database entries returned in the reporting of the results, and B is the number of 
sequence alignments returned in the reporting of the results (Altschul et al, J. MoL Biol, 
215:403-410 (1990)). 

7. The sequence was translated into protein for all six reading frames, and the 
15 protein sequences were compared to a non-redundant protein database compiled from 

Genpept Swissprot PIR (National Center for Biotechnology Information, National Library of 
Medicine, 38A, 8N905, 8600 Rockville Pike, Bethesda, MD 20894; www.ncbi.nlm.nih.gov). 
The parameters for this search were E=0.05, V=50, B= 50, where E, V, and B are defined as 
above. 

20 8. The BAC DNA sequence was compared to the database of the cDNA clones 

derived from direct selection experiments (described below) using blastn2 (Altschul et al, 
NucL Acids. Res., 25:3389-3402 (1997)). The parameters for this search were E=0.05, V-250, 
B=250, where E, V, and B are defined as above. 

9. The BAC sequence was compared to the sequences of all other BACs from the 
25 HBM region on chromosome 1 lql2-13 using blastn2 (Altschul et al, NucL Acids. Res., 

25:3389-3402 (1997)). The parameters for this search were E=0.05, V=50, B=50, where E, V, 
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and B arc defined as above. 

10. The BAC sequence was compared to the sequences derived from the ends of 
BACs from the HBM region on chromosome 1 1 q 1 2- 1 3 using blastn2 (Altschul et al, NucL 
Acids. Res., 25:3389-3402 (1997)). The parameters for this search were E=0.05, V=50, B=50, 

5 where E, V, and B are defined as above. 

1 1 . The BAC sequence was compared to the Genbank database (National Center 
for Biotechnology Information, National Library of Medicine, 38A, 8N905, 8600 Rockville 
Pike, Bethesda, MD 20894; www.ncbi.nlm.nih.gov) using blastn2 (Altschul et al, NucL 
Acids. Res., 25:3389-3402 (1997)). The parameters for this search were E=0.05, V=50, B=50, 

10 where E, V, and B are defined as above. 

12. The BAC sequence was compared to the STS division of Genbank database 
(National Center for Biotechnology Information, National Library of Medicine, 38A, 8N905, 
8600 Rockville Pike, Bethesda, MD 20894; www.ncbi.nlm.nih.gov) using blastn2 (Altschul 
et al, 1997). The parameters for this search were E=0.05, V=50, B= 50, where E, V, and B are 

15 defined as above. 

13. The BAC sequence was compared to the Expressed Sequence (EST) Tag 
Genbank database (National Center for Biotechnology Information, National Library of 
Medicine, 38A, 8N905, 8600 Rockville Pike, Bethesda, MD 20894; www.ncbi.nlm.nih.gov) 
using blastn2 (Altschul et al, NucL Acids. Res., 25:3389-3402 (1997)). The parameters for 

20 this search were E=0.05, V=250, B=250, where E, V, and B are defined as above. 
IX. Gene Identification by Direct cDNA Selection 

Primary linkered cDNA pools were prepared from bone marrow, calvarial bone, 
femoral bone, kidney, skeletal muscle, testis and total brain. Poly (A) + RNA was prepared 
from calvarial and femoral bone tissue (Chomczynski et al, Anal. Biochem., 162:156-159 
25 (1987); D'Alessio et al, Focus, 9:1-4 (1987)) and the remainder of the mRNA was purchased 
from Clontech (Palo Alto, California). In order to generate oligo(dT) and random primed 
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cDNA pools from the same tissue, 2.5 ug mRNA was mixed with oligo(dT) primer in one 
reaction and 2.5 ug mRNA was mixed with random hexamers in another reaction, and both 
were converted to first and second strand cDNA according to manufacturers 
recommendations (Life Technologies, Bethesda, MD). Paired phosphorylated cDNA linkers 

5 (see sequence below) were annealed together by mixing in a 1:1 ratio ( 1 0 ug each) incubated 
i . 

p\ at 65 Y for r,ve ' 1linulcs antl allowed to cool to room temperature. 
Paired linkers oligol/2 

OL1GO 1 : 5'CTG AGC GG A ATT CGT GAG ACC3' (SEQ ID NO: 1 2) 
OL1GO 2: 5'TTG GTC TCA CGT ATT CCG CTC GA3' (SEQ ID NO: 1 3) 



10 



15 



20 



Paired linkers oligo3/4 

OLIGO 3: 5'CTC GAG AAT TCT GGA TCC TC3' (SEQ ID NO: 14) 
OL1GO 4: 5'TTG AGG ATC CAG AAT TCT CGA G3' (SEQ ID NO: 1 5) 



Paired linkers oligo5/6 

OLIGO 5 : 5'TGT ATG CGA ATT CGC TGC GCG3' (SEQ ID NO: 1 6) 
25 OLIGO 6: 5'TTC GCG CAG CGA ATT CGC ATA CA3' (SEQ ID NO: 1 7) 

Paired linkers oligo7/8 
30 OLIGO 7: 5'GTC CAC TGA ATT CTC AGT GAG3' (SEQ ID NO: 1 8) 

OLIGO 8 : 5'TTG TCA CTG AGA ATT CAG TGG AC3' (SEQ ID NO: 1 9) 

35 Paired linkers oligol 1/12 

OLIGO 1 1 : 5'GAA TCC GAA TTC CTG GTC AGC3' (SEQ ID NO:20) 
OLIGO 12: 5'TTG CTG ACC AGG AAT TCG GAT TC3' (SEQ ID NO:2 1 ) 

40 

Linkers were ligated to all oligo(dT) and random primed cDNA pools (see below) according 
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to manufacturers instructions (Life Technologies, Bethesda, MD). 

Oligo 1/2 was ligated to oligo(dT) and random primed cDNA pools prepared from 
bone marrow . Oligo 3 4 was ligated to oligo(dT) and random primed cDNA pools prepared 
from calvarial bone. Oligo 5/6 was ligated to oligo(dT) and random primed cDNA pools 
5 prepared from brain and skeletal muscle. Oligo 7/8 was ligated to oligo(dT) and random 

primed cDNA pools prepared from kidney. Oligo 11/12 was ligated to oligo(dT) and random 
primed cDNA pools prepared from femoral bone. 

The cDNA pools were evaluated for length distribution by PCR amplification using 1 
jLtl of a 1:1, 1:10, and 1 : 100 dilution of the ligation reaction, respectively. PCR reactions were 
10 performed in a Perkin Elmer 9600, each 25 \i\ volume reaction contained 1 jul of DNA, 10 

mM Tns-HCl (pH 8.3), 50 mM KC1, 1.5 mM MgC12, 0.001% gelatin, 200 mM each dNTPs, 

10 |nM primer and 1 unit Taq DNA polymerase (Perkin Elmer) and was amplified under the 

0 O ^ 

(\ following conditions: 30 seconds at 94 jC, 30 seconds at 60 $f and 2 minutes at 72^ for 30 

cycles. The length distribution of the amplified cDNA pools were evaluated by 
15 electrophoresis on a 1 % agarose gel. The PCR reaction that gave the best representation of 
the random primed and oligo(dT) primed cDNA pools was scaled up so that -2-3 jig of each 
cDNA pool was produced. The starting cDNA for the direct selection reaction comprised of 
0.5 jig of random primed cDNAs mixed with 0.5 jig of oligo(dT) primed cDNAs. 

The DNA from the 54 BACs that were used in the direct cDNA selection procedure 
20 was isolated using Nucleobond AX columns as described by the manufacturer (The Nest 
Group, Inc.). 

The BACs were pooled in equimolar amounts and 1 |ig of the isolated genomic DNA 
w as labelled with biotin 16-UTP by nick translation in accordance with the manufacturers 
instructions (Boehringer Mannheim). The incorporation of the biotin was monitored by 
25 methods that could be practiced by one skilled in the art (Del Mastro and Lovett, Methods in 
Molecular Biology, Humana Press Inc., NJ (1996)). 

- 51 - 




Atty. Dkt. No.: 032796-014 
Direct cDNA selection was performed using methods that could be practiced by one 
skilled in the art (Del Mastro and Lovett, Methods in Molecular Biology, Humana Press Inc., 
NJ (1996)). Briefly, the cDNA pools were multiplexed in two separate reactions: In one 
reaction cDNA pools from bone marrow, calvarial bone, brain and testis were mixed, and in 

5 the other cDNA pools from skeletal muscle, kidney and femoral bone were mixed. 

Suppression of the repeats, yeast sequences and plasmid in the cDNA pools was performed to 
a Cot of 20. 100 ng of biotinylated BAC DNA was mixed with the suppressed cDNAs and 
hybridized in solution to a Cot of 200. The biotinylated DNA and the cognate cDNAs was 
captured on streptavidin-coated paramagnetic beads. The beads were washed and the primary 

10 selected cDNAs were eluted. These cDNAs were PCR amplified and a second round of 
direct selection was performed. The product of the second round of direct selection is 
referred to as the secondary selected material. A Galanin cDNA clone, previously shown to 
map to 1 lql2-13 (Evans, Genomics, 18:473-477 (1993)), was used to monitor enrichment 
during the two rounds of selection. 

15 The secondary selected material from bone marrow, calvarial bone, femoral bone, 

kidney, skeletal muscle, testis and total brain was PCR amplified using modified primers of 

oligos 1, 3, 5, 7 and 1 1, shown below, and cloned into the UDG vector pAMPIO (Life 

Technologies, Bethesda, MD), in accordance with the manufacturer's recommendations. 

Modified primer sequences: 

20 Oligol-CUA: 5'CUA CUA CUA CUA CTG AGC GGA ATT CGT GAG ACC3' (SEQ ID 
NO:22) 

01igo3-CUA: 5'CUA CUA CUA CUA CTC GAG AAT TCT GGA TCC TC3' (SEQ ID 
NO:23) 

25 

Ohgo5-CUA: 5'CUA CUA CUA CUA TGT ATG CGA ATT CGC TGC GCG3' (SEQ ID 
NO:24) 

Ohgo7-CUA: 5'CUA CUA CUA CUA GTC CAC TGA ATT CTC AGT GAG3' (SEQ ID 
30 NO:25) 
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Ol.gol 1-CUA: 5'CUA CUA CUA CUA GAA TCC GAA TTC CTG GTC AGC3' (SEQ ID 
NO:26) 

5 The cloned secondary selected material, from each tissue source, was transformed into 

MAX Efficiency DH5a Competent Cells (Life Technologies, Bethesda, MD) as 
recommended by the manufacturer. 384 colonies were picked from each transformed source 
and arrayed into four 96 well microliter plates. 

All secondarily selected cDNA clones were sequenced using Ml 3 dye primer 
10 terminator cycle sequencing kit (Applied Biosystems), and the data collected by the ABI 377 
automated fluorescence sequencer (Applied Biosystems). 

All sequences were analyzed using the BLASTN, BLASTX and FASTA programs 
(Altschul et al, J. Mul. Biol., 215:403-410 (1990), Altschul et al, Nucl. Acids. Res., 25:3389- 
3402 (1997)). The cDNA sequences were compared to a database containing sequences 
15 derived from human repeats, mitochondrial DNA, ribosomal RNA, E. coli DNA to remove 
background clones from the dataset using the program cross match. A further round of 
comparison was also performed using the program BLASTN2 against known genes 
(Genbank) and the BAC sequences from the HBM region. Those cDNAs that were >90% 
homologous to these sequences were filed according to the result and the data stored in a 
20 database for further analysis. cDNA sequences that were identified but did not have 

significant similarity to the BAC sequences from the HBM region or were eliminated by 
cross match were hybridized to nylon membranes which contained the BACs from the HBM 
region, to ascertain whether they hybridized to the target. 

Hybridization analysis was used to map the cDNA clones to the BAC target that 
25 selected them. The BACs that were identified from the HBM region were arrayed and grown 
into a 96 well microtiter plate. LB agar containing 25 ug/ml kanamycin was poured into 96 
well microtiter plate lids. Once the agar had solidified, pre-cut Hybond N+ nylon membranes 



-53- 



Atty. Dkt. No.: 032796-014 

(Amersham) were laid on top of the agar and the BACs were stamped onto the membranes in 
duplicate using a hand held 96 well replica plater (V&P Scientific, Inc.). The plates were 
f\ incubated overnight at 37.4*. The membranes were processed according to the manufacturers 
recommendations. 

5 The cDNAs that needed to be mapped by hybridization were PCR amplified using the 

relevant primer (oligos 1, 3, 5, 7 and 1 1 ) that would amplify that clone. For this PCR 
amplification, the primers were modified to contain a linkered digoxigenin molecule at the 5' 
of the oligonucleotide. The PCR amplification was performed under the same conditions as 
described in Preparation of cDNA Pools (above). The PCR products were evaluated for 
10 quality and quantity by electrophoresis on a 1% agarose gel by loading 5 ul of the PCR 
reaction. The nylon membranes containing the stamped BACs were individually pre- 
hybridized in 50 ml conical tubes containing 10 ml of hybridization solution (5x SSPE, 0.5x 
Blotto 7 5% SDS and 1 mM EDTA (pH 8.0)). The 50 ml conical tubes were placed in a 
rotisserie oven (Robbins Scientific) for 2 hours at 65 £. 25 ng of each cDNA probe was 
15 denatured and added into individual 50 ml conical tubes containing the nylon membrane and 
f\ hybridization solution. The hybridization was performed overnight at 65 / £. The filters were 
ft washed for 20 minutes at 65X in each of the following solutions: 3x SSPE, 0.1% SDS; lx 
SSPE, 0.1% SDS and O.lx SSPE, 0.1% SDS. 

The membranes were removed from the 50 ml conical tubes and placed in a dish. 
20 Acetate sheets were placed between each membrane to prevent them from sticking to each 

other. The incubation of the membranes with the Anti-DIG-AP and CDP-Star was performed 
according to manufacturers recommendations (Boehringer Mannheim). The membranes were 
wrapped in Saran wrap and exposed to Kodak Bio-Max X-ray film for 1 hour. 
X. cDNA Cloning and Expression Analysis 
25 To characterize the expression of the genes identified by direct cDNA selection and 

genomic DNA sequencing in comparison to the publicly available databases, a series of 
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experiments were performed lo further characterize the genes in the HBM region. First, 
oligonucleotide primers were designed for use in the polymerase chain reaction (PCR) so that 
portions of a cDNA, EST, or genomic DNA could be amplified from a pool of DNA 
molecules (a cDNA library) or RNA population (RT-PCR and RACE). The PCR primers 
5 were used in a reaction containing genomic DNA to verify that they generated a product of 
the size predicted based on the genomic (BAC) sequence. A number of cDNA libraries were 
then examined for the presence of the specific cDNA or EST. The presence of a fragment of 
a transcription unit in a particular cDNA library indicates a high probability that additional 
portions of the same transcription unit will be present as well. 
10 A critical piece of data that is required when characterizing novel genes is the length, 

in nucleotides, of the processed transcript or messenger RNA (mRNA). One skilled in the art 
primarily determines the length of an mRNA by Northern blot hybridization (Sambrook et al, 
Molecular Cloning: A Laboratory Manual, Cold Spring Harbor Laboratory, Cold Spring 
Harbor NY (1989)). Groups of ESTs and direct-selected cDNA clones that displayed 
15 significant sequence similarity to sequenced BACs in the critical region were grouped for 
convenience into approximately 30 kilobase units. Within each 30 kilobase unit there were 
from one up to fifty ESTs and direct-selected cDNA clones which comprised one or more 
independent transcription units. One or more ESTs or direct-selected cDNAs were used as 
hybridization probes to determine the length of the mRNA in a variety of tissues, using 
20 commercially available reagents (Multiple Tissue Northern blot; Clontech, Palo Alto, 
California) under conditions recommended by the manufacturer. 

Directionally cloned cDNA libraries from femoral bone, and calvarial bone 
tissue were constructed by methods familiar to one skilled in the art (for example, Soares in 
Automated DNA Sequencing and Analysis, Adams, Fields and Venter, Eds., Academic Press, 
25 NY, pages 110-114 (1 994)). Bones were initially broken into fragments with a hammer, and 
the small pieces were frozen in liquid nitrogen and reduced to a powder in a tissue pulverizer 
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(Spectrum Laboratory Products). RNA was extracted from the powdered bone by 
homogenizing the powdered bone 

with a standard Acid Guanidinium Thiocyanate-Phenol-Chloroform extraction 
butler (e.g. Chomczynski and Sacchi, Anal. Bioclwm., 162:156-159 (1987)) using a polytron 
5 homogenizer (Bnnkman Instruments). Additionally, human brain and lung total RNA was 
purchased from Clontech. PolyA RNA was isolated from total RNA using dynabeads-dT 
according to the manufacturer's recommendations (Dynal, Inc.). 

First strand cDNA synthesis was initiated using an oligonucleotide primer with the 
sequence: 5'-AACTGGAAGAATTC GCGGCCGC AGGAATTTTTTTTTT 
10 TTTTTTTT-3' (SEQ ID NO:27). This primer introduces a NotI restriction site (underlined) 
at the 3' end of the cDNA. First and second strand synthesis were performed using the 
"one-tube" cDNA synthesis kit as described by the anufacturer (Life Technologies, Bethesda, 
MD). Double stranded cDNAs were treated with T4 polynucleotide kinase to ensure that the 
ends of the molecules were blunt (Soares in Automated DNA Sequencing and Analysis, 
15 Adams, Fields and Venter, Eds., Academic Press, NY, pages 110-1 14 (1994)), and the blunt 
ended cDNAs were then size selected by a Biogel column (Huynh et al in DNA Cloning, Vol. 
1, Glover, Ed., IRL Press, Oxford, pages 49-78 (1985)) or with a size-sep 400 sepharose 
column (Pharmacia, catalog # 27-5105-01). Only cDNAs of 400 base pairs or longer were 
used in subsequent steps. EcoRI adapters (sequence: 5' OH-AATTCGGCACGAG-OH 3' 
20 (SEQ ID NO:28), and 5' p-CTCGTGCCG-OH 3' (SEQ ID NO:29)) were then ligated to the 
double stranded cDNAs by methods familiar to one skilled in the art (Soares, 1994). The 
EcoRI adapters were then removed from the 3' end of the cDNA by digestion with NotI 
(Soares, 1994). The cDNA was then ligated into the plasmid vector pBluescript II KS+ 
(Stratagene, La Jolla, California), and the ligated material was transformed into E. coli host 
25 DH 1 0B or DH 1 2S by electroporation methods familiar to one skilled in the art (Soares, 

1994). After growth overnight at 37 <f, DNA was recovered from the E. coli colonies after 
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scraping the plates by processing as directed for the Mega-prep kit (Qiagen, Chatsworth, 
California). The quality of the cDNA libraries was estimated by counting a portion of the 
total numbers of primary transformants and determining the average insert size and the 
percentage of plasmids with no cDNA insert. Additional cDNA libraries (human total brain, 
5 heart, kidney, leukocyte, and fetal brain) were purchased from Life Technologies, Bethesda, 
MD. 

cDNA libraries, both oligo (dT) and random hexamer (N 6 ) primed, were used for 
isolating cDNA clones transcribed within the HBM region: human bone, human brain, human 
kidney and human skeletal muscle (all cDNA libraries were made by the inventors, except for 
10 skeletal muscle (dT) and kidney (dT) cDNA libraries). Four 10 x 10 arrays of each of the 
cDNA libraries were prepared as follows: the cDNA libraries were titered to 2.5 x 10" using 
primary transformants. The appropriate volume of frozen stock was used to inoculate 2 L of 
LB/ampicillin (100 mg/ml). This inoculated liquid culture was aliquotted into 400 tubes of 4 
ml each. Each tube contained approximately 5000 cfu. The tubes were incubated at 30/: 
15 overnight with gentle agitation. The cultures were grown to an OD of 0.7-0.9. Frozen stocks 
were prepared for each of the cultures by aliquotting 100 ul of culture and 300 ul of 80% 
pH glycerol. Stocks were frozen in a dry ice/ethanol bath and stored at -70^ The remaining 
culture was DNA prepared using the Qiagen (Chatsworth, CA) spin miniprep kit according to 
the manufacturer's instructions. The DNAs from the 400 cultures were pooled to make 80 
20 column and row pools. The cDNA libraries were determined to contain HBM cDNA clones 
of interest by PCR. Markers were designed to amplify putative exons. Once a standard PCR 
optimization was performed and specific cDNA libraries were determined to contain cDNA 
clones of interest, the markers were used to screen the arrayed library. Positive addresses 
indicating the presence of cDNA clones were confirmed by a second PCR using the same 
25 markers. 

Once a cDNA library was identified as likely to contain cDNA clones corresponding 
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to a specific transcript of interest from the HBM region, it was manipulated to isolate the 
clone or clones containing cDNA inserts identical to the EST or direct-selected cDNA of 
interest. This was accomplished by a modification of the standard "colony screening" method 
(Sambrook et al, Molecular Cloning: A Laboratory Manual, Cold Spring Harbor Laboratory, 
5 Cold Spring Harbor NY (1989)). Specifically, twenty 150 mm LB+ampicillin agar plates 

were spread with 20,000 colony forming units (cfu) of cDNA library and the colonies allowed 
^ to grow overnight at 37 £. Colonies were transferred to nylon filters (Hybond from 

Amersham, or equivalent) and duplicates prepared by pressing two filters together essentially 
as described (Sambrook et al, Molecular Cloning: A Laboratory Manual, Cold Spring Harbor 
10 Laboratory', Cold Spring Harbor NY (1989)). The "master" plate was then incubated an 
additional 6-8 hours to allow the colonies to grow back. The DNA from the bacterial 
colonies was then affixed to the nylon filters by treating the filters sequentially with 
denaturing solution (0.5 N NaOH, 1 .5 M NaCl) for two minutes, neutralization solution (0.5 
M Tns-CT pH 8.0, 1 .5 M NaCl) for two minutes (twice). The bacterial colonies were 
15 removed from the filters by washing in a solution of 2X SSC/0.1% SDS for one minute while 
rubbing with tissue paper. The filters were air dried and baked under vacuum at SoJ^for 1-2 
hours. 

A cDNA hybridization probe was prepared by random hexamer labeling (Fineberg 
and Vogelstein, Anal. Biochem., 132:6-13 (1983)) or by including gene-specific primers and 
20 no random hexamers in the reaction (for small fragments). Specific activity was calculated 
and was >5X10 8 cpm/10* ug of cDNA. The colony membranes were then prewashed in 10 
Q mM Tris-Cl pH 8.0, 1 M NaCl, 1 mM EDTA, 0.1% SDS for 30 minutes at 55 J?. Following 
the prewash, the filters were prehybridized in > 2 ml/filter of 6X SSC, 50 % deionized 
formamide, 2% SDS, 5X Denhardt's solution, and 100 mg/ml denatured salmon sperm DNA, 
25 (\ at 42 ^for 30 minutes. The filters were then transferred to hybridization solution (6X SSC, 
2% SDS, 5X Denhardt's, 100 mg/ml denatured salmon sperm DNA) containing denatured 
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P, 'P-dCTP-labelled cDN A probe and incubated at 42 for 16-18 hours. 

After the 16-18 hour incubation, the filters were washed under constant agitation in 
pr 2X SSC 2% SDS at room temperature for 20 minutes, followed by two washes at 65 £* for 15 
minutes each. A second wash was performed in 0.5 X SSC, 0.5% SDS for 15 minutes at 
5 A 65 jif. Filters were then wrapped in plastic wrap and exposed to radiographic film for several 
hours to overnight. After film development, individual colonies on plates were aligned with 
the autoradiograph so that they could be picked into a 1 ml solution of LB Broth containing 
^ ampicillin. After shaking at 37 jffor 1-2 hours, aliquots of the solution were plated on 150 
mm plates for secondary screening. Secondary screening was identical to primary screening 
10 (above) except that it was performed on plates containing -250 colonies so that individual 
colonies could be clearly identified for picking. 

After colony screening with radiolabeled probes yielded cDNA clones, the clones 
were characterized by restriction endonuclease cleavage, PCR, and direct sequencing to 
confirm the sequence identity between the original probe and the isolated clone. To obtain 
15 the full-length cDNA, the novel sequence from the end of the clone identified was used to 
probe the library again. This process was repeated until the length of the cDNA cloned 
matches that estimated to be full-length by the northern blot analysis. 

RT-PCR was used as another method to isolate full length clones. The cDNA was 
synthesized and amplified using a "Superscript One Step RT-PCR" kit (Life Technologies, 
20 Gaithersburg, MD). The procedure involved adding 1.5 fig of RNA to the following: 25 [x\ of 
reaction mix provided which is a proprietary buffer mix with MgS0 4 and dNTP's, 1 )il sense 
primer (10 \iM) and 1 |al anti-sense primer (10 |aM), 1 \i\ reverse transcriptase and Taq DNA 
polymerase mix provided and autoclaved water to a total reaction mix of 50 |il. The reaction 

C\ was then placed in a thermocycler for 1 cycle at Sofefov 15 to 30 minutes, then 9A0^ox 15 

c 

2? ft seconds, 55-60^for 30 seconds and 68-72 f for 1 minute per kilobase of anticipated product 
^ and finally 1 cycle of 72^for 5-10 minutes. The sample was analyzed on an agarose gel. 
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The product was excised from the gel and purified from the gel (GcncClean, Bio 101). The 
purified product was cloned in pCTNR (General Contractor DNA Cloning System, 5 Prime - 
3 Prime, Inc.) and sequenced to verify that the clone was specific to the gene of interest. 
Rapid Amplification of cDNA ends (RACE) was performed following the 
5 manufacturer's instructions using a Marathon cDNA Amplification Kit (Clontech, Palo Alto, 
CA) as a method for cloning the 5* and 3' ends of candidate genes. cDNA pools were 
prepared from total RNA bv performing first strand synthesis, where a sample of total RNA 
ft sample was mixed with a modified oligo (dT) primer, heated to 70 cooled on ice and 
followed by the addition of: 5x first strand buffer, 10 mM dNTP mix, and AMV Reverse 
10 ft Transcriptase (20 U/ul). The tube was incubated at 42 ,T for one hour and then the reaction 
tube was placed on ice. For second strand synthesis, the following components were added 
directly to the reaction tube: 5x second strand buffer, 10 mM dNTP mix, sterile water, 20x 
second strand enzyme cocktail and the reaction tube was incubated at 16 C for 1.5 hours. T4 
ft DNA Polymerase was added to the reaction tube and incubated at \6^for 45 minutes. The 
15 second-strand synthesis was terminated with the addition of an EDTA/Glycogen mix. The 
sample was subjected to a phenol/chloroform extraction and an ammonium acetate 
precipitation. The cDNA pools were checked for quality by analyzing on an agarose gel for 
size distribution. Marathon cDNA adapters (Clontech) were then ligated onto the cDNA 
ends. The specific adapters contained priming sites that allowed for amplification of either 5' 
20 or 3' ends, depending on the orientation of the gene specific primer (GSP) that was chosen. 
An aliquot of the double stranded cDNA was added to the following reagents: 10 uM 
Marathon cDNA adapter, 5x DNA ligation buffer, T4 DNA ligase. The reaction was 
ft incubated at 16 ^overnight. The reaction was heat inactivated to terminate the reaction. PCR 
was performed by the addition of the following to the diluted double stranded cDNA pool: 
25 lOx cDNA PCR reaction buffer, 10 uM dNTP mix, 10 uM GSP, 10 uM API primer (kit), 
ft 50x Advantage cDNA Polymerase Mix. Thermal Cycling conditions were 94^for 30 
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[\ seconds, 5 cycles of 94 /for 5 seconds. 72 jf for 4 minutes, 5 cycles of 94 € for 5 seconds, 
ft 70 C^for 4 minutes, 23 cycles of 94 j??or 5 seconds, 68 j^for 4 minutes. After the first round 
of PCR was performed using the GSP to extend to the end of the adapter to create the adapter 
primer binding site, exponential amplification of the specific cDNA of interest was observed. 
5 Usually a second nested PCR is performed to confirm the specific cDNA. The RACE 
product was analyzed on an agarose gel and then excised and purified from the gel 
(GeneClean, BIO 101 ). The RACE product was then cloned into pCTNR (General 
Contractor DNA Cloning System, 5' - 3', Inc.) and the DNA sequence determined to verify 
that the clone is specific to the gene of interest. 
10 XI. Mutation Analysis 

Comparative genes were identified using the above procedures and the exons from 
each gene were subjected to mutation detection analysis. Comparative DNA sequencing was 
used to identify polymorphisms in HBM candidate genes from chromosome 1 lql2-13. DNA 
sequences for candidate genes were amplified from patient lymphoblastoid cell lines. 
15 The inventors developed a method based on analysis of direct DNA sequencing of 

PCR products amplified from candidate regions to search for the causative polymorphism. 
The procedure consisted of three stages that used different subsets of HBM family to find 
segregating polymorphisms and a population panel to assess the frequency of the 
polymorphisms. The family resources result from a single founder leading to the assumption 
20 that all affected individuals will share the same causative polymorphism. 

Candidate regions were first screened in a subset of the HBM family consisting of the 
proband, daughter, and her mother, father and brother. Monochromosomal reference 
sequences were produced concurrently and used for comparison. The mother and daughter 
carried the HBM polymorphism in this nuclear family, providing the ability to monitor 
25 polymorphism transmission. The net result is that two HBM chromosomes and six non-HBM 
chromosomes were screened. This allowed exclusion of numerous frequent alleles. Only 
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alleles exclusively present in the affected individuals passed to the next level of analysis. 

Polymorphisms that segregated exclusively with the HBM phenotypc in this original 
family were then re-examined in an extended portion of the HBM pedigree consisting of two 
additional nuclear families. These families consisted of five HBM and three unaffected 
5 individuals. The HBM individuals in this group included the two critical crossover 
individuals, providing the centromenc and telomeric boundaries of the critical region. 
Tracking the heredity of polymorphisms between these individuals and their affected parents 
allowed for further refining of the critical region. This group brought the total of HBM 
chromosomes screened to seven and the total of non-HBM chromosomes to seventeen. 
10 When a given polymorphism continued to segregate exclusively with the HBM 

phenotype in the extended group, a population panel was then examined. This panel of 84 
persons consisted of 42 individuals known to have normal bone mineral density and 42 
individuals known to be unrelated but with untyped bone mineral density. Normal bone 
mineral density is within two standard deviations of BMD Z score 0. The second group was 
15 from the widely used CEPH panel of individuals. Any segregating polymorphisms found to 
be rare in this population were subsequently examined on the entire HBM pedigree and a 
larger population. 

Polymerase chain reaction (PCR) was used to generate sequencing templates from the 
HBM family's DNA and monochromosomal controls. Enzymatic amplification of genes 

20 within the HBM region on 1 1 q 1 2-1 3 was accomplished using the PCR with oligonucleotides 
flanking each exon as well as the putative 5' regulatory elements of each gene. The primers 
were chosen to amplify each exon as well as 15 or more base pairs within each intron on 
either side of the splice. All PCR primers were made as chimeras to facilitate dye primer 
sequencing. The M13-21F (5'- GTA CGA CGG CCA GT -3') (SEQ ID NO:30) and -28REV 

25 (5'- AAC AGC TAT GAC CAT G -3') (SEQ ID NO:31) primer binding sites were built on to 
the 5' end of each forward and reverse PCR primer, respectively, during synthesis. 150 ng of 
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genomic DNA was used in a 50 ul PCR with 2UAmpliTaq, 500 nM primer and 125 uM 
dNTP. Buffer and cycling conditions were specific to each primer set. TaqStart antibody 
(Clontech) was used for hot start PCR to minimize primer dimer formation. 10% of the 
product was examined on an agarose gel. The appropriate samples were diluted 1 :25 with 
5 deioni/ed water before sequencing. 

Each PCR product was sequenced according to the standard Energy Transfer primer 
(Amersham) protocol. All reactions took place in 96 well trays. 4 separate reactions, one 
each for A, C, G and T were performed for each template. Each reaction included 2 ul of the 
sequencing reaction mix and 3 ul of diluted template. The plates were then heat sealed with 
10 foil tape and placed in a thermal cycler and cycled according to the manufacturer's 

recommendation. After cycling, the 4 reactions were pooled. 3 ul of the pooled product was 
transferred to a new 96 well plate and 1 ul of the manufacturer's loading dye was added to 
each well. All 96 well pipetting procedures occurred on a Hydra 96 pipetting station 
(Robbins Scientific, USA). 1 ul of pooled material was directly loaded onto a 48 lane gel 
15 running on an AB1 377 DNA sequencer for a 1 0 hour, 2.4 kV run. 

Polyphred (University of Washington) was used to assemble sequence sets for 
viewing with Consed (University of Washington). Sequences were assembled in groups 
representing all relevant family members and controls for a specified target region. This was 
done separately for each of the three stages. Forward and reverse reads were included for 
20 each individual along with reads from the monochromosomal templates and a color annotated 
reference sequence. Polyphred indicated potential polymorphic sites with a purple flag. Two 
readers independently viewed each assembly and assessed the validity of the purple-flagged 
sites. 

A total of 23 exons present in the mature mRNA and several other portions of the 
25 primary transcript were evaluated for heterozygosity in the nuclear family of two HBM- 
affected and two unaffected individuals. 25 SNPs were identified, as shown in the table 
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below. 



TABLE 4: Single Nucleotide Polymorphisms in the Zmaxl Gene and E 



Exon Name 


Location 


Ddsc V 11 <i ii 


b2()()e21-h_Contigl_l.nt 


69169 (MJ9UJ 


CIA 


b2()()e2 1 -h Contig4 1 2.nt 


27402 (309U) 


A/Ci 

r\J vJ 


b200e21-h_Contig4_13.nt 


27841 (309G) 


TIC 


b200e2 1 -h_Contig4_l 6.nt 


35600 (309G) 


A ICt 


b200e2 1 -h_Contig4_2 1 .nt 


45619 (309G) 




b200e2 1 -h_Contig4_22.nt-a 


46018 (309G) 


TIC 


b200e21-h__Contig4_22.nt-b 


46093 (309G) 


TIC 


b200e21-h Contig4_22.nt-c 


46190 (309G) 




b200e2 1 -h_Contig4_24.nt-a 


50993 (309G) 




b200e2 1 -h_Contig4_24.nt-b 


51 124 (309G) 


C IT 


b200e2 1 -h_Contig4_25 .nt 


55461 (309G) 


W 1 


b200e21-h_Contig4__33.nt-a 


63645 (309G) 


C 1 A 


b200e2 1 -h_C ontig4_33 .nt-b 


63646 (309G) 


A/C 


b200e2 1 -h_Contig4_6 1 .nt 


24809 (309G) 


T/G 


b200e2 1 -h Contig4_62.nt 


27837 (309G) 


T/C 


b200e2 1 -h_Contig4_63 .nt-a 


31485 (309G) 


C/T 
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b200c2 1 -h_Contig4_63.nt-b 


31683 (309G) 


A/G 


b2()0e2 1 -h_Contig4_9.nt 


24808 (309G) 


T/G 


b527d 1 2-h _Contig030g 1 .nt-a 


31340 (308G) 


T/C 


b527dl 2-h_Contigl)30g_l .nt-b 


32538 (308G) 


A/G 


b527dl2-h_Contig080C_2.nt 


13224 (308G) 


A/G 


b527dl2-h_Contig087C_l.nt 


21119 (308G) 


C/A 


b527dl2-h_Contig087C_4.nt 


30497 (308G) 


G/A 


b527dl 2-h_Contig088C_4.nt 


24811 (309G) 


A/C 


b527dl2-h_Contig089_lHP.nt 


68280 (309G) 


G/A 



In addition to the polymorphisms presented in Table 4, two additional 
polymorphisms can also be present in SEQ ID NO:2. These is a change at position 2002 of 
SEQ ID NO:2. Either a guanine or an adenine can appear at this position. This 
polymorphism is silent and is not associated with any change in the amino acid sequence. 
The second change is at position 4059 of SEQ ID NO:2 corresponding in a cytosine (C) to 
thymine (T) change. This polymorphism results in a corresponding amino acid change from 
a valine (V) to an alanine (A). Other polymorphisms were found in the candidate gene exons 
and adjacent intron sequences. Any one or combination of the polymorphisms listed in Table 
4 or the two discussed above could also have a minor effect on bone mass when present in 
SEQ ID NO:2. 

The present invention encompasses the nucleic acid sequences having the nucleic acid 
sequence of SEQ ID NO: 1 with the above-identified point mutations. 
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Preferably, the present invention encompasses the nucleic acid of SEQ ID NO: 2. 
Specifically, a base-pair substitution changing G to T at position 582 in the coding sequence 
of Zmaxl (the HBM gene) was identified as heterozygous in all HBM individuals, and not 
found in the unaffected individuals (i.e., b527dl2-h_Contig087C_l .nt). Fig. 5 shows the 
5 order of the contigs in B527D12. The direction of transcription for the HBM gene is from 
left to right. The sequence of contig308G of B527D12 is the reverse complement of the 
coding region to the HBM gene. Therefore, the relative polymorphism in contig 308G shown 
in Table 4 as a base change substitution of C to A is the complement to the G to T 
substitution in the HBM gene. This mutation causes a substitution of glycine 171 with valine 
10 (G171V). 

The HBM polymorphism was confirmed by examining the DNA sequence of different 
groups of individuals. In all members of the HBM pedigree (38 individuals), the HBM 
polymorphism was observed in the heterozygous form in affected (i.e., elevated bone mass) 
individuals only (N=18). In unaffected relatives (N=20) (BMDZ<2.0) the HBM 
15 polymorphism was never observed. To determine whether this gene was ever observed in 
individuals outside of the HBM pedigree, 297 phenotyped individuals were characterized at 
the site of the HBM gene. None were heterozygous at the site of the HBM polymorphism. In 
an unphenotyped control group, 1 of 42 individuals was observed to be heterozygous at 
position 582. Since this individual is deceased, their bone mineral density could not be 
20 obtained Taken together, these data prove that the polymorphism observed in the kindred 
(V displaying the high bone mass phenotype is strongly correlated with the GfT polymorphism at 
position 582 of Zmaxl . Taken together, these results establish that the HBM polymorphism 
genetically segregates with the HBM phenotype, and that both the HBM polymorphism and 
phenotype are rare in the general population. 
25 XII. Allele Specific Oligonucleotide (ASO) Analysis 

The amplicon containing the HBM1 polymorphism was PCR amplified using primers 
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specific for the exon of interest. The appropriate population of individuals was PCR 
amplified in 96 well microliter plates as follows. PCR reactions (20 ul) containing IX 
Promega PCR buffer (Cat. # Ml 883 containing 1.5 mM MgCf ), 100 M dNTP, 200 nM PCR 
primers ( 1 863F: CCAAGTTCTGAGAAGTCC and 1 864R: AATACCTGAAACCATAC 
5 CTG), 1 U Amplitaq, and 20 ng of genomic DNA were prepared and amplified under the 
following PCR conditions: 94°C, 1 minute, (94"C, 30 sec; 58"C, 30 sec; 72"C, 1 min.) X35 
cycles), 72"C, 5', 4"C, hold. Loading dye was then added and 10 ul of the products was 
clectrophoresed on 1.5% agarose gels containing 1 ug/ml ethidium bromide at 100-150 V for 
5-10 minutes. Gels were treated 20 minutes in denaturing solution (1.5 M NaCl, 0.5 N 
10 NaOH), and rinsed briefly with water. Gels were then neutralized in 1 M Tris-HCl, pH 7.5, 
1 .5 M NaCl, for 20 minutes and rinsed with water. Gels were soaked in 10 X SSC for 20 
minutes and blotted onto nylon transfer membrane (Hybond N+- Amersham) in 10X SSC 
overnight. Filters were the rinsed in 6X SSC for 10 minutes and UV crosslinked. 

The allele specific oligonucleotides (ASO) were designed with the polymorphism 
15 approximately in the middle. Oligonucleotides were phosphate free at the 5'end and were 
purchased from Gibco BRL. Sequences of the oligonucleotides are: 

2326 Zmaxl .ASO.g: AGACTGGGGTGAGACGC 

2327 Zmaxl. ASO.t: CAGACTGGGTTGAGACGCC 

The polymorphic nucleotides are underlined. To label the oligos, 1 .5 ul of 1 ixg/ul ASO 

20 oligo (2326.Zmaxl. ASO.g or 2327.Zmaxl .ASO.t), 1 1 ul ddH 2 0, 2 ul 10X kinase forward 

buffer, 5 ul y- 32 P-ATP (6000 Ci/mMole), and 1 ul T4 polynucleotide kinase (10 U/ul) were 

mixed and the reaction incubated at 37T for 30-60 minutes. Reactions were then placed at 

At 

Pr 95 'C for 2 minutes and 30M-UO was added. The probes were purified us.ng a G25 microspin 
column (Pharmacia). 

25 Blots were prehybridized in 10 ml 5X SSPE, 5X Denhardt's. 2% SDS, and 100 ug/ml, 

denatured, sonicated salmon sperm DNA at 40"C for 2 hr. The entire reaction mix of kinased 
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ol.go was then added lo 10 ml fresh hybridization buffer (5X SSPE, 5X Denhardts, 2% SDS) and 

hvbridi/ed at 40"C for at least 4 hours to overnight. 

All washes done in 5X SSPE, 0. 1 % SDS. The first wash was at 45"C for 1 5 minutes; 

the solution was then changed and the filters washed 50T for 15 minutes. Filters were then 
5 exposed lo Kodak biomax film with 2 intensifying screens at -70"C for 15 minutes to 1 hr. If 

necessary the filters were washed at 55"C for 15 minutes and exposed to film again. Filters 

were stripped by washing in boiling 0. 1 X SSC, 0. 1 % SDS for 1 0 minutes at least 3 times. 

The two films that best captured the allele specific assay with the 2 ASOs were converted 

into digital images by scanning them into Adobe PhotoShop. These images were overlaid 
10 against each other in Graphic Converter and then scored and stored in FileMaker Pro 4.0 (see 

Fig. 9). 

XIII. Cellular Localization of Zmaxl 

A. Gene Expression in Rat tibia by non isotopic In Situ Hybridization 

In situ hybridization was conducted by Pathology Associates International (PAI), 
15 Frederick, MD. This study was undertaken to determine the specific cell types that express 
the Zmaxl gene in rat bone with particular emphasis on areas of bone growth and 
remodeling. Zmaxl probes used in this study were generated from both human (HuZmaxl) 
and mouse (MsZmaxl ) cDNAs, which share an 87% sequence identity. The homology of 
human and mouse Zmaxl with rat Zmaxl is unknown. 
20 For example, gene expression by non-isotopic in situ hybridization was performed as 

follows, but other methods would be known to the skilled artisan. Tibias were collected from 
two 6 to 8 week old female Sprague Dawley rats euthanized by carbon dioxide asphyxiation. 
Distal ends were removed and proximal tibias were snap frozen in OCT embedding medium 
with liquid nitrogen immediately following death. Tissues were stored in a -80 C freezer. 
25 Probes for amplifying PCR products from cDNA were prepared as follows. The 

primers to amplify PCR products from a cDNA clone were chosen using published sequences 
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of both human LRP5 (Genbank Accession No. AB017498) and mouse LRP5 (Genbank 
Accession No. AF064984). In order to minimize cross reactivity with other genes in the 
I DL receptor family, the PCR products w ere derived from an intracellular portion of the 
protein coding region. PCR was performed in a 50 ul reaction volume using cDNA clone as 
5 template. PCR reactions contained 1 .5 mM MgCK, 1 unit Amplitaq, 200 uM dNTPs and 2 
uM each primer. PCR cycling conditions were 94°C for 1 min., followed by 35 cycles of 
94°C for 30 seconds, 55°C for 30 seconds, 72°C for 30 seconds; followed by a 5 minute 
extension at 72°C. The reactions were then run on a 1 .5 % agarose Tris-Acetate gel. DNA 
was eluted from the agarose, ethanol precipitated and resuspended in 10 mM Tris, pH 8.0. 
10 Gel purified PCR products were prepared for both mouse and human cDNAs and supplied to 
Pathology Associates International for in situ hybridizations. 

The sequence of the human and mouse PCR primers and products were as follows: 
Human Zmax 1 sense pri mer (HBM1253) 

CCCGTGTGCTCCGCCGCCCAGTTC 
15 Human Zmax 1 antisense primer (HBM1465) 

GGCTCACGGAGCTCATCATGGACTT 

Human Zmaxl PCR product 

CCCGTGTGCTCCGCCGCCCAGTTCCCCTGCGCGCGGGGTCAGTGTGTGGACCTGC 
GCCTGCGCTGCGACGGCGAGGCAGACTGTCAGGACCGCTCAGACGAGGTGGACT 
GTGACGCCATCTGCCTGCCCAACCAGTTCCGGTGTGCGAGCGGCCAGTGTGTCCT 
CATCAAACAGCAGTGCGACTCCTTCCCCGACTGTATCGACGGCTCCGACGAGCTC 
ATGTGTGAAATCACCAAGCCGCCCTCAGACGACAGCCCGGCCCACAGCAGTGCC 
ATCGGGCCCGTCATTGGCATCATCCTCTCTCTCTTCGTCATGGGTGGTGTCTATTT 
TGTGTGCCAGCGCGTGGTGTGCCAGCGCTATGCGGGGGCCAACGGGCCCTTCCCG 
CACGAGTATGTCAGCGGGACCCCGCACGTGCCCCTCAATTTCATAGCCCCGGGCG 



20 



25 



-69- 



Atty. Dkt. No.: 032796-014 

GTTCCCAGCATGGCCCCTTCACAGGCATCGCATGCGGAAAGTCCATGATGAGCTC 
CGTGAGCC 

Mouse Zmax 1 Sense prim er (HBM1655) 
AGCGAGGCCACCATCCACAGG 
5 Mouse /.max 1 antisense primer (HBM 1656) 
TCGCTGGTCGGCATAATCAAT 
Mouse Zmaxl PGR product 

AGCAGAGCCACCATCCACAGGATCTCCCTGGAGACTAACAACAACGATGTGGCT 
ATCCCACTCACGGGTGTCAAAGAGGCCTCTGCACTGGACTTTGATGTGTCCAACA 

10 ATCACATCTACTGGACTGATGTTAGCCTCAAGACGATCAGCCGAGCCTTCATGAA 
TGGGAGCTCAGTGGAGCACGTGATTGAGTTTGGCCTCGACTACCCTGAAGGAATG 
GCTGTGGACTGGATGGGCAAGAACCTCTATTGGGCGGACACAGGGACCAACAGG 
ATTGAGGTGGCCCGGCTGGATGGGCAGTTCCGGCAGGTGCTTGTGTGGAGAGAC 
CTTGACAACCCCAGGTCTCTGGCTCTGGATCCTACTAAAGGCTACATCTACTGGA 

15 CTGAGTGGGGTGGCAAGCCAAGGATTGTGCGGGCCTTCATGGATGGGACCAATT 
GTATGACACTGGTAGACAAGGTGGGCCGGGCCAACGACCTCACCATTGATTATG 

CCGACCAGCGA 

Riboprobes were synthesized as follows. The PCR products were reamplified with 
chimeric primers designed to incorporate either a T3 promoter upstream, or a T7 promoter 

20 downstream of the reamplif.cation products. The resulting PCR products were used as 
template to synthesize digoxigenin-labeled riboprobes by in vitro transcription (IVT). 
Antisense and sense riboprobes were synthesized using T7 and T3 RNA polymerases, 
respectively, in the presence of digoxigenin-1 1-UTP (Boehringer-Mannheim) using a 
MAXIscnpt IVT kit (Ambion) according to the manufacturer. The DNA was then degraded 

25 with Dnase- 1 , and unincorporated digoxigenin was removed by ultrafiltration. Riboprobe 
integrity was assessed by electrophoresis through a denaturing polyacrylamide gel. 
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Molecular size was compared with the electrophoretic mobility of a 100-1000 base pair (bp) 
RNA ladder (Ambion). Probe yield and labeling was evaluated by blot immunochemistry. 
Riboprobes were stored in 5 ul aliquots at 80 C. 

The //; situ hybridization was performed as follows. Frozen rat bone was cut into 5 
5 uM sections on a Jung CM3000 cryostat (Leica) and mounted on adhesive slides 

(Instrumedics). Sections were kept in the cryostat at 20 C until all the slides were prepared 
in order to prevent mRNA degradation prior to post-fixation for 15 minutes in 4% 
paraformaldehyde. Following post- fixation, sections were incubated with 1 ng/ul of either 
antisense or sense riboprobe in Pathology Associates International (PA1) customized 
10 hybridization buffer for approximately 40 hours at 58 C. Following hybridization, slides 
were subjected to a series of post-hybridization stringency washes to reduce nonspecific 
probe binding. Hybridization was visualized by immunohistochemistry with an anti- 
digoxigenin antibody (FAB fragment) conjugated to alkaline phosphatase. Nitroblue 
tetrazolium chlonde/bromochloroindolyl phosphate (Boehringer-Mannheim), a precipitating 
15 alkaline phosphatase substrate, was used as the chromogen to stain hybridizing cells purple to 
nearly black, depending on the degree of staining. Tissue sections were counter-stained with 
nuclear fast red. Assay controls included omission of the probe, omission of probe and anti- 

digoxigenin antibody. 

Specific cell types were assessed for demonstration of hybridization with antisense 

20 probes by visualizing a purple to black cytoplasmic and/or peri-nuclear staining indicating a 
positive hybridization signal for mRNA. Each cell type was compared to the replicate 
sections, which were hybridized with the respective sense probe. Results were considered 
positive if staining was observed with the antisense probe and no staining or weak 
background with the sense probe. 

25 The cellular localization of the hybridization signal for each of the study probes is 

summarized in Table 5. Hybridization for Zmaxl was primarily detected in areas of bone 
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involved in remodeling, including the endostcum and trabecular bone within the metaphysis. 
Hybridization in selected bone lining cells of the periosteum and epiphysis were also 
observed. Positive signal was also noted in chondrocytes within the growth plate, particularly 
in the proliferating chondrocytes. See Figs. 10, 11 and 12 for representative 
5 photomicrographs of /'/; situ hybridization results. 



TABLE 5 

Summary of Zmaxl in situ hybridization in rat tibia 
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These studies confirm the positional expression of Zmaxl in cells involved in bone 
5 remodeling and bone formation. Zmaxl expression in the zone of proliferation and in the 
osteoblasts and osteoclasts of the proximal metaphysis, suggests that the Zmaxl gene is 
involved in the process of bone growth and mineralization. The activity and differentiation of 
osteoblasts and osteoclasts are closely coordinated during development as bone is formed and 
during growth as well as in adult life as bone undergoes continuous remodeling. The 
10 formation of internal bone structures and bone remodeling result from the coupling of bone 

resorption by activated osteoclasts with subsequent deposition of new material by osteoblasts. 
Zmaxl is related to the LDL receptor gene, and thus may be a receptor involved in 
mechanosensation and subsequent signaling in the process of bone remodeling. Therefore, 
changes in the level of expression of this gene could impact on the rate of remodeling and 
15 degree of mineralization of bone. 
XIV. Antisense 

Antisense oligonucleotides are short synthetic nucleic acids that contain 
complementary base sequences to a targeted RNA. Hybridization of the RNA in living cells 
with the antisense oligonucleotide interferes with RNA function and ultimately blocks protein 
20 expression. Therefore, any gene for which the partial sequence is known can be targeted by 
an antisense oligonucleotide. 

Antisense technology is becoming a widely used research tool and will play an 
increasingly important role in the validation and elucidation of therapeutic targets identified 
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by genomic sequencing efforts. 

Antisense technology was developed to inhibit gene expression by utilizing an 
ol.gonucleot.de complementary to the mRNA that encodes the target gene. There are several 
possible mechanisms for the inhibitory effects of antisense oligonucleotides. Among them, 
5 degradation of mRNA by RNase H is considered to be the major mechanism of inhibition of 
protein function. This technique was originally used to elucidate the function of a target 
gene, but may also have therapeutic applications, provided it is designed carefully and 
properly. 

An example of materials and methods for preparing antisense oligonucleotides can be 
10 performed as follows. Preliminary studies have been undertaken in collaboration with 

Sequiter (Natick, MA) using the antisense technology in the osteoblast-like murine cell line, 
MC3T3. These cells can be triggered to develop along the bone differentiation sequence. An 
initial proliferation period is characterized by minimal expression of differentiation markers 
and initial synthesis of collagenous extracellular matrix. Collagen matrix synthesis is 
15 required for subsequent induction of differentiation markers. Once the matrix synthesis 
begins, osteoblast marker genes are activated in a clear temporal sequence: alkaline 
phosphatase is induced at early times while bone sialoprotien and osteocalcin appear later in 
the differentiation process. This temporal sequence of gene expression is useful in 
monitoring the maturation and mineralization process. Matrix mineralization, which does not 
20 begin until several days after maturation has started, involves deposition of mineral on and 
within collagen fibrils deep within the matrix near the cell layer-culture plate interface. The 
collagen fibril-associated mineral formed by cultured osteoblasts resembles that found in 
woven bone in vivo and therefore is used frequently as a study reagent. 

MC3T3 cells were transfected with antisense oligonucleotides for the first week of the 
25 differentiation, according to the manufacturer's specifications (U.S. Patent No. 5,849,902). 
The oligonucleotides designed for Zmaxl are given below: 
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10875: AGUACAGCUUCUUGCCAACCCAGUC 
10876: UCCUCCAGGUCGAUGGUCAGCCCAU 
10877: GUCUGAGUCCGAGUUCAAAUCCAGG 

Figure 13 shows the results of antisense inhibition of Zmaxl in MC3T3 cells. The three 
5 oligonucleotides shown above were transfected into MC3T3 and RNA was isolated according 
to standard procedures. Northern analysis clearly shows markedly lower steady state levels 
of the Zmaxl transcript while the control gene GAPDH remained unchanged. Thus, 
antisense technology using the primers described above allows for the study of the role of 
Zmaxl expression on bone biology. 
10 The protein encoded by Zmax 1 is related to the Low Density Lipoprotein receptor 

(LDL receptor). See, Goldstein et al, Ann. Rev. Cell Biology, 1 : 1-39 (1985); Brown et al, 
Science, 232:34-47 (1986). The LDL receptor is responsible for uptake of low density 
l.poprotein, a hpid-protein aggregate that includes cholesterol. Individuals with a defect in 
the LDL receptor are deficient in cholesterol removal and tend to develop artherosclerosis. In 
15 addition, cells with a defective LDL receptor show increased production of cholesterol, in 
part because of altered feedback regulation of cholesterol synthetic enzymes and in part 
because of increased transcription of the genes for these enzymes. In some cell types, 
cholesterol is a precursor for the formation of steroid hormones. 

Thus, the LDL receptor may, directly or indirectly, function as a signal transduction 
20 protein and may regulate gene expression. Because Zmax 1 is related to the LDL receptor, 
this protein may also be involved in signaling between cells in a way that affects bone 
remodeling. 

The glycine 171 amino acid is likely to be important for the function of Zmaxl 
because this amino acid is also found in the mouse homologue of Zmaxl. The closely related 
25 LRP6 protein also contains glycine at the corresponding position (Brown et al, Biochemical 
and Biophysical Research Comm., 248:879-888 (1988)). Amino acids that are important in a 
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protein's structure or function tend to be conserved between species, because natural selection 
prevents mutations with altered amino acids at important positions from arising. 

In addition, the extracellular domain of Zmaxl contains four repeats consisting of five 
YWT motifs followed by an EFG motif. This SYWTtEGF repeat is likely to form a distinct 
5 folded protein domain, as this repeat is also found in the LDL receptor and other LDL 
receptor-related proteins. The first three 5YWT+EGF repeats are very similar in their 
structure, while the fourth is highly divergent. Glycine 171 occurs in the central YWT motif 
of the first 5YWT+EGF repeat in Zmaxl . The other two similar 5YWT+EGF repeats of 
Zmaxl also contain glycine at the corresponding position, as does the 5YWT+EGF repeat in 
10 the LDL receptor protein. However, only 1 7.6% of the amino acids are identical among the 
first three 5YWT+EGF repeats in Zmaxl and the single repeat in the LDL receptor. These 
observations indicate that glycine 171 is essential to the function of this repeat, and mutation 
of glycine 171 causes a functional alteration of Zmaxl. The cDNA and peptide sequences are 
shown in Figs. 6A-6E. The critical base at nucleotide position 582 is indicated in bold and is 
15 underlined. 

Northern blot analysis (Figs. 7A-B) reveals that Zmaxl is expressed in human bone 
tissue as well as numerous other tissues. A multiple-tissue Northern blot (Clontech, Palo 
Alto, CA) was probed with exons from Zmaxl . As shown in Fig. 7A, the 5.5 kb Zmaxl 
transcript was highly expressed in heart, kidney, lung, liver and pancreas and is expressed at 
20 lower levels in skeletal muscle and brain. A second northern blot, shown in Fig. 7B, 

confirmed the transcript size at 5.5 kb, and indicated that Zmaxl is expressed in bone, bone 
marrow, calvaria and human osteoblastic cell lines. 

Taken together, these results indicate that the HBM polymorphism in the Zmaxl gene 
is responsible for the HBM phenotype, and that the Zmaxl gene is important in bone 
25 development. In addition, because mutation of Zmaxl can alter bone mineralization and 
development, it is likely that molecules that bind to Zmaxl may usefully alter bone 
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development. Such molecules may include, for example, small molecules, proteins, RNA 
aptamers, peptide aptamers. and the like. 

XV . Preparation of Nucleic Acids, Vectors, Transformations and Host Cells 

Large amounts of the nucleic acids of the present invention may be produced by 
5 replication in a suitable host cell. Natural or synthetic nucleic ac.d fragments coding for a 
desired fragment will be incorporated into recombinant nucleic acid constructs, usually DNA 
constructs, capable of introduction into and replication in a prokaryotic or eukaryotic cell. 
Usually the nucleic acid constructs will be suitable for replication in a unicellular host, such 
as yeast or bacteria, but may also be intended for introduction to (with and without integration 
10 within the genome) cultured mammalian or plant or other eukaryotic cell lines. The 

purification of nucleic acids produced by the methods of the present invention is described, 
for example, in Sambrook et al, Molecular Cloning. A Laboratory Manual, 2nd Ed. (Cold 
Spnng Harbor Laboratory, Cold Spring Harbor, NY (1989) or Ausubel et al, Current 
Protocols in Molecular Biology, J. Wiley and Sons, NY (1992). 
15 The nucleic acids of the present invention may also be produced by chemical 

synthesis, e.g., by the phosphoramidite method described by Beaucage et al, Tetra. Letts., 
22:1859-1862 (1981) or the triester method according to Matteucci,et al, J. Am. Chem. Soc, 
103:3185 (1981), and may be performed on commercial, automated oligonucleotide 
synthesizers. A double-stranded fragment may be obtained from the single-stranded product 
20 of chemical synthesis either by synthesizing the complementary strand and annealing the 
strands together under appropriate conditions or by adding the complementary strand using 
DNA polymerase with an appropriate primer sequence. 

Nucleic acid constructs prepared for introduction into a prokaryotic or eukaryotic host 
may comprise a replication system recognized by the host, including the intended nucleic acid 
25 fragment encoding the des.red protein, and will preferably also include transcription and 

translations initiation regulatory sequences operably linked to the protein encoding segment. 
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Expression vectors may include, for example, an origin of replication or autonomously 
replicating sequence ( ARS) and expression control sequences, a promoter, an enhancer and 
necessary processing information sites, such as ribosome-binding sites, RNA splice sites, 
polyadenylation s.tes, transcriptional terminator sequences, and mRNA stabilizing sequences. 
5 Secret.on signals may also be included where appropriate, whether from a native HBM or 
Zmaxl protein or from other receptors or from secreted proteins of the same or related 
species, which allow the protein to cross and/or lodge in cell membranes, and thus attain its 
functional topology, or be secreted from the cell. Such vectors may be prepared by means of 
standard recombinant techniques well known in the art and discussed, for example, in 
10 Sambrook et al, Molecular Cloning. A Laboratory Manual, 2nd Ed. (Cold Spring Harbor 

Laboratory, Cold Spring Harbor, NY (1989) or Ausubel et al, Current Protocols in Molecular 
Biology, J. Wiley and Sons, NY (1992). 

An appropriate promoter and other necessary vector sequences will be selected so as 
to be functional in the host, and may include, when appropriate, those naturally associated 
15 with Zmaxl or HBM genes. Examples of workable combinations of cell lines and expression 
vectors are described in Sambrook et al, Molecular Cloning. A Laboratory Manual, 2nd Ed. 
(Cold Spring Harbor Laboratory, Cold Spring Harbor, NY (1989) or Ausubel et al, Current 
Protocols in Molecular Biology, J. Wiley and Sons, NY (1992). Many useful vectors are 
known in the art and may be obtained from such vendors as Stratagene, New England 
20 BioLabs, Promega Biotech, and others. Promoters such as the trp, lac and phage promoters, 
tRNA promoters and glycolytic enzyme promoters may be used in prokaryotic hosts. Useful 
yeast promoters include promoter regions for metallothionein, 3-phosphoglycerate kinase or 
other glycolytic enzymes such as enolase or glyceraldehyde-3-phosphate dehydrogenase, 
enzymes responsible for maltose and galactose utilization, and others. Vectors and promoters 
25 suitable for use in yeast expression are further described in EP 73,675A. Appropriate non- 
native mammalian promoters might include the early and late promoters from SV40 (Fiers et 
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al, Xanax; 273:1 13 (1978)) or promoters derived from murine Moloney leukemia virus, 
mouse tumor virus, avian sarcoma viruses, adenovirus II, bovine papilloma virus or polyoma. 
In addition, the construct may be joined to an amplifiable gene (e.g., DHFR) so that multiple 
copies ofthe gene may be made. For appropriate enhancer and other expression control 
5 sequences, see also Enhancers and Eukaryotic Gene Expression, Cold Spring Harbor Press, 
Cold Spring Harbor, NY (1983). 

While such expression vectors may replicate autonomously, they may also replicate 
by being inserted into the genome ofthe host cell, by methods well known in the art. 
Expression and cloning vectors will likely contain a selectable marker, a gene 
10 encoding a protein necessary for survival or growth of a host cell transformed with the vector. 
The presence of this gene ensures growth of only those host cells which express the inserts. 
Typical selection genes encode proteins that a) confer resistance to antibiotics or other toxic 
substances, e.g. ampicillin, neomycin, methotrexate, etc.; b) complement auxotrophic 
deficiencies, or c) supply critical nutrients not available from complex media, e.g., the gene 
15 encoding D-alanine racemase for Bacilli. The choice of the proper selectable marker will 

depend on the host cell, and appropriate markers for different hosts are well known in the art. 

The vectors containing the nucleic acids of interest can be transcribed in vitro, and the 
resulting RNA introduced into the host cell by well-known methods, e.g., by injection (see, 
Kubo et al, FEBS Letts. 241 :1 19 (1988)), or the vectors can be introduced directly into host 
20 cells by methods well known in the art, which vary depending on the type of cellular host, 
including electroporation; transfection employing calcium chloride, rubidium chloride, 
calcium phosphate, DEAE-dextran, or other substances; microprojectile bombardment; 
lipofection; infection (where the vector is an infectious agent, such as a retroviral genome); 
and other methods. See generally, Sambrook et al.. 1989 and Ausubel et al., 1992. The 
25 introduction of the nucleic acids into the host cell by any method known in the art, including 
those described above, will be referred to herein as "transformation." The cells into which 
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have been introduced nucleic acids described above are meant to also include the progeny of 
such cells. 

Large quantities of the nucleic acids and proteins of the present invention may be 
prepared by expressing the Zmaxl or HBM nucleic acids or portions thereof in vectors or 
5 other expression vehicles in compatible prokaryotic or eukaryotic host cells. The most 

commonly used prokaryotic hosts are strains of Escherichia coli, although other prokaryotes, 
such as Bacillus subulis or Pseiulomonas may also be used. 

Mammalian or other eukaryotic host cells, such as those of yeast, filamentous fungi, 
plant, insect, or amphibian or avian species, may also be useful for production of the proteins 
10 of the present invention. Propagation of mammalian cells in culture is per se well known. 
See, Jakoby and Pastan (eds.), Cell Culture. Methods in Enzymology, volume 58, Academic 
Press, Inc., Harcourt Brace Jovanovich, NY, (1979)). Examples of commonly used 
mammalian host cell lines are VERO and HeLa cells, Chinese hamster ovary (CHO) cells, 
and W138, BHK, and COS cell lines, although it will be appreciated by the skilled 
15 practitioner that other cell lines may be appropriate, e.g., to provide higher expression 
desirable glycosylation patterns, or other features. 

Clones are selected by using markers depending on the mode of the vector 
construction. The marker may be on the same or a different DNA molecule, preferably the 
same DNA molecule. In prokaryotic hosts, the transformant may be selected, e.g., by 
20 resistance to ampicillin, tetracycline or other antibiotics. Production of a particular product 
based on temperature sensitivity may also serve as an appropriate marker. 

Prokaryotic or eukaryotic cells transformed with the nucleic acids of the present 
invention will be useful not only for the production of the nucleic acids and proteins of the 
present invention, but also, for example, in studying the characteristics of Zmaxl or HBM 
25 proteins. 

Antisense nucleic acid sequences arc useful in preventing or diminishing the 
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expression of Zmaxl or HBM, as will be appreciated by one skilled in the art. For example, 
nucleic ac.d vectors containing all or a portion of the Zmaxl or HBM gene or other sequences 
from the Zmaxl or HBM region may be placed under the control of a promoter in an 
antisense orientation and introduced into a cell. Expression of such an antisense construct 
5 within a cell will interfere with Zmaxl or HBM transcription and/or translation and/or 
replication. 

The probes and primers based on the Zmaxl and HBM gene sequences disclosed 
herein are used to identify homologous Zmaxl and HBM gene sequences and proteins in 
other species. These Zmaxl and HBM gene sequences and proteins are used in the 
10 diagnostic/prognostic, therapeutic and drug screening methods described herein for the 
species from which they have been isolated. 
XVI. Protein Expression and Purification 

Expression and purification of the HBM protein of the invention can be performed 
essentially as outlined below. To facilitate the cloning, expression and purification of 
15 membrane and secreted protein from the HBM gene, a gene expression system, such as the 
pET System (Novagen), for cloning and expression of recombinant proteins in E. coli was 
selected. Also, a DNA sequence encoding a peptide tag, the His-Tap, was fused to the 3' end 
of DNA sequences of interest to facilitate purification of the recombinant protein products. 
The 3' end was selected for fusion to avoid alteration of any 5' terminal signal sequence. 
20 Nucleic acids chosen, for example, from the nucleic acids set forth in SEQ ID NOS: 1 , 

3 and 5-12 for cloning HBM were prepared by polymerase chain reaction (PCR). Synthetic 
oligonucleotide primers specific for the 5' and 3' ends of the HBM nucleotide sequence were 
designed and purchased from Life Technologies (Gaithersburg, MD). All forward primers 
(specific for the 5' end of the sequence) were designed to include an Ncol cloning site at the 5' 
25 terminus. These primers were designed to permit initiation of protein translation at the 

methionine residue encoded within the Ncol site followed by a valine residue and the protein 
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encoded by the HBM DNA sequence. All reverse primers (specific for the 3' end of the 
sequence) included an EcoRl site at the 5' terminus to permit cloning of the HBM sequence 
into the reading frame of the pET-28b. The P ET-28b vector provided a sequence encoding an 
additional 2() carboxyl-termmal amino acids including s.x histid.ne residues (at the C- 
5 terminus), which comprised the histidine affinity tag. 

Genomic DNA prepared from the HBM gene was used as the source of template DNA 
for PCR amplification (Ausubel et al, Current Protocols in Molecular Biology', John Wiley & 
Sons (1994)). To amplify a DNA sequence containing the HBM nucleotide sequence, 
genomic DNA (50 ng) was introduced into a reaction vial containing 2 mM MgCU, 1 M 
10 synthetic oligonucleotide primers (forward and reverse primers) complementary to and 

flanking a defined HBM, 0.2 mM of each of deoxynucleotide triphosphate, dATP, dGTP, 
dCTP, dTTP and 2.5 units of heat stable DNA polymerase (Amplitaq, Roche Molecular 
Systems, Inc., Branchburg, NJ) in a final volume of 100 microliters. 

Upon completion of thermal cycling reactions, each sample of amplified DNA was 
15 purified using the Qiaquick Spin PCR purification kit (Qiagen, Gaithersburg, MD). All 

amplified DNA samples were subjected to digestion with the restriction endonucleases, e.g., 
Ncol and EcoRI (New England BioLabs, Beverly, MA) (Ausubel et al, Current Protocols in 
Molecular Biology, John Wiley & Sons, Inc. (1994)). DNA samples were then subjected to 
electrophoresis on 1.0% NuSeive (FMC BioProducts, Rockland, ME) agarose gels. DNA 
20 was visualized by exposure to ethidium bromide and long wave UV irradiation. DNA 

contained in slices isolated from the agarose gel was purified using the Bio 101 GeneClean 
Kit protocol (Bio 101, Vista, CA). 

The pET-28b vector was prepared for cloning by digestion with restriction 
endonucleases, e.g., Ncol and EcoRI (New England BioLabs, Beverly, MA) (Ausubel et al, 
25 Current Protocols in Molecular Biology, John Wiley & Sons, Inc. ( 1 994)). The pET-28a 
vector, which encodes the histidine affinity tag that can be fused to the 5' end of an inserted 
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gene, was prepared by digestion with appropriate restriction endonucleases. 

Following digestion, DNA inserts were cloned (Ausubel et al, Current Protocols in 
Molecular Biology: John Wiley & Sons. Inc. (1994)) into the previously digested pET-28b 
expression vector. Products of the ligation reaction were then used to transform the BL21 
5 stra.n of E. coli (Ausubel et al, Current Protocols in Molecular Btologx; John Wiley & Sons, 
Inc. (1994)) as described below. 

Competent bacteria, E. coli strain BL21 or E. coli strain BL21 (DE3), were 
transformed with recombinant P ET expression plasmids carrying the cloned HBM sequence 
according to standard methods (Ausubd et al, Current Protocols in Modular Biology, John 
10 ft Wiley & Sons, Inc. (1994)). Bnefly^fof ligation reaction was mixed wiuWf of 
electrocompetent cells and subjected to a high voltage pulse, after which samples were 
incubated in 0.45 ml SOC medium (0.5% yeast extract, 2.0% tryptone, 10 mM NaCl, 2.5 mM 
ft KC1. 10 mM MgCL, 10 mM MgS0 4 and 20 mM glucose) at 37^vith shaking for 1 hour. 
Samples were then spread on LB agar plates containing 25 ug/ml kanamycin sulfate for 
15 growth overnight. Transformed colonies of BL21 were then picked and analyzed to evaluate 

cloned inserts, as described below. 

Individual BL21 clones transformed with recombinant pET-28b HBM nucleotide 
sequences were analyzed by PCR amplification of the cloned inserts using the same forward 
and reverse primers specific for the HBM sequences that were used in the original PCR 
20 amplification cloning reactions. Successful amplification verifies the integrate of the HBM 
sequence in the expression vector (Ausubel et al, Current Protocols in Molecular Biology, 

John Wiley & Sons, Inc. (1994)). 

Individual clones of recombinant pET-28b vectors carrying properly cloned HBM 
nucleotide sequences were picked and incubated in 5 ml of LB broth plus 25 ug/ml 
25 kanamycin sulfate overnight. The following day plasmid DNA was isolated and purified 
using the Qiagen plasmid purification protocol (Qiagen Inc., Chatsworth, CA). 
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The pET vector can be propagated in any E. coli K-12 strain, e.g., HMS174, HB101, 
JM 1 09. DH5 and the like, for purposes of cloning or plasmid preparation. Hosts for 
expression include E. coli strains containing a chromosomal copy of the gene for T7 RNA 
polymerase. These hosts were lysogens of bacteriophage DE3, a lambda derivative that 
5 carr.es the lacl gene, the lacUVS promoter and the gene for T7 RNA polymerase. T7 RNA 
^. polymerase was induced by addition ot^soprepyl-^D-th^alactosrcfc (IPTG), and the T7 
RNA polymerase transcribes any target plasmid containing a functional T7 promoter, such as 
pET-28b. carrying its gene of interest. Strains include, for example, BL21(DE3) (Studier et 
al, Mcth. Enzymol, 185:60-89 (1990)). 
10 To express the recombinant HBM sequence, 50 ng of plasmid DNA are isolated as 

described above to transform competent BL21(DE3) bacteria as described above (provided by 
\o\ agen as part of the pET expression kit). The lacZ gene (JVgalactosidase) is expressed in 
the pET-System as described for the HBM recombinant constructions. Transformed cells 
were cultured in SOC medium for 1 hour, and the culture was then plated on LB plates 
15 containing 25 ug/ml kanamycin sulfate. The following day, the bacterial colonies were 
pooled and grown in LB medium containing kanamycin sulfate (25 ug/ml) to an optical 
density at 600 nM of 0.5 to 1.0 O.D. units, at which point 1 mM IPTG was added to the 
culture for 3 hours to induce gene expression of the HBM recombinant DNA constructions. 
After induction of gene expression with IPTG, bacteria were collectedly 
20 ft centrifugation in a Sorvall RC-3B centrifuge at 3500 x g for 15 minutes at Pellets were 
resuspended in 50 ml of cold mM Tris-HCl, pH 8.0, 0.1 M NaCl and 0.1 mM EDTA (STE 
(V buffer) Cells were then centnfuged at 2000 x g for 20 minutes at 4^ Wet pellets were 
k weighed and frozen at -80^until ready for protein purification. 

A variety of methodologies known in the art can be used to purify the isolated 
25 proteins (Coligan et al, Current Protocols in Protein Science, John Wiley & Sons (1995)). 
For example, the fro/en cells can be thawed, resuspended in buffer and ruptured by several 
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passages through a small volume microfluidizer (Model M-l 10S, Microfluidics International 
Corp., Newton. MA). The resultant homogenate is centrifuged to yield a clear supernatant 
(crude extract) and, following filtration, the crude extract is fractioned over columns. 
Fractions are monitored by absorbar.ee at OD_, M , nm and peak fractions may be analyzed by 

5 SDS-PAGH. 

The concentrations of purified protein preparations are quantified 
spec.rophotometncally using absorbance coefficients calculated from amino acid content 
(Perkins, Eur. J. Biochem., 157:169-180 (1986)). Protein concentrations are also measured 
by the method of Bradford, Anal. Biochem., 72:248-254 (1976) and Lowry et al, J. Biol. 
10 Chem.. 193:265-275 (1951) using bovine serum albumin as a standard. 

SDS-polyacrylamide gels of various concentrations were purchased from BioRad 
(Hercules, CA), and stained with Coomassie blue. Molecular weight markers may include 
rabbit skeletal muscle myosin (200 kDa), E. coli (5-galactosidase (116 kDa), rabbit muscle 
phosphorylase B (97.4 kDa), bovine serum albumin (66.2 kDa), ovalbumin (45 kDa), bovine 
15 carbonic anyhdrase (31 kDa), soybean trypsin inhibitor (21 .5 kDa), egg white lysozyme (14.4 
kDa) and bovine aprotinin (6.5 kDa). 

Once a sufficient quantity of the desired protein has been obtained, it may be used for 
vanous purposes. A typical use is the production of antibodies specific for binding. These 
antibodies may be either polyclonal or monoclonal, and may be produced by in vitro or in 
20 vivo techniques well known in the art. Monoclonal antibodies to epitopes of any of the 
peptides identified and isolated as described can be prepared from murine hybndomas 
(Kohler, Nature, 256:495 (1975)). In summary, a mouse is inoculated with a few micrograms 
of HBM protein over a period of two weeks. The mouse is then sacrificed. The cells that 
produce antibodies are then removed from the mouse's spleen. The spleen cells are then 
25 fused with polyethylene glycol with mouse myeloma cells. The successfully fused cells are 
diluted in a microliter plate and growth of the culture is continued. The amount of antibody 
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per well is measured by immunoassay methods such as EL1SA (Engvall, Meth. Enzymoi, 
70:419 (1980)). Clones producing antibody can be expanded and further propagated to 
produce HBM antibodies. Other suitable techniques involve in vitro exposure of 
lymphocytes to the antigenic polypeptides, or alternatively, to selection of libraries of 
S antibodies in phage or s.milar vectors. See Huse et a>, Science. 246:1275-1281 (1989). For 
additional information on antibody production see Davis et al, Basic Methods in Molecular 
Biology. Elsevier, NY, Section 21-2 (1989). 
XVII. Methods of Use: Gene Therapy 

In recent years, significant technological advances have been made in the area of gene 
10 therapy for both genetic and acquired diseases. (Kay et al, Proc. Natl. Acad Sci. USA, 

94:12744-12746 (1997)) Gene therapy can be defined as the deliberate transfer of DNA for 
therapeutic purposes. Improvement in gene transfer methods has allowed for development of 
gene therapy protocols for the treatment of diverse types of diseases. Gene therapy has also 
taken advantage of recent advances in the identification of new therapeutic genes, 
15 improvement in both viral and nonviral gene delivery systems, better understanding of gene 
regulation, and improvement in cell isolation and transplantation. 

The preceding experiments identify the HBM gene as a dominant mutation conferring 
elevated bone mass. The fact that this mutation is dominant indicates that expression of the 
HBM protein causes elevated bone mass. Older individuals carrying the HBM gene, and, 
20 therefore expressing the HBM protein, do not suffer from osteoporosis. These individuals are 
equivalent to individuals being treated with the HBM protein. These observations are a 
strong experimental indication that therapeutic treatment with the HBM protein prevents 
osteoporosis. The bone mass elevating activity of the HBM gene is termed "HBM function." 
Therefore, according to the present invention, a method is also provided of supplying 
25 HBM function to mesenchymal stem cells (Onyia et al, J. Bone Miner. Res., 13:20-30(1998); 
Ko et al, Cancer Res., 56:4614-4619 (1996)). Supplying such a function provides protection 
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against osteoporosis. The HBM gene or a part of the gene may be introduced into the cell in 
a vector such that the gene remains extrachromosomal. In such a situation, the gene will be 
expressed by the cell from the extrachromosomal location. 

Vectors for introduction of genes both for recombination and for extrachromosomal 
maintenance are known in the art, and any suitable vector may be used. Methods for 
introducing DNA into cells such as electroporation, calcium phosphate co-precipitation, and 
viral transduction are known in the art, and the choice of method is within the competence of 
one skilled in the art (Robbins, Ed., Gene Therapy Protocols, Human Press, NJ (1997)). 
Cells transformed with the HBM gene can be used as model systems to study osteoporosis 
and drug treatments that promote bone growth. 

As generally discussed above, the HBM gene or fragment, where applicable, may be 
used in gene therapy methods in order to increase the amount of the expression products of 
such genes in mesenchymal stem cells. It may be useful also to increase the level of 
expression of a given HBM protein, or a fragment thereof, even in those cells in which the 
wild type gene is expressed normally. Gene therapy would be carried out according to 
generally accepted methods as described by, for example, Friedman, Therapy for Genetic 
Diseases, Friedman, Ed., Oxford University Press, pages 105-121 (1991). 

A virus or plasmid vector containing a copy of the HBM gene linked to expression 
control elements and capable of replicating inside mesenchymal stem cells, is prepared. 
Suitable vectors are known and described, for example, in U.S. Patent No. 5,252,479 and WO 
93/07282, the disclosures of which are incorporated by reference herein in their entirety. The 
vector is then injected into the patient, either locally into the bone marrow or systemically (in 
order to reach any mesenchymal stem cells located at other sites, i.e., in the blood). If the 
transfected gene is not permanently incorporated into the genome of each of the targeted 
cells, the treatment may have to be repeated periodically. 

Gene transfer systems known in the art may be useful in the practice of the gene 
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therapy methods of the present invention. These include viral and non-viral transfer methods. 
A number of viruses have been used as gene transfer vectors, including polyoma, i.e., SV40 
(Mad/ak et al,./. Gen. Virol., 73:1533-1536 (1992)), adenovirus (Berkner. Curr. Top. 
Microbiol. Immunol.. 158:39-61 (1992); Berkner et al. Bio Techniques, 6:6 1 6-629 ( 1 988); 
5 Gor/iglia et al, J. Virol.. 66:4407-4412 (1992); Quantin et al, Proc. Natl. Acad. Sci. USA, 
89:2581-2584 (1992); Rosenfeld et al. Cell, 68:143-155 (1992); Wilkinson et al, Nucl. Acids 
Res., 20:2233-2239 (1992); Stratford-Perricaudet et al, Hum. Gene Then, 1:241-256 (1990)), 
vaccinia virus (Mackett et al, Biotechnology, 24:495-499 (1992)), adeno-associated virus 
(Muzyczka, Curr. Top. Microbiol. Immunol., 158:91-123 (1992); Ohi et al, Gene, 89:279-282 
10 ( 1 990)), herpes viruses including HSV and EBV (Margolskee, Curr. Top. Microbiol. 

Immunol , 158:67-90 (1992); Johnson et al, J. Virol., 66:2952-2965 (1992); Fink et al, Hum. 
Gene Ther., 3:11-19 (1992); Breakfield et al, Mol. Neurobioi, 1 :337-371 (1987;) Fresse et al, 
Biochem. Pharmacol., 40:2189-2199 (1990)), and retroviruses of avian (Brandyopadhyay et 
al, Mol. Cell Biol., 4:749-754 (1984); Petropouplos et al, J. Virol., 66:3391-3397 (1992)), 
15 murine (Miller, Curr. Top. Microbiol. Immunol., 158:1-24 (1992); Miller et al, Mol. Cell 
Biol., 5:431-437 (1985); Sorge et al, Mol. Cell Biol., 4:1730-1737 (1984); Mann et al, J. 
Virol., 54:401-407 (1985)), and human origin (Page et al, J. Virol., 64:5370-5276 (1990); 
Buchschalcher et al, J. Virol, 66:273 1-2739 (1992)). Most human gene therapy protocols 
have been based on disabled murine retroviruses. 
20 Non-viral gene transfer methods known in the art include chemical techniques such as 

calcium phosphate coprecipitation (Graham et al, Virology, 52:456-467 (1973); Pellicer et al, 
Science, 209:1414-1422 (1980)), mechanical techniques, for example microinjection 
(Anderson et al, Proc. Natl. Acad. Sci. USA, 77:5399-5403 (1980); Gordon et al, Proc. Natl. 
Acad. Sci. USA, 77:7380-7384 (1980); Brinster et al, Cell, 27:223-231 (1981); Constantini et 
25 al. Nature. 294:92-94 (1981 )), membrane fusion-mediated transfer via liposomes (Feigner et 
al, Proc. Natl. Acad. Sci. USA, 84:7413-7417 (1987); Wang et al, Biochemistry, 28:9508- 
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9514 (1989); Kaneda et al, J. Biol. Chew., 264:12126-12129 (1989); Stewart et al, Hum. 
Gene The,:, 3:267-275 (1992); Nabel et al, Science, 249:1285-1288 (1990); Lim et al, 
Circulation. 83:2007-201 1 (1992)), and direct DNA uptake and receptor-mediated DNA 
transfer (Wolff et al. Science. 247:1465-1468 (1990); Wu et al, BioTeclmiques, 1 1 :474-485 
5 (1991 ); Zenke et al, Proc. Natl. Acad. Sci. USA, 87:3655-3659 (1990); Wu et al, J. Biol. 

Chen,., 264:16985-16987 (1989); Wolff et al, BioTeclmiques, 1 1:474-485 (1991); Wagner et 
al, 1990; Wagner et al, Proc. Natl. Acad. Sci. USA, 88:4255-4259 (1991 ); Cotten et al, Proc. 
Natl. Acad. Sci. USA, 87:4033-4037 (1990); Curiel et al, Proc. Natl. Acad. Sci. USA, 
88:8850-8854 (1991); Curiel et al. Hum. Gene The,:, 3:147-154(1991)). Viral-mediated 
10 gene transfer can be combined with direct in vivo vectors to the mesenchymal stem cells and 
not into the surrounding cells (Romano et al, In Vivo, 12(l):59-67 (1998); Gonez et al, Hum. 
Mol. Genetics, 7(12):1913-9 (1998)). Alternatively, the retroviral vector producer cell line 
can be injected into the bone marrow (Culver et al, Science, 256:1550-1552 (1992)). 
Injection of producer cells would then provide a continuous source of vector particles. This 
15 technique has been approved for use in humans with inoperable brain tumors. 

In an approach which combines biological and physical gene transfer methods, 
plasmid DNA of any size is combined with a poly ly sine-conjugated antibody specific to the 
adenovirus hexon protein, and the resulting complex is bound to an adenovirus vector. The 
trimolecular complex is then used to infect cells. The adenovirus vector permits efficient 
20 binding, internalization, and degradation of the endosome before the coupled DNA is 
damaged. 

Liposome/DNA complexes have been shown to be capable of mediating direct in vivo 
gene transfer. While in standard liposome preparations the gene transfer process is non- 
specific, localized in vivo uptake and expression have been reported in tumor deposits, for 
25 example, following direct in situ administration (Nabel, Hum. Gene Ther., 3:399-410 (1992)). 
XV III. Methods of Use: Transformed Hosts, Development of Pharmaceuticals and 
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Research Tools 

Cells and animals that carry the HBM gene can be used as model systems to study and 
test for substances that have potential as therapeutic agents (Onyia et al, ./. Bone Miner. Res., 
13:20-30 (1998); Broder et al. Bone, 2 1 :225-235 ( 1 997)). The cells are typically cultured 
5 mesenchymal stem cells. These may be isolated from individuals with somatic or germline 
HBM genes. Alternatively, the cell line can be engineered to carry the HBM gene, as 
described above. After a test substance is applied to the cells, the transformed phenotype of 
the cell is determined. Any trait of transformed cells can be assessed, including formation of 
bone matrix in culture (Broder et al. Bone, 21:225-235 (1997)), mechanical properties (Kizer 
10 et al, Proe. Natl. Acad. Sci. USA, 94:1013-1018 (1997)), and response to application of 
putative therapeutic agents. 

Animals for testing therapeutic agents can be selected after treatment of germline cells 
or zygotes. Such treatments include insertion of the Zmaxl gene, as well as insertion of the 
HBM gene and disrupted homologous genes. Alternatively, the inserted Zmaxl gene(s) 
15 and/or HBM gene(s) of the animals may be disrupted by insertion or deletion mutation of 
other genetic alterations using conventional techniques, such as those described by, for 
example, Capechi, Science, 244:1288 (1989); Valancuis et al, Mol. Cell Biol, 1 1:1402 
(1991); Hasty et al, Nature, 350:243 (1991); Shinkai et al, Cell, 68:855 (1992); Mombaerts et 
al. Cell, 68:869 (1992); Philpott et al, Science, 256:1448 (1992); Snouwaert et al, Science, 
20 257:1083 (1992); Donehower et al, Nature, 356:215 (1992). After test substances have been 
administered to the animals, the growth of bone must be assessed. If the test substance 
enhances the growth of bone, then the test substance is a candidate therapeutic agent. These 
animal models provide an extremely important vehicle for potential therapeutic products. 

Individuals carrying the HBM gene have elevated bone mass. The HBM gene causes 
25 this phenotype by altering the activities, levels, expression patterns, and modification states of 
other molecules involved in bone development. Using a variety of established techniques, it 
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is possible to identify molecules, preferably proteins or mRNAs, whose activities, levels, 
expression patterns, and modification states are different between systems containing the 
Zmax 1 gene and systems containing the HBM gene. Such systems can be, for example, cell- 
free extracts, cells, tissues or living organisms, such as mice or humans. For a mutant form of 
5 Zmaxl , a complete deletion of Zmax 1 . mutations lacking the extracellular or intracellular 
portion of the protein, or any other mutation in the Zmaxl gene may be used. It is also 
possible to use expression of antisense Zmaxl RNA or oligonucleotides to inhibit production 
of the Zmaxl protein. For a mutant form of HBM, a complete deletion of HBM, mutations 
lacking the extracellular or intracellular portion of the HBM protein, or any other mutation in 
10 the HBM gene may be used. It is also possible to use expression of antisense HBM RNA or 
oligonucleotides to inhibit production of the HBM protein. 

Molecules identified by comparison of Zmaxl systems and HBM systems can be used 
as surrogate markers in pharmaceutical development or in diagnosis of human or animal bone 
disease. Alternatively, such molecules may be used in treatment of bone disease. See, 
15 Schena et al, Science, 270:467-470 (1995). 

For example, a transgenic mouse carrying the HBM gene in the mouse homologue is 
constructed. A mouse of the genotype HBM/+ is viable, healthy and has elevated bone mass. 
To identify surrogate markers for elevated bone mass, HBM/+ (i.e., heterozygous) and 
isogenic +/+ (i.e., wild-type) mice are sacrificed. Bone tissue mRNA is extracted from each 
20 animal, and a "gene chip" corresponding to mRNAs expressed in the +/+ individual is 

constructed. mRNA from different tissues is isolated from animals of each genotype, reverse- 
transcribed, fluorescently labeled, and then hybridized to gene fragments affixed to a solid 
support. The ratio of fluorescent intensity between the two populations is indicative of the 
relative abundance of the specific mRNAs in the +/+ and HBM/+ animals. Genes encoding 
25 mRNAs over- and under-expressed relative to the wild-type control are candidates for genes 
coordinately regulated by the HBM gene. 
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One standard procedure for identification of new proteins that are part of the same 
signaling cascade as an already-discovered protein is as follows. Cells are treated with 
radioactive phosphorous, and the already-discovered protein is manipulated to be more ore 
less active. The phosphorylation state of other proteins in the cell is then monitored by 
5 polyacrylamide gel electrophoresis and autoradiography, or similar techniques. Levels of 
activity of the known protein may be manipulated by many methods, including, for example, 
comparing wild-type mutant proteins using specific inhibitors such as drugs or antibodies, 
simply adding or not adding a known extracellular protein, or using antisense inhibition of the 
expression of the known protein (Tamura et al, Science, 280(5369): 1 614-7 (1998); Meng, 
10 EMBOJ., 1 7(1 5):439 1-403 (1998); Cooper et al, Cell, 1:263-73 (1982)). 

In another example, proteins with different levels of phosphorylation are identified in 
TE85 osteosarcoma cells expressing either a sense or antisense cDNA for Zmaxl. TE85 cells 
normally express high levels of Zmaxl (Dong et al, Biochem. & Biophys. Res. Comm., 
251:784-790 (1998)). Cells containing the sense construct express even higher levels of 
15 Zmaxl , while cells expressing the antisense construct express lower levels. Cells are grown 
in the presence of 32 P, harvested, lysed, and the lysates run on SDS polyacrylamide gels to 
separate proteins, and the gels subjected to autoradiography (Ausubel et al, Current Protocols 
in Molecular Biology, John Wiley & Sons (1997)). Bands that differ in intensity between the 
sense and antisense cell lines represent phosphoproteins whose phosphorylation state or 
20 absolute level varies in response to levels of Zmaxl. As an alternative to the 32 P-labeling, 

unlabeled proteins may be separated by SDS-PAGE and subjected to immunoblotting, using 
the commercially available anti-phosphotyrosine antibody as a probe (Thomas et al, Nature, 
376(6537):267-71 (1995)). As an alternative to the expression of antisense RNA, transfection 
with chemically modified antisense oligonucleotides can be used (Woolf et al, Nucleic Acids 
25 Res., 18(7):1763-9 (1990)). 

Many bone disorders, such as osteoporosis, have a slow onset and a slow response to 
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treatment. It is therefore useful to develop surrogate markers for bone development and 
mmerali/ation. Such markers can be useful in developing treatments for bone disorders, and 
for diagnosing patients who may be at risk for later development of bone disorders. 
Examples of preferred markers are N- and C-lerminal telopeptidc markers described, for 
5 example, in U.S. Patent Nos. 5,455,179, 5,641,837 and 5,652,1 12, the disclosures of which 
are incorporated by reference herein in their entirety. In the area of HIV disease, CD4 counts 
and viral load are useful surrogate markers for disease progression (Vlahov et al, JAMA, 
279(1 ):35-40 (1998)). There is a need for analogous surrogate markers in the area of bone 
disease. 

10 A surrogate marker can be any characteristic that is easily tested and relatively 

insensitive to non-specific influences. For example, a surrogate marker can be a molecule 
such as a protein or mRNA in a tissue or in blood serum. Alternatively, a surrogate marker 
may be a diagnostic sign such as sensitivity to pain, a reflex response or the like. 

In yet another example, surrogate markers for elevated bone mass are identified using 
15 a pedigree of humans carrying the HBM gene. Blood samples are withdrawn from three 

individuals that carry the HBM gene, and from three closely related individuals that do not. 
Proteins in the serum from these individuals are electrophoresed on a two dimensional gel 
system, in which one dimension separates proteins by size, and another dimension separates 
proteins by isoelectric point (Epstein et al, Electrophoresis, 17(1 1):1655-70 (1996)). Spots 
20 corresponding to proteins are identified. A few spots are expected to be present in different 
amounts or in slightly different positions for the HBM individuals compared to their normal 
relatives. These spots correspond to proteins that are candidate surrogate markers. The 
identities of the proteins are determined by microsequencing, and antibodies to the proteins 
can be produced by standard methods for use in diagnostic testing procedures. Diagnostic 
25 assays for HBM proteins or other candidate surrogate markers include using antibodies 
described in this invention and a reporter molecule to detect HBM in human body fluids, 
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membranes, bones, cells, tissues or extracts thereof. The antibodies can be labeled by joining 
them covalently or noncovalently with a substance that provides a detectable signal. In many 
scientific and patent literature, a variety of reporter molecules or labels are described 
including radionuclides, enzymes, fluorescent, chemi-luminescent or chromogcnic agents 
5 (U.S. Patent Nos. 3,817,837; 3,850,752; 3.939,350; 3,996,345; 4,277,437; 4,275,149; and 
4.366,241). 

Using these antibodies, the levels of candidate surrogate markers are measured in 
normal individuals and in patients suffering from a bone disorder, such as osteoporosis, 
osteoporosis pseudoglioma, Engelmann's disease, Ribbing's disease, hyperphosphatasemia, 
10 Van Buchem's disease, melorheostosis, osteopetrosis, pychodysostosis, sclerosteosis, 

osteopoikilosis, acromegaly, Pagefs disease, fibrous dysplasia, tubular stenosis, osteogenesis 
imperfecta, hypoparathyroidism, pseudohypoparathyroidism, 

pseudopseudohypoparathyroidism, primary and secondary hyperparathyroidism and 
associated syndromes, hypercalciuria, medullary carcinoma of the thyroid gland, 
15 osteomalacia and other diseases. Techniques for measuring levels of protein in serum in a 

clinical setting using antibodies are well established. A protein that is consistently present in 
higher or lower levels in individuals carrying a particular disease or type of disease is a useful 
surrogate marker. 

A surrogate marker can be used in diagnosis of a bone disorder. For example, 
20 consider a child that present to a physician with a high frequency of bone fracture. The 

underlying cause may be child abuse, inappropriate behavior by the child, or a bone disorder. 
To rapidly test for a bone disorder, the levels of the surrogate marker protein are measured 
using the antibody described above. 

Levels of modification states of surrogate markers can be measured as indicators of 
25 the likely effectiveness of a drug that is being developed. It is especially convenient to use 
surrogate markers in creating treatments for bone disorders, because alterations in bone 
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development or mineralization may require a long time to be observed. For example, a set of 
bone mRNAs, termed the "HBM-inducible mRNA set" is found to be overexpressed in 
HBM/+ mice as compared to +/+ mice, as described above. Expression of this set can be 
used as a surrogate marker. Specifically, if treatment of +/+ mice with a compound results in 
5 overexpression of the HBM-inducible mRNA set, then that compound is considered a 
promising candidate for further development. 

This invention is particularly useful for screening compounds by using the Zmaxl or 
HBM protein or binding fragment thereof in any of a variety of drug screening techniques. 

The Zmaxl or HBM protein or fragment employed in such a test may either be free in 
10 solution, affixed to a solid support, or borne on a cell surface. One method of drug screening 
utilizes eukaryotic or prokaryotic host cells which are stably transformed with recombinant 
nucleic acids expressing the protein or fragment, preferably in competitive binding assays. 
Such cells, either in viable or fixed form, can be used for standard binding assays. One may 
measure, for example, for the formation of complexes between a Zmaxl or HBM protein or 
15 fragment and the agent being tested, or examine the degree to which the formation of a 

complex between a Zmaxl or HBM protein or fragment and a known ligand is interfered with 
by the agent being tested. 

Thus, the present invention provides methods of screening for drugs comprising 
contacting such an agent with a Zmaxl or HBM protein or fragment thereof and assaying (i) 
20 for the presence of a complex between the agent and the Zmax 1 or HBM protein or fragment, 
or (ii) for the presence of a complex between the Zmaxl or HBM protein or fragment and a 
ligand, by methods well known in the art. In such competitive binding assays the Zmaxl or 
HBM protein or fragment is typically labeled. Free Zmaxl or HBM protein or fragment is 
separated from that present in a proteimprotein complex, and the amount of free (i.e., 
25 uncomplexed) label is a measure of the binding of the agent being tested to Zmaxl or HBM 
or its interference with Zmaxl or HBM: ligand binding, respectively. 
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Another technique for drug screening provides high throughput screening for 
compounds having suitable binding affinity to the Zmaxl or HBM proteins and is described 
in detail in WO 84/03564. Briefly stated, large numbers of different small peptide test 
compounds are synthesized on a solid substrate, such as plastic pins or some other surface. 
5 The peptide test compounds are reacted with Zmaxl or HBM proteins and washed. Bound 
Zmaxl or HBM protein is then detected by methods well known in the art. Purified Zmaxl 
or HBM can be coated directly onto plates for use in the aforementioned drug screening 
techniques. However, non-neutralizing antibodies to the protein can be used to capture 
antibodies to immobilize the Zmaxl or HBM protein on the solid phase. 
10 This invention also contemplates the use of competitive drug screening assays in 

which neutralizing antibodies capable of specifically binding the Zmaxl or HBM protein 
compete with a test compound for binding to the Zmaxl or HBM protein or fragments 
thereof. In this manner, the antibodies can be used to detect the presence of any peptide that 
shares one or more antigenic determinants of the Zmaxl or HBM protein. 
15 A further technique for drug screening involves the use of host eukaryotic cell lines or 

cells (such as described above) that have a nonfunctional Zmaxl or HBM gene. These host 
cell lines or cells are defective at the Zmaxl or HBM protein level. The host cell lines or 
cells are grown in the presence of drug compound. The rate of growth of the host cells is 
measured to determine if the compound is capable of regulating the growth of Zmaxl or 

20 HBM defective cells. 

The goal of rational drug design is to produce structural analogs of biologically active 
proteins of interest or of small molecules with which they interact (e.g., agonists, antagonists, 
inhibitors) in order to fashion drugs which are, for example, more active or stable forms of 
the protein, or which, e.g., enhance or interfere with the function of a protein in vivo. See, 
25 e.g.. Hodgson, Bio/Technolog\\ 9:19-21 (1991 ). In one approach, one first determines the 
three-dimensional structure of a protein of interest (e.g., Zmaxl or HBM protein) or, for 
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example, of the Zmaxl- or HBM -receptor or ligand complex, by x-ray crystallography, by 
computer modeling or most typically, by a combination of approaches. Less often, useful 
information regarding the structure of a protein may be gained by modeling based on the 
structure of homologous proteins. An example of rational drug design is the development of 
5 HIV protease inhibitors (Erickson et al, Science, 249:527-533 (1990)). In addition, peptides 
(e.g., Zmaxl or HBM protein) are analyzed by an alanine scan (Wells, Methods in Enzymol., 
202:390-41 1 ( 1991 )). In this technique, an amino acid residue is replaced by Ala, and its 
effect on the peptide's activity is determined. Each of the amino acid residues of the peptide is 
analyzed in this manner to determine the important regions of the peptide. 
10 It is also possible to isolate a target-specific antibody, selected by a functional assay, 

and then to solve its crystal structure. In principle, this approach yields a pharmacore upon 
which subsequent drug design can be based. It is possible to bypass protein crystallography 
altogether by generating anti-idiotypic antibodies (anti-ids) to a functional, pharmacologically 
active antibody. As a mirror image of a mirror image, the binding site of the anti-ids would 
15 be expected to be an analog of the original receptor. The anti-id could then be used to 

identify and isolate peptides from banks of chemically or biologically produced banks of 
peptides. Selected peptides would then act as the pharmacore. 

Thus, one may design drugs which have, e.g., improved Zmaxl or HBM protein 
activity or stability or which act as inhibitors, agonists, antagonists, etc. of Zmaxl or HBM 
20 protein activity. By virtue of the availability of cloned Zmaxl or HBM sequences, sufficient 
amounts of the Zmaxl or HBM protein may be made available to perform such analytical 
studies as x-ray crystallography. In addition, the knowledge of the Zmaxl or HBM protein 
sequence provided herein will guide those employing computer modeling techniques in place 
of, or in addition to x-ray crystallography. 
25 XIX. Methods of Use: Avian and Mammalian Animal Husbandry 

The Zmaxl DNA and Zmaxl protein and/or the HBM DNA and HBM protein can be 
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used for vertebrate and preferably human therapeutic agents and for avian and mammalian 
veterinary agents, including for livestock breeding. Birds, including, for example, chickens, 
roosters, hens, turkeys, ostriches, ducks, pheasants and quails, can benefit from the 
identification of the gene and pathway for high bone mass. In many examples cited in 
5 literature (for example, McCoy et al, Res. Vet. ScL, 60(2): 1 85- 1 86 ( 1 996)), weakened bones 
due to husbandry conditions cause cage layer fatigue, osteoporosis and high mortality rates. 
Additional therapeutic agents to treat osteoporosis or other bone disorders in birds can have 
considerable beneficial effects on avian welfare and the economic conditions of the livestock 
industry, including, for example, meat and egg production. 
10 XX. Methods of use: Diagnostic assays using Zmaxl-specific oligonucleotides for 
detection of genetic alterations affecting bone development. 

In cases where an alteration or disease of bone development is suspected to involve an 
alteration of the Zmaxl gene or the HBM gene, specific oligonucleotides may be constructed 
and used to assess the level of Zmaxl mRNA or HBM mRNA, respectively, in bone tissue or 
15 in another tissue that affects bone development. 

For example, to test whether a person has the HBM gene, which affects bone density, 
polymerase chain reaction can be used. Two oligonucleotides are synthesized by standard 
methods or are obtained from a commercial supplier of custom-made oligonucleotides. The 
length and base composition are determined by standard criteria using the Oligo 4.0 primer 
20 Picking program (Wojchich Rychlik, 1992). One of the oligonucleotides is designed so that 
it will hybridize only to HBM DNA under the PCR conditions used. The other 
oligonucleotide is designed to hybridize a segment of Zmaxl genomic DNA such that 
amplification of DNA using these oligonucleotide primers produces a conveniently identified 
DNA fragment. For example, the pair of primers CCAAGTTCTGAGAAGTCC (SEQ ID 
25 NO:32) and AATACCTGAAACCA TACCTG (SEQ ID NO:33) will amplify a 530 base pair 
jf\ DNA fragment from a DNA sample when the following conditions are used: step 1 at 95/^for 
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120 seconds; step 2 at 95 j^1br 30 seconds; step 3 at 58^for 30 seconds; step 4 at 72^for 
120 seconds; where steps 2-4 are repeated 35 times. Tissue samples may be obtained from 
hair follicles, whole blood, or the buccal cavity. 

I he fragment generated by the above procedure is sequenced by standard techniques. 
Individuals heterozygous for the HBM gene will show an equal amount of G and T at the 
second position in the codon for glycine 171. Normal or homozygous wild-type individuals 
will show only G at this position. 

Other amplification techniques besides PCR may be used as alternatives, such as 
ligalion-mediated PCR or techniques involving Q-beta replicase (Cahill et al, Clin. Chem., 
37(9):1482-5 (1991)). For example, the oligonucleotides AGCTGCTCGTAGCT 
G TCTCTCCCTGGATCACGGGTACATGTACTGGACAGACTGGGT (SEQ ID NO:34) 
and TGAGACGCCCCGGATTGAGCGGGCAGGGATAGCTTATTCCCTGT 
GCCGCATTACGGC (SEQ ID NO:35) can be hybridized to a denatured human DNA 
sample, treated with a DNA ligase, and then subjected to PCR amplification using the primer 
oligonucleotides AGCTGCTCGTAG CTGTCTCTCCCTGGA (SEQ ID NO:36) and 
GCCGTAATGCGGCACAGGGAATAAGCT (SEQ ID NO:37). In the first two 
oligonucleotides, the outer 27 bases are random sequence corresponding to primer binding 
sites, and the inner 30 bases correspond to sequences in the Zmaxl gene. The T at the end of 
the first oligonucleotide corresponds to the HBM gene. The first two oligonucleotides are 
ligated only when hybridized to human DNA carrying the HBM gene, which results in the 
formation of an amplifiable 1 14 bp DNA fragment. 

Products of amplification can be detected by agarose gel electrophoresis, quantitative 
hybridization, or equivalent techniques for nucleic acid detection known to one skilled in the 
art of molecular biology (Sambrook et al, Molecular Cloning: A Laboratory Manual, Cold 
Spring Harbor Laboratory, Cold Spring, NY (1989)). 

Other alterations in the Zmaxl gene or the HBM gene may be diagnosed by the same 
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type of amplification-detection procedures, by using oligonucleotides designed to identify 
those alterations. These procedures can be used in animals as well as humans to identify 
alterations in Zmaxl or HBM that affect bone development. 

Expression of Zmaxl or HBM in bone tissue may be accomplished by fusing the 
5 cDNA of Zmaxlor HBM, respectively, to a bone-specific promoter in the context of a vector 
for genetically engineering vertebrate cells. DNA constructs are introduced into cells by 
packaging the DNA into virus capsids, by the use of cationic liposomes, electroporation, or 
by calcium phosphate transfection. Transfected cells, preferably osteoblasts, may be studied 
in culture or may be introduced into bone tissue in animals by direct injection into bone or by 
10 intravenous injection of osteoblasts, followed by incorporation into bone tissue (Ko et al, 
Cancer Research, 56(20):4614-9 (1996)). For example, the osteocalcin promoter, which is 
specifically active in osteoblasts, may be used to direct transcription of the Zmaxl gene or the 
HBM gene. Any of several vectors and transfection methods may be used, such as retroviral 
vectors, adenovirus vectors, or vectors that are maintained after transfection using cationic 
15 liposomes, or other methods and vectors described herein. 

Alteration of the level of functional Zmaxl protein or HBM protein affects the level 
of bone mineralization. By manipulating levels of functional Zmaxl protein or HBM protein, 
it is possible to affect bone development and to increase or decrease levels of bone 
mineralization. For example, it may be useful to increase bone mineralization in patients with 
20 osteoporosis. Alternatively, it may be useful to decrease bone mineralization in patients with 
osteopetrosis or Paget's disease. Alteration of Zmaxl levels or HBM levels can also be used 
as a research tool. Specifically, it is possible to identify proteins, mRNA and other molecules 
whose level or modification status is altered in response to changes in functional levels of 
Zmaxl or HBM. The pathology and pathogenesis of bone disorders is known and described, 
25 for example, in Rubin and Farber (Eds.), Pathology; 2nd Ed., S.B. Lippincott Co., 
Philadelphia, PA (1994). 
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A variety of techniques can be used to alter the levels of functional Zmaxl or HBM. 
For example, intravenous or intraosseous injection of the extracellular portion of Zmaxl or 
mutations thereof, or HBM or mutations thereof, will alter the level of Zmaxl activity or 
HBM activity, respectively, in the body of the treated human, animal or bird. Truncated 
5 versions of the Zmax 1 protein or HBM protein can also be injected to alter the levels of 
functional Zmaxl protein or HBM protein, respectively. Certain forms of Zmaxl or HBM 
enhance the activ ity of endogenous protein, while other forms are inhibitory. 

In a preferred embodiment, the HBM protein is used to treat osteoporosis. In a further 
preferred embodiment, the extracellular portion of the HBM protein is used. This HBM 
10 protein may be optionally modified by the addition of a moiety that causes the protein to 
adhere to the surface of cells. The protein is prepared in a pharmaceutical^ acceptable 
solution and is administered by injection or another method that achieves acceptable 
pharmacokinetics and distribution. 

In a second embodiment of this method, Zmaxl or HBM levels are increased or 
15 decreased by gene therapy techniques. To increase Zmaxl or HBM levels, osteoblasts or 

another useful cell type are genetically engineered to express high levels of Zmaxl or HBM 
as described above. Alternatively, to decrease Zmaxl or HBM levels, antisense constructs 
that specifically reduce the level of translatable Zmaxl or HBM mRNA can be used. In 
general, a tissue-nonspecific promoter may be used, such as the CMV promoter or another 
20 commercially available promoter found in expression vectors (Wu et al, Toxicol. Appl. 
Pharmacol., 141(l):330-9 (1996)). In a preferred embodiment, a Zmaxl cDNA or its 
antisense is transcribed by a bone-specific promoter, such as the osteocalcin or another 
promoter, to achieve specific expression in bone tissue. In this way, if a Zmaxl -expressing 
DNA construct or HBM-expressing construct is introduced into non-bone tissue, it will not be 
25 expressed. 

In a third embodiment of this method, antibodies against Zmaxl or HBM are used to 
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inhibit its function. Such antibodies are identified herein. 

In a fourth embodiment of this method, drugs that inhibit Zmaxl function or HBM 
function are used. Such drugs are described herein and optimized according to techniques of 
medicinal chemistry well known to one skilled in the art of pharmaceutical development. 
5 Zmaxl and HBM interact with several proteins, such as ApoE. Molecules that inhibit 

the interaction between Zmaxl or HBM and ApoE or another binding partner are expected to 
alter bone development and mineralization. Such inhibitors may be useful as drugs in the 
treatment of osteoporosis, osteopetrosis, or other diseases of bone mineralization. Such 
inhibitors may be low molecular weight compounds, proteins or other types of molecules. 

10 See, Kim et al, J. Biochem. (Tokyo), 124(6):1072-1076 (1998). 

Inhibitors of the interaction between Zmaxl or HBM and interacting proteins may be 
isolated by standard drug-screening techniques. For example, Zmaxl protein, (or a fragment 
thereof) or HBM protein (or a fragment thereof) can be immobilized on a solid support such 
as the base of microliter well. A second protein or protein fragment, such as ApoE is 

15 derivatized to aid in detection, for example with fluorescein. Iodine, or biotin, then added to 
the Zmaxl or HBM in the presence of candidate compounds that may specifically inhibit this 
protein-protein domain of Zmaxl or HBM, respectively, and thus avoid problems associated 
with its transmembrane segment. Drug screens of this type are well known to one skilled in 
the art of pharmaceutical development. 

20 Because Zmaxl and HBM are involved in bone development, proteins that bind to 

Zmaxl and HBM are also expected to be involved in bone development. Such binding 
proteins can be identified by standard methods, such as co-immunoprecipitation, co- 
fractionation, or the two-hybrid screen ( Ausubel et al, Current Protocols in Molecular 
Biology, John Wiley & Sons (1997)). For example, to identify Zmaxl -interacting proteins or 

25 HBM-interacting proteins using the two-hybrid system, the extracellular domain of Zmaxl or 
HBM is fused to LexA and expressed for the yeast vector pEG202 (the "bait") and expressed 
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in the yeast strain EGY48. The yeast strain is transformed with a "prey" library in the 
appropriate vector, which encodes a galactose-inducible transcription-activation sequence 
fused to candidate interacting proteins. The techniques for initially selecting and 
subsequently verifying interacting proteins by this method are well known to one skilled in 
5 the art of molecular biology ( Ausubel et al. Current Protocols in Molecular Biology, John 
Wiley & Sons (1997)). 

In a preferred embodiment, proteins that interact with HBM, but not Zmaxl, are 
identified using a variation of the above procedure (Xu et al, Proc. Natl. Acad. Sci. USA, 
94(23): 12473-8 (Nov. 1997)). This variation of the two-hybrid system uses two baits, and 
10 Zmaxl and HBM are each fused to Lex A and TetR, respectively. Alternatively, proteins that 
interact with the HBM but not Zmaxl are also isolated. These procedures are well known to 
one skilled in the art of molecular biology, and are a simple variation of standard two-hybrid 
procedures. 

As an alternative method of isolating Zmaxl or HBM interacting proteins, a 
15 biochemical approach is used. The Zmaxl protein or a fragment thereof, such as the 

extracellular domain, or the HBM protein or a fragment thereof, such as the extracellular 
domain, is chemically coupled to Sepharose beads. The Zmaxl - or HBM-coupled beads are 
poured into a column. An extract of proteins, such as serum proteins, proteins in the 
supernatant of a bone biopsy, or intracellular proteins from gently lysed TE85 osteoblastic 
20 cells, is added to the column. Non-specifically bound proteins are eluted, the column is 

washed several times with a low-salt buffer, and then tightly binding proteins are eluted with 
a high-salt buffer. These are candidate proteins that bind to Zmaxl or HBM, and can be 
tested for specific binding by standard tests and control experiments. Sepharose beads used 
for coupling proteins and the methods for performing the coupling are commercially available 
25 (Sigma), and the procedures described here are well known to one skilled in the art of protein 
biochemistry. 
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As a variation of the above procedure, proteins that are eluted by high salt from the 
Zmaxl- or HBM-Sepharose column are then added to an HBM-Zmaxl-sepharose column. 
Proteins that flow through without sticking are proteins that bind to Zmaxl but not to HBM. 
Alternatively, proteins that bind to the HBM protein and not to the Zmaxl protein can be 
5 isolated by reversing the order in which the columns are used. 

XXI. Method of Use: Transformation-Associated Recombination (TAR) Cloning 
Essential for the identification of novel allelic variants of Zmaxl is the ability to 
examine the sequence of both copies of the gene in an individual. To accomplish this, two 
"hooks," or regions of significant similarity, are identified within the genomic sequence such 
10 that they flank the portion of DNA that is to be cloned. Most preferably, the first of these 

hooks is derived from sequences 5' to the first exon of interest and the second is derived from 
sequences 3' to the last exon of interest. These two "hooks" are cloned into a bactenal/yeast 
shuttle vector such as that described by Larionov et al, Proc. Natl. Acad. Sci. USA, 94:7384- 
7387 (1997). Other similar vector systems may also be used. To recover the entire genomic 
15 copy of the Zmaxl gene, the plasmid containing the two "hooks" is linearized with a 

restriction endonuclease or is produced by another method such as PCR. This linear DNA 
fragment is introduced into yeast cells along with human genomic DNA. Typically, the yeast 
Saccharomyc.es cerevisiae is used as a host cell, although Larionov et al (in press) have 
reported using chicken host cells as well. During and after the process of transformation, the 
20 endogenous host cell converts the linear plasmid to a circle by a recombination event 

whereby the region of the human genomic DNA homologous to the "hooks" is inserted into 
the plasmid. This plasmid can be recovered and analyzed by methods well known to one 
skilled in the art. Obviously, the specificity for this reaction requires the host cell machinery 
to recognize sequences similar to the "hooks" present in the linear fragment. However, 100% 
25 sequence identity is not required, as shown by Kouprina et al, Genomics, 53(l):21-28 

(October 1998), where the author describes using degenerate repeated sequences common in 
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the human genome to recover fragments of human DNA from a rodent/human hybrid cell 
line. 

In another example, only one "hook" is required, as described by Larionov ct al, Proc. 
Sail. Acad. Sci. USA, 95(8):4469-74 (April 1998). For this type of experiment, termed 
5 "radial TAR cloning," the other region of sequence similarity to drive the recombination is 
derived from a repeated sequence from the genome. In this way, regions of DNA adjacent to 
the Zmaxl gene coding region can be recovered and examined for alterations that may affect 
function. 

XXII. Methods of Use: Genomic Screening 

10 The use of polymorphic genetic markers linked to the HBM gene or to Zmaxl is very 

useful in predicting susceptibility to osteoporosis or other bone diseases. Koller et al, Amer. 
J. Bone Min. Res., 13: 1903-1908 (1998) have demonstrated that the use of polymorphic 
genetic markers is useful for linkage analysis. Similarly, the identification of polymorphic 
genetic markers within the high bone mass gene will allow the identification of specific 
15 allelic variants that are in linkage disequilibrium with other genetic lesions that affect bone 
development. Using the DNA sequence from the BACs, a dinucleotide CAn repeat was 
identified and two unique PCR primers that will amplify the genomic DNA containing this 
repeat were designed, as shown below: 

B200E21C16_L: GAGAGGCTATATCCCTGGGC (SEQ ID NO:38) 
20 B200E21C16__R: ACAGCACGTGTTTAAAGGGG (SEQ ID NO:39) 

and used in the genetic mapping study. 

This method has been used successfully by others skilled in the art (e.g., Sheffield et 
al. Genet., 4:1837-1844 (1995); LeBlanc-Straceski et al. Genomics, 19:341-9 (1994); Chen et 
al. Genomics, 25:1-8 (1995)). Use of these reagents with populations or individuals will 
25 predict their risk for osteoporosis. Similarly, single nucleotide polymorphisms (SNPs). such 
as those shown in Table 4 above, can be used as well to predict risk for developing bone 
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diseases or resistance to osteoporosis in the case of the HBM gene. 
XXIII. Methods of Use: Modulators of Tissue Calcification 

The calcification of tissues in the human body is well documented. Towler et al, ./. 
Biol. Clwm., 273:30427-34 (1998) demonstrated that several proteins known to regulate 
5 calcification of the developing skull in a model system are expressed in calcified aorta. The 
expression of Msx2. a gene transcribed in osteoprogenitor cells, in calcified vascular tissue 
indicates that genes which are important in bone development are involved in calcification of 
other tissues. Treatment with HBM protein, agonists or antagonists is likely to ameliorate 
calcification (such as the vasculature, dentin and bone of the skull visera) due to its 
10 demonstrated effect on bone mineral density. In experimental systems where tissue 

calcification is demonstrated, the over-expression or repression of Zmaxl activity permits the 
identification of molecules that are directly regulated by the Zmaxl gene. These genes are 
potential targets for therapeutics aimed at modulating tissue calcification. For example, an 
animal, such as the LDLR -/-, mouse is fed a high fat diet and is observed to demonstrate 
15 expression of markers of tissue calcification, including Zmaxl . These animals are then 

treated with antibodies to Zmaxl or HBM protein, antisense oligonucleotides directed against 
Zmaxl or HBM cDNA, or with compounds known to bind the Zmaxl or HBM protein or its 
binding partner or ligand. RNA or proteins are extracted from the vascular tissue and the 
relative expression levels of the genes expressed in the tissue are determined by methods well 
20 known in the art. Genes that are regulated in the tissue are potential therapeutic targets for 
pharmaceutical development as modulators of tissue calcification. 

The nucleic acids, proteins, peptides, amino acids, small molecules or other 
pharmaceutically useful compounds of the present invention that are to be given to an 
individual may be administered in the form of a composition with a pharmaceutically 
25 acceptable carrier, excipient or diluent, which are well known in the art. The individual may 
be a mammal or a bird, preferably a human, a rat, a mouse or bird. Such compositions may 
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be administered to an individual in a pharmaceutical^' effective amount. The amount 
administered will vary depending on the condition being treated and the patient being treated. 
The compositions may be administered alone or in combination with other treatments. 

EXAMPLES 

5 The present invention is described by reference to the following Examples, which are 

offered by way of illustration and are not intended to limit the invention in any manner. 
Standard techniques well known in the art or the techniques specifically described below 
were utilized. 
Example 1 

10 The propositus was referred by her physicians to the Creighton Osteoporosis Center 

for evaluation of what appeared to be unusually dense bones. She was 1 8 years old and came 
to medical attention two years previous because of back pain, which was precipitated by an 
auto accident in which the car in which she was riding as a passenger was struck from behind. 
Her only injury was soft tissue injury to her lower back that was manifested by pain and 
15 muscle tenderness. There was no evidence of fracture or subluxation on radiographs. The 
pain lasted for two years, although she was able to attend school full time. By the time she 
was seen in the Center, the pain was nearly resolved and she was back to her usual activities 
as a high school student. Physical exam revealed a normal healthy young woman standing 66 
inches and weighing 128 pounds. Radiographs of the entire skeleton revealed dense looking 
20 bones with thick cortices. All bones of the skeleton were involved. Most importantly, the 

shapes of all the bones were entirely normal. The spinal BMC was 94.48 grams in Ll-4, and 
the spinal BMD was 1 .667 gm/cm 2 in Ll-4. BMD was 5.62 standard deviations (SD) above 
peak skeletal mass for women. These were measured by DXA using a Hologic 2000-. Her 
mother was then scanned and a lumbar spinal BMC of 58.05 grams and BMD of 1.500 
25 gm/cm 2 were found. Her mother's values place her 4. 1 2 SD above peak mass and 4.98 SD 
above her peers. Her mother was 5 1 years old, stood 65 inches and weighed 140 pounds. 
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Her mother w as in excellent health with no history of musculoskeletal or other symptoms. 
Her father's lumbar BMC was 75.33 grams and his BMD was 1.118 gm/enr. These values 
place him 0.25 SD above peak bone mass for males. He was in good health, stood 72 inches 
tall, and weighed 187 pounds. 
5 These clinical data suggested that the propositus inherited a trait from her mother, 

w hich resulted in very high bone mass, but an otherwise normal skeleton, and attention was 
focused on the maternal kindred. In U.S. Patent No. 5,691,153, twenty- two of these 
members had measurement of bone mass by DXA. In one case, the maternal grandfather of 
the propositus, was deceased, however, medical records, antemortem skeletal radiographs and 
10 a gall bladder specimen embedded in paraffin for DNA genotyping were obtained. His 

radiographs showed obvious extreme density of all of the bones available for examination 
including the femur and the spine, and he was included among the affected members. In this 
invention, the pedigree has been expanded to include 37 informative individuals. These 
additions are a significant improvement over the original kinship (Johnson et al, Am. J. Hum. 
15 Genet., 60:1326-1332 (1997)) because, among the fourteen individuals added since the 

original study, two individuals harbor key crossovers. X-linkage is ruled out by the presence 
of male-to-male transmission from individual 12 to 14 and 15. 
Example 2 

The present invention describes DNA sequences derived from two BAC clones from 
20 the HBM gene region, as evident in Table 6 below, which is an assembly of these clones. 

Clone b200e21-h (ATCC No. 980812; SEQ ID NOS: 10-1 1) was deposited at the American 
Type Culture Collection (ATCC), 10801 University Blvd., Manassas, VA 201 10-2209 
U.S.A., on December 30, 1997. Clone b527dl2-h (ATCC No. 980720; SEQ ID NOS: 5-9) 
was deposited at the American Type Culture Collection (ATCC), 10801 University Blvd., 
25 Manassas, VA 201 10-2209 U.S.A., on October 2, 1998. These sequences are unique reagents 
that can be used by one skilled in the art to identify DNA probes for the Zmaxl gene, PCR 
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primers to amplify the gene, nucleotide polymorphisms in the Zmaxl gene, or regulatory 
elements of the Zmaxl gene. 

TABLE 6 



C ontig 


ATCC No. 


SEQ ID 
NO. 


Length (base 
pairs) 


b527dl2-h_contig302G 


980720 


5 


3096 


b527dl2-h_contig306G 


980720 


6 


26928 


b527dl2-h_contig307G 


980720 


7 


29430 


b527dl2-h_contig308G 


980720 


8 


33769 


b527dl2-h_contig309G 


980720 


9 


72049 


b200e21-h_contigl 


980812 


10 


8705 


b200e2 1 -h_contig4 


980812 


11 


66933 



The disclosure of each of the patents, patent applications and publications cited in the 
specification is hereby incorporated by reference herein in its entirety. 

Although the invention has been set forth in detail, one skilled in the art will recognize 
that numerous changes and modifications can be made, and that such changes and 
modifications may be made without departing from the spirit and scope of the invention. 
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